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1. FOREWORD 1 

Among the many physics processes at TeV hadron colliders, we look most eagerly for those that display 
signs of the Higgs boson or of new physics. We do so however amid an abundance of processes that 
proceed via Standard Model (SM) and in particular Quantum Chromodynamics (QCD) interactions, 
and that are interesting in their own right. Good knowledge of these processes is required to help us 
distinguish the new from the known. Their theoretical and experimental study teaches us at the same 
time more about QCD/SM dynamics, and thereby enables us to further improve such distinctions. This is 
important because it is becoming increasingly clear that the success of finding and exploring Higgs boson 
physics or other New Physics at the Tevatron and LHC will depend significantly on precise understanding 
of QCD/SM effects for many observables. 

To improve predictions and deepen the study of QCD/SM signals and backgrounds was therefore 
the ambition for our QCD/SM working group at this Les Houches workshop. Members of the working 
group made significant progress towards this on a number of fronts. A variety of tools were further devel- 
oped, from methods to perform higher order perturbative calculations or various types of resummation, 
to improvements in the modelUng of underlying events and parton showers. Furthermore, various precise 
studies of important specific processes were conducted. 

A signficant part of the activities in Les Houches revolved around Monte Carlo simulation of 
collision events. A number of contributions in this report reflect the progress made in this area. At present 
a large number of Monte Carlo programs exist, each written with a different purpose and employing 
different techniques. Discussions in Les Houches revealed the need for an accessible primer on Monte 
Carlo programs, featuring a listing of various codes, each with a short description, but also providing a 
low-level explanation of the underlying methods. This primer has now been compiled and a synopsis of 
it is included here as the first contribution to this report (see below for where to obtain the full document). 

This report reflects the hard and creative work by the many contributors which took place in the 
working group. After the MC guide description, the next contributions report on progress in describing 
multiple interactions, important for the LHC, and underlying events. An announcement of a Monte 
Carlo database, under construction, is followed by a number of contributions improving parton shower 
descriptions. Subsequently, a large number of contributions address resummations in various forms, after 
which follow studies of QCD effects in pion pair, top quark pair and photon pair plus jet production. After 
a study of electroweak corrections to hadronic precision observables, the report ends by presenting recent 
progress in methods to compute finite order corrections at one-loop with many legs, and at two-loop. 

2. LES HOUCHES GUIDEBOOK TO MONTE CARLO GENERATORS FOR HADRON COL- 
LIDER PHYSICS 

Editors: M. Dobbs, S. Frixione, E. Laenen, K. Tollefson 

Contributing Authors: H. Baer, E. Boos, B. Cox, M. Dobbs, R. Engel, S. Frixione, W. Giele, J. Huston, 
S. Ilyin, B. Kersevan, F. Krauss, Y. Kurihara, E. Laenen, L. Ldnnblad, F. Maltoni, M. Mangano, S. Odaka, 
P. Richardson, A. Ryd, T. Sjdstrand, P. Skands, Z. Was, B.R. Webber, D. Zeppenfeld 
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Abstract 

Recently the collider physics community has seen significant advances in the 
formaUsms and implementations of event generators. This review is a primer 
of the methods commonly used for the simulation of high energy physics 
events at particle colliders. We provide brief descriptions, references, and links 
to the specific computer codes which implement the methods. The aim is to 
provide an overview of the available tools, allowing the reader to ascertain 
which tool is best for a particular application, but also making clear the limita- 
tions of each tool. 

Due to its long length and stand-alone nature, the Monte Carlo Guidebook en- 
try in the Les Houches proceedings has been pubUshed as a separate document 
(hep-ph / 0403045). The table of contents follows. 

1. Introduction 

2. The Simulation of Hard Processes 

3. Tree Level Matrix Element Generators 

4. Higher Order Corrections - Perturbative QCD Computations 

5. Parton Distribution Functions 

6. Higher Order Corrections - Showering and Hadronization Event Generators 

7. Resummation 

8. Combining Matrix Elements with Showering 

9. Conclusions 

10. Acknowledgments 

3. MULTIPLE INTERACTIONS AND BEAM REMNANTS ^ 
3.1 Introduction 

Hadrons are composite systems of quarks and gluons. A direct consequence is the possibility to have 
hadron-hadron collisions in which several distinct pairs of partons collide with each other, i.e. multiple 
interactions, a.k.a. multiple scatterings. At first glance, the divergence of the perturbative ^-channel one- 
gluon-exchange graphs in the p± ^ limit implies an infinity of interactions per event. However, the 
perturbative framework does not take into account screening from the fact that a hadron is in an overall 
colour singlet state. Therefore an effective cutoff p_Lmin of the order of one to a few GeV is introduced, 
representing an inverse colour correlation distance inside the hadron. For realistic p±ram values most 
inelastic events in high-energy hadronic collisions should then contain several perturbatively calculable 
interactions, in addition to whatever nonperturbative phenomena may be present. 

Although most of this activity is not hard enough to play a significant role in the description of 
high-px jet physics, it can be responsible for a large fraction of the total multiplicity (and large fluctu- 
ations in it) for semi-hard (mini-)jets in the event, for the details of jet profiles and for the jet pedestal 
effect, leading to random as well as systematic shifts in the jet energy scale. Thus, a good understanding 
of multiple interactions would seem prerequisite to carrying out precision studies involving jets and/or 
the underlying event in hadronic collisions. 

^Contributed by: T. Sjostrand and P. Skands 
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In an earlier study [1], it was argued that all the underlying event activity is triggered by the 
multiple interactions mechanism. However, while the origin of underlying events is thus assumed to be 
perturbative, many nonperturbative aspects still need to be considered and understood: 

( i) What is the detailed mechanism and functional form of the dampening of the perturbative cross section 
at small p±1 (Certainly a smooth dampening is more realistic than a sharp p±inm cutoff.) 

(ii) Which energy dependence would this mechanism have? 

(Hi) How is the internal structure of the proton reflected in an impact-parameter-dependent multiple 
interactions rate, as manifested e.g. in jet pedestal effects? 

(iv) How can the set of colliding partons from a hadron be described in terms of correlated multiparton 
distribution functions of flavours and longitudinal momenta? 

(v) How does a set of initial partons at some low perturbative cutoff scale, 'initiators', evolve into such a 
set of colliding partons? (Two colliding partons could well have a common initiator.) Is standard DGLAP 
evolution sufficient, or must BFKL/CCFM effects be taken into account? 

(vi) How would the set of initiators correlate with the flavour content of, and the longitudinal momentum 
sharing inside, the left-behind beam remnant? 

(vii) How are the initiator and remnant partons correlated by confinement effects ('primordial k±')7 

(viii) How are all produced partons, both the interacting and the beam-remnant ones, correlated in colour? 
Is the large number-of-colours limit relevant, wherein partons can be hooked up into strings (with quarks 
as endpoints and gluons as intermediate kinks) representing a linear confinement force [2]? 

(ix) How is the original baryon number of an incoming proton reflected in the colour topology? 

(x) To what extent would a framework with independently fragmenting string systems, as defined from 
the colour topology, be modified by the space-time overlap of several strings? 

Needless to say, we should not expect to find a perfect solution to any of these issues, but only 
successively improved approximations. The framework in [1] is very primitive in a number of respects. 
Nevertheless, it has turned out to be quite successful. Thus the Pythia Tune A of R.D. Field [3] is 
capable of describing a host of jet and minimum-bias event data at the Tevatron. The model appears 
inadequate to fully describe correlations and fluctuations, however, and we would expect a poor perfor- 
mance for several topics not yet studied experimentally. 

In particular, only very simple beam remnant structures could technically be dealt with in [1]. One 
recent development was the extension of the standard Lund string framework [2] to include a junction 
fragmentation description [4] that allows the hadronization of nontrivial colour topologies containing 
non-zero baryon number. In the context of multiple interactions, this improvement means that almost 
arbitrarily complicated baryon beam remnants may now be dealt with, hence many of the restrictions 
present in the old model are no longer necessary. 

Here, we report on the development of a new model for the flavour-, colour-, and momentum- 
correlated partonic structure involved in a hadron-hadron collision, i.e. partly addressing several of the 
points above. We first present the main work on flavour and momentum space correlations, and thereafter 
separately the very thorny issue of colour correlations, before concluding. A more complete description 
of the model, also including references to experimental data and other theoretical ideas, and with com- 
ments on all the issues, may be found in [5]. A toy model study of the first two points is found in [6]. 
The Pythia manual [7] contains some complementary information. 

3.2 Correlated Parton Densities 

Consider a hadron undergoing multiple interactions in a collision. Such an object should be described 
by multi-parton densities, giving the joint probability of simultaneously finding n partons with flavours 
/i) • • • ) /n> carrying momentum fractions xi, . . . , x„ inside the hadron, when probed by interactions at 
scales Qi, . . . ,Q^- However, we are nowhere near having sufficient experimental information to pin 
down such distributions. Therefore, and wishing to make maximal use of the information that we do 
have, namely the standard one-parton-inclusive parton densities, we propose the following strategy. 
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As described in [1], the interactions may be generated in an ordered sequence of falling p±. For 
the hardest interaction, all smaller pj_ scales may be effectively integrated out of the (unknown) fully 
correlated distributions, leaving an object described by the standard one-parton distributions, by defini- 
tion. For the second and subsequent interactions, again all lower-p^ scales can be integrated out, but the 
correlations with the first cannot, and so on. Thus, we introduce modified parton densities, that correlate 
the z'th interaction and its shower evolution to what happened in the z — 1 previous ones. 

The first and most trivial observation is that each interaction i removes a momentum fraction Xi 
from the hadron remnant. Already in [1] this momentum loss was taken into account by assuming a 
simple scaling ansatz for the parton distributions, f{x) f{x/X)/X, where X = I — Y17=i 
momentum remaining in the beam hadron after the n first interactions. Effectively, the PDF's are simply 
'squeezed' into the range x G [0, X]. 

Next, for a given baryon, the valence distribution of flavour / after n interactions, qfvn{x, Q^), 
should integrate to the number Nf^^ of valence quarks of flavour / remaining in the hadron remnant. 
This rule may be enforced by scahng the original distribution down, by the ratio of remaining to original 
valence quarks Nfy^/Nf^Q, in addition to the x scahng mentioned above. 

Also, when a sea quark is knocked out of a hadron, it must leave behind a corresponding antisea 
parton in the beam remnant. We call this a companion quark. In the perturbative approximation the sea 
quark qs and its companion qc come from a gluon branching g ^ qs + qc (it is implicit that if qs is a 
quark, qc is its antiquark). Starting from this perturbative ansatz, and neglecting other interactions and 
any subsequent perturbative evolution of the qc, we obtain the Qc distribution from the probabihty that a 
sea quark qg, carrying a momentum fraction Xg, is produced by the branching of a gluon with momentum 
fraction y, so that the companion has a momentum fraction x = y — Xg, 

qo{x; Xs) oc / g{y) Pg^q.q, (z) 5{xs - zy) dz = ^^J'^^^ -Pg-q^qc f ' (1) 

Jo Xg -r X \Xs -\- X y 

with Pg^q^q^ the usual DGLAP gluon sphtting kernel. A simple ansatz ^(a;) oc (1 — x)"'/x is here used 
for the gluon. Normalizations are fixed so that a sea quark has exactly one companion. Quahtatively, 

xqc{x; Xs) is peaked around x « Xg, by virtue of the symmetric Pg^q^q^ splitting kernel. 

Without any further change, the reduction of the valence distributions and the introduction of com- 
panion distributions, in the maimer described above, would result in a violation of the total momentum 
sum rule, that the x-weighted parton densities should integrate to X: by removing a valence quark from 
the parton distributions we also remove a total amount of momentum corresponding to {xfy), the average 
momentum fraction carried by a valence quark of flavour /, 

_ xqfyn{x,Q^) dx , > 

Jo y/vnv-^'^ ; ux 

and by adding a companion distribution we add an analogously defined momentum fraction. 

To ensure that the momentum sum rule is still respected, we assume that the sea-i-gluon normal- 
izations fluctuate up when a valence distribution is reduced and down when a companion distribution is 
added, by a multiplicative factor 

l-E/^/vo(a;/vo) 

The requirement of a physical x range is of course still maintained by 'squeezing' all distributions into 
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^ remnant parton Q 




Fig. 1: (a) The initial state of a baryon, with the valence quarks colour-connected via a central string junction J. (b) Example of 
a topology with initiators connected at random, (c) Alternative with the junction in the remnant. 



the interval x G [0, X]. The full parton distributions after n interactions thus take the forms 



1 

X 



N 



fvO 



9n{x) = —go [^,Q ) , 



iV/vn / X 

QfvO 



Q^)+aqf,o{^,Q^)+Y,qf. 



X' 



- 1 



(4) 
(5) 



where gjvO (9/so) denotes the original valence (sea) distribution of flavour /, and the index j on the 
companion distributions qjc^ counts different companion quarks of the same flavour /. 

After the perturbative interactions have each taken their fraction of longitudinal momentum, the 
remaining momentum is to be shared between the beam remnant partons. Here, valence quarks receive 
an X picked at random according to a small-Q^ valence-like parton density, while sea quarks must be 
companions of one of the initiator quarks, and hence should have an x picked according to the qdx; Xg) 
distribution introduced above. In the rare case that no valence quarks remain and no sea quarks need 
be added for flavour conservation, the beam remnant is represented by a gluon, carrying all of the beam 
remnant longitudinal momentum. 

Further aspects of the model include the possible formation of composite objects in the beam 
remnants (e.g. diquarks) and the addition of non-zero primordial k± values to the parton shower initiators. 
Especially the latter introduces some complications, to obtain consistent kinematics. Details on these 
aspects are presented in [5]. 



3.3 Colour Correlations 

The initial state of a baryon may be represented by three valence quarks, connected antisymmetrically 
in colour via a central junction, which acts as a switchyard for the colour flow and carries the net baryon 
number. Fig. la. 

The colour-space evolution of this state into the initiator and remnant partons actually found in 
a given event is not predicted by perturbation theory, but is crucial in determining how the system 
hadronizes; in the Lund string model [2], two colour-connected final state partons together define a 
string piece, which hadronizes by successive non-perturbative breakups along the string. Thus, the 
colour flow of an event determines the topology of the hadronizing strings, and consequently where 
and how many hadrons will be produced. The question can essentially be reduced to one of choosing 
a fictitious sequence of gluon emissions off the initial valence topology, since sea quarks together with 
their companion partners are associated with parent gluons, by construction. 
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The simplest solution is to assume that gluons are attached to the initial quark lines in a random 
order, see Fig. lb. If so, the junction would rarely be colour-connected directly to two valence quarks in 
the beam remnant, and the initial-state baryon number would be able to migrate to large p± and small 
xp values. While such a mechanism should be present, there are reasons to believe that a purely random 
attachment exaggerates the migration effects. Hence a free parameter is introduced to suppress gluon 
attachments onto colour lines that lie entirely within the remnant, so that topologies such as Fig. Ic 
become more likely. 

This still does not determine the order in which gluons are attached to the colour line between 
a valence quark and the junction. We consider a few different possibilities: 1) random, 2) gluons are 
ordered according to the rapidity of the hard scattering subsystem they are associated with, and 3) gluons 
are ordered so as to give rise to the smallest possible total string lengths in the final state. The two latter 
possibilities correspond to a tendency of nature to minimize the total potential energy of the system, 
i.e. the string length. Empirically such a tendency among the strings formed by multiple interactions is 
supported e.g. by the observed rapid increase of {p±) with ncharged- It appears, however, that a string 
minimization in the initial state is not enough, and that also the colours inside the initial-state cascades 
and hard interactions may be nontrivially correlated. These studies are still ongoing, and represent the 
major open issues in the new model. 

3.4 Conclusion 

A new model for the underlying event in hadron-hadron collisions [5] has been introduced. This model 
extends the multiple interactions mechanism proposed in [1] with the possibility of non-trivial flavour 

and momentum correlations, with initial- and final-state showers for all interactions, and with several 
options for colour correlations between initiator and remnant partons. Many of these improvements rely 
on the development of junction fragmentation in [4]. 

This is not the end of the line. Rather we see that many issues remain to understand better, such 
as colour correlations between partons in interactions and beam remnants, whereas others have not yet 
been studied seriously, such as the extent to which two interacting partons stem from the same initiator. 
Theoretical advances alone cannot solve all problems; guidance will have to come from experimental 
information. The increased interest in such studies bodes well for the future. 

4. DESCRIBING MINIMUM BIAS AND THE UNDERLYING EVENT AT THE LHC IN 
PYTHIA AND PHOJET ^ 

4.1 Introduction 

Our ability to describe parton scatterings through QCD depends on the amount of transverse momenta 
with respect to the collision axis (pt) involved in a given scattering [8]. QCD has been fairly successful 
in describing quark, anti-quark and gluon scatterings involving large amounts of transverse momenta 
(pt >> Aqcd)» also known as "hard" interactions. On the other hand, QCD simply caimot be applied to 
interactions with small transverse momenta (or "soft" interactions) because the strong coupling constant, 
"s (Q^)' becomes too large for perturbation theory to be applied and QCD models suffer from diver- 
gent cross sections as p^ — > [8]. Most high-energy hadron collisions are dominated by soft partonic 
interactions. 

A full picture of high-energy hadron collisions will typically combine perturbative QCD to explain 
parton interactions where it is applicable (high-p^ scatterings), with an alternative phenomenological 
approach to describe soft processes. Examples of these are the Dual Parton Model (DPM) [9] and 
modified versions of QCD in which the divergencies presented by the running coupUng constant are 
phenomenologically corrected to reproduce experimental observations [1]. 

^Contributed by: A. Moraes, C. Buttar, and I. Dawson 
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In this article we investigate two Monte Carlo (MC) event generators, PYTHIA6.214 [10, 11] 
and PHOJET1.12 [12, 13], focusing on their models for soft interactions in hadron-hadron collisions. 
Aiming to check the consistency of these models, we compare their predictions to wide range of data 
for minimum bias and the underlying event. A tuning for PYTH1A6.214 is presented and examples of 
its predictions are compared to those generated with PH0JET1.12 . Predictions for levels of particle 
production and event activity at the LHC for interactions dominated by soft processes such as minimum 
bias interactions and the underlying event associated to jet production are also discussed. 

4.2 PYTHIA Model for Hadron Collisions 

A comprehensive description of PYTHIA can be found at [11] and references therein. The evolution 
of a hadronic event generated by PYTHIA is based on parton-parton scatterings [10, 11]. In this model 
the total rate of parton interactions, Nparton-parton^ as a function of the transverse momentum scale 
pt, is assumed to be given by perturbative QCD. At reasonably large p^ values (pi > 2 GeV) parton 
scatterings can be correctly described by the standard perturbative QCD, but to extend the parton-parton 
scattering framework to the low-pt region a regularisation to correct the divergence in the cross-section 
is introduced. 

In order to deal with low-pt interactions, PYHTIA introduces a cut-off parameter pt^j„ given by 



which can be interpreted as the inverse of some colour screening length in the hadron [6]. There are two 
strategies, or scenarios, to implement the cut-off parameter defined by equation 6. 

In the first one, labelled "simple" scenario, an effective cut-off is established at p^ , which means 



that da /dpi = for p^ < p^^^^. This model assumes that different pairwise interactions take place 
essentially independent of each other, and that therefore the number of interactions in an event is given 
by a Poissonian distribution [1]. In the second approach, called the 'complex' scenario, the probability 
associated with each interacting parton depends on the assumed matter distribution inside the colliding 
hadrons. In the 'complex' scenario an impact parameter dependent approach is therefore introduced [1]. 

The parameters defining pf^^^ are PARP(81), PARP(82), PARP(89) and PARP(90). The factor 
1.9 GeV is defined in the simple scenario by PARP(81) and by PARP(82) in the complex scenario. The 
energy scale 1 TeV is defined by PARP(89) and is included in equation (1) to be a convenient tuning 
parameter rather than a parameter with physical meaning. PARP(90) gives the power with which ^tmin 
varies with the centre of mass energy, -y/i. The default option is set as PARP(90)=0.16 [10, 11]. 

4.3 PHOJET 

The physics model used in the MC event generator PHOJET combines the ideas of the DPM [9] with 
perturbative QCD [8] to give an almost complete picture of high-energy hadron collisions [12-14]. 

PHOJET is formulated as a two-component model containing contributions from both soft and 
hard interactions. The DPM is used describe the dominant soft processes and perturbative QCD is applied 

to generate hard interactions [13]. 

The model employed by PHOJET is based on the calculation of scattering amplitudes, taking into 
account the unitarization principle. Comparisons between the calculated results for cross-sections and 
the available data are used to determine the unknown model parameters (couplings, Pomeron intercepts 
and slope parameters), which are needed to generate multiparticle final states produced in inelastic inter- 
actions [12, 13]. 

The soft, (Jsofu and hard, ahard^ cross sections are inclusive cross sections and the average multi- 




(6) 
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plicities of soft and hard scatterings in an inelastic event are 



(us) = — —, {rih) = , (7) 

respectively. The hard scatterings are mostly independent of each other, being related only by the sharing 
of energy and momentum of the incoming protons. These multiplicities increase with the colliding 
centre-of-mass energy. For pp collisions at y/s = 14 TeV a considerable part of interactions is expected 
to have more than one hard or soft scattering. 

4.4 Minimum Bias Interactions 

Throughout this article, we will associate minimum bias events with non-single diffractive inelas- 
tic (NSD) interactions, following the experimental definition used in [15-19]. In the language of 
the MC event generators used in this work, this means that subprocesses 94 and 95 are switched 
on in PYTHIA6.214 (MSUB(94)=1 and MSUB(95)=1), and processes IPR0N(1,1), IPR0N(4,1) and 
IPR0N(7,1) in PH0JET1.12. For both generators, we also adapt the MC distributions to the data by 
setting TT^ , Kg and A*' as stable particles. 

4.41 KNO Distribution 

The KNO distribution [20] has been widely used as an important tool for studying multiple particle 
production in inelastic and NSD events. The observed KNO scaling violation for pp collisions at energies 
higher than those achieved at ISR [15, 16] has been explained by the rising number of multiple parton 
scatterings as s —>■ cxd [18,19]. KNO distributions are therefore good tools to exploit how well hadronic 
models can describe the event properties associated to multiple parton scattering. 

By default PYTHIA is set to use multiple parton interactions. Nevertheless, one still has to define 
how the divergency for scatterings with p^ < p^ will be treated by the event generator. PYTHIA 
allows two different phenomenological approaches: simple (MSTP(82)=1) and complex scenarios 
(MSTP(82)=2, 3 or 4). Selecting the complex scenario one has also the choice of selecting different 
matter distributions for the colliding hadrons: uniform (MSTP(82)=2), single Gaussian (MSTP(82)=3) 
and double Gaussian (MSTP(82)=4) matter distributions. 
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Fig. 2: KNO distributions for NSD pp collisions at ^/s = 546 GeV: (a) simple scenario and (b) complex scenario distributions 
compared to data. 
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Figure 2 shows KNO distributions for NSD pp collisions at ^/s = 546 GeV. We compare distribu- 
tions generated by PYTHlA's simple and complex scenarios to UA5 data [16]. Apart from the mentioned 
changes in the setting MSTP(82), all other parameters are set to use PYTHlA's default options, as de- 
scribed in [11]. Figure 2(a) shows that using the simple scenario (MSTP(82)=1), which is the default 
PYTH1A6.214 option [11], the generated distributions fail to reproduce the data, especially in the region 
of high z (z>1.5). This is the region of events with particle multiplicities several times greater than the 
average multiplicity. Distributions generated using the complex scenario vary with the hadronic matter 
distribution selected for each case, as can be seen in figure 2(b) . 

The comparisons of KNO distributions shown in figure 2(b) indicate that the matter distribution 
used to describe the colliding hadrons does affect the probability of particle production in minimum bias 
events. Although in the comparisons shown in figure 2 the best agreement to the data was obtained by 
selecting the complex scenario with the single Gaussian matter distribution option, we shall adopt the 
complex scenario with a double Gaussian matter distribution (MSTP(82)=4) as our preferred choice. 

This is done because by choosing the double Gaussian option, the user is able to control some 
of the properties of this matter distribution. Hadrons described by this distribution have a small core 
region of radius a2 containing a fraction f3 of the total hadronic matter. This core is embedded in a larger 
volume of radius ai containing the remaining fraction of matter, i.e., (1 - /3) of the total hadronic matter. 
The parameter PARP(83) controls the portion f3 of the total hadronic matter assigned to the core of the 
hadron. The ratio a2/ai is given by the parameter PARP(84). By default, PYTHIA sets PARP(83)=0.5 
and PARP(84)=0.2 describing any given hadron as a body with half of its matter concentrated within a 
core which is limited by a radius a.2 = 20% of the hadron radius ai [1 1]. 




As shown in figure 3(a) considerable changes in the high-z tale of the KNO distributions are 
observed as the core radius varies from 20% to 50% and 80% of the radius of the colliding hadrons. 
As the core is made harder and denser (smaller core-radius) the overlap between two colliding cores 
makes high-pt partonic scatterings more likely, yielding higher multiplicity events more often. When 
two relatively softer cores (larger radius) overlap in a collision, the generated activity will be smaller and 
softer, hence producing high-multiplicity events less frequently. 

Figure 3(b) shows a comparison between PHO JET 1.12 and the KNO distribution measured by 
E735 [18, 19] for NSD pp collisions at ^s = 1.8 TeV. Describing hadron collisions using the multiple 
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Pomeron exchange mechanism proposed by the DPM [9, 13] and the QCD picture for high-p^ interac- 
tions, PH0JET1.12 is in good agreement to the data. 



4.42 Pseudorapidity Distribution 

The rate of parton-parton scattering in a hadronic collision is strongly correlated to the observed particle 
multiplicity and the pseudorapidity distribution of produced particles. This happens because multiple 
parton interactions convert part of the collision energy that would otherwise be carried by the fast moving 
system of beam-remnants in the forward regions, into low-p^ particles which populate the central region. 

In PYTHIA, one of the main parameters used to regulate the rate of parton-paiton interactions 
is pt„j„ given by equation 6. Low values of pt,ni,i imply in high rates of parton-parton scatterings and 
hence in high levels of particle multiplicity. Increasing pt^^^ the opposite is expected. 



z 

■= 5 



■ PARP(82)=1.7 
4 PARP(82)=1.9 
o PARP(82)=2.1 



* UA5 - 900 GeV 
NSD interactions 



1.5 



5 - 



2 - 



6 B! 



1.5 



u 1 



S 0.5 



PHOJET1.12 



» . 



• UA5 - 900 GeV 
NSD interactions 



(a) 



(b) 



Fig. 4: Charged particle density distributions.dNch/dry, for NSD pp at a/s = 900 GeV collisions comparing the data [21] to (a) 
PYTHIA6.214 with various pt,„,^ and (b) PHOJET1.12. 



As can be seen in figure 4(a), increasing PARP(82) from 1.7 to 1.9 and 2.1, which effectively 
increases the pt„j„ used by PYTHIA6.214, the charged particle density, dNch/d??, decreases. Notice that 
relatively small changes in PARP(82) (~ 10%) can cause significant variations in the plateau of dNc/i/d?]. 

In PHOJET, multiple Pomeron exchanges predicted by the DPM control the plateau of dNc/j/d?7. 
Similarly to PYTHIA, this model also depends on a p™t-off which is used to connect the soft and hard 
components of a hadronic interaction. PHOJET 1.12 has its default options tuned for p™t-ott — 2.5 GeV. 
Figure 4(b) shows dN ch/di] generated by PHOJET 1.12 with its default cuts, compared to UA5 data [21]. 
There is a good agreement between PHOJET 1.12 predictions and the data. 

4.5 The Underlying Event 

In a hadronic event containing jets, the underlying event (UE) consists of all event activity except the 
two outgoing hard scattered jets [22]. As for minimum bias events, soft interactions and the mechanism 
of multiple parton interaction play an important role in the structure of the underlying event and ought to 
be carefully considered by any model attempting to describe the underlying event. 

The conditions applied to particle selection and to the event region to be investigated are described 
in Ref. [22]. The region transverse to the leading jet is used to study the UE and is defined by 60° < 
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I A(/)| < 120°, where the angular difference in the azimuthal angle is given by Ac/) = ^particle — </'ijet- 

Figure 5(a) shows PYTHIA6.214 - MSTP(82)=4 distributions generated with different values of 
PARP(82), i.e. different Pt„i„, compared to the data for the average charged particle multiplicity in the 
transverse region. Increasing Pt^^^, which corresponds to a decrease on the rate of semi-hard parton 
scatterings, < A^chg > decreases. This effect is similar to the one observed in figure 4(a) for minimum 
bias charged particle density distributions dNch/drj. 
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Fig. 5: Average charged particles multiplicity in the transverse region showing PYTHIA6.214 - MSTP(82)=4 with (a) different 
values of PARP(82) (i.e. pt„„,) and (b) different values of PARP(84) (core-size). 



As shown in figure 5(b), depending on the core size variation (PARP(84)) the plateau level of 
< Nchg > can suffer severe changes. For example, changing PARP(84) from 0.2 to 0.5 reduces the 
plateau of < Nchg > by nearly a factor of two, while a further increase in PARP(84) from 0.5 to 0.8 only 
reduces the plateau by ~ 15%. 

Jets are likely to be produced when there is a core overlap in the hadronic collision. Smaller 
and dense cores imply that events with a core overlap have also a large overlap of less dense matter 
regions which surround the core, and when overlapped generate high rates of soft interactions causing 
the higher plateaus observed in the < N^hg > distributions shown in figure 5(b). Larger cores also imply 
in smaller soft surrounding regions in the colliding hadrons, hence producing lower multiplicity (and p^) 
distributions in the UE. 

Figure 6 shows PHO JET 1.12 predictions compared to data for: (a) average multiplicity in the 
transverse region and (b) average p^^^ in the transverse region. PHO JET 1.12 reproduces reasonably 
well the data for the UE multiplicity distribution, as displayed in figure 6 (a). However, it underestimates 
the average p^ distribution (figure 6(b) ). The measured < p^ > distribution is underestimated by 
PHOJETl. 12 by ~ 20%. 

4.51 UE vs. Minimum Bias 

The CDF measurement shows that the underlying event multiplicity forms a plateau for events with 
^tijct ~ ^ GeV at < A'chg >~ 2.3. Supposing that the transverse region in events with P^i^^.^ > 5 GeV 
is uniform in azimuthal angle and in pseudorapidity rj, this multiplicity corresponds to 3.45 particle 
per unit pseudorapidity. Further corrections to detector effects and low-pt extrapolation [22] implies that 
there are roughly 10 charged particles per pseudorapidity unit with pt > GeV in the underlying event. 
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Fig. 6: PHOJET1.12 predictions compared to CDF data for: (a) average multiplicity in the transverse region and (b) average 
Pj in the transverse region. 



In pp collision at 1.8 TeV, the minimum bias density, which has also been measured by CDF, gives 
dNcfc/dr/ ~ 4 for \7]\ < 1 [17], while the equivalent density for the underlying event is at least a factor 
of two larger. This comparison, though not highly accurate due to the uncertainties in estimating the 
particle density for the underlying event (i.e. extrapolation to low-p^ and several assumptions made on 
the particle distribution in cf) and rj), clearly shows that the underlying event in hard scattering processes 
(Piijet ^ ^ GeV) has much more activity than an average minimum bias event. 

4.6 PYTHIA6.214 - Tuned VS. PHOJET1.12 

Combining the effects of variations in p(_^^.^ and in the core-size we obtained a set of PYTHIA6.214 
parameters which considerably improves PYTHIA's description of minimum bias and underlying event 
distributions. Our tuned parameters for PYTHIA6.214 are displayed in table 1. 



PYTHIA6.214 - tuned 


ISUB: 11,12,13,28,53,68 


QCD 2 — > 2 partonic scattering 


94,95,96 


+ non-diffractive + double diffractive 


MSTP(51)=7 


CTEQ5L - selected p.d.f. 


MSTP(81)=1 


multiple interactions 


MSTP(82)=4 


complex scenaiio 




+ double Gaussian matter distribution 


PARP(82)=1.8 


pt„,i,^ parameter 


PARP(84)=0.5 


core radius: 50% of the 




hadronic radius 



Table 1: PYTHIA6.214 tuned parameters for minimum bias and the underlying event. 



Figure 7 shows predictions generated by PYTHIA6.214 - tuned and default, and PH0JET1.12 
compared to some minimum bias and underlying event distributions. The description of both minimum 
bias and underlying event distributions is improved by using PYTHIA6.214 - tuned compared to the 
predictions generated by PYTHIA's default settings. Notice that PYTHIA6.214 - tuned and PH0JET1.12 
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Fig. 7: (a) KNO distributions for NSD pp collisions at ^/s = 900 GeV; (b) dNch/d?; for NSD pp at \/s = 1-8 TeV; (c) dNch/d?? 
at = for a wide range of -^s; and (d) < Nchg > in the transverse region. 
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can generate very different predictions when extrapolated to higher energies, as shown in fig. 7(c). 
4.7 LHC Predictions 
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Fig. 8: (a) Charged particle density distributions, dNch/dTj, for NSD pp collisions at ^/i = 14 TeV; (b) Average multiplicity in 
the underlying event for jet events in pp collisions at -^s = 14 TeV. 



Figure 8(a) shows dNc/i/dr/ distributions for minimum bias pp collisions at 14 TeV generated 
by PH0JET1.12 and PYTHIA6.214 - tuned. The charged particle density generated by PH0JET1.12 
and PYTHIA6.214 - tuned at = is 5.13 and 6.82, respectively. In the central region (|?7| < 2.5) 
dNch/dr] is ~ 5.5 and ~ 7, respectively for PH0JET1.12 and PYTHIA6.214 - tuned. Contrasting to the 
agreement shown for pp collisions at = 1.8 TeV in figure 7(b), at the LHC PYTHIA6.214 - tuned 
generates ~ 27% more charged particle density in the central region than PH0JET1.12. 

Compared to the charged particle density dNc^/dr/ measured by CDF at 1.8 TeV (figure 7(b) ), 
PYTHIA6.214 - tuned indicates a plateau rise of 70% at the LHC in the central region while PHO- 
JET1.12 suggests a smaller rise of ~ 35%. 

Figure 8(b) displays PYTHIA6.214 - tuned and PH0JET1.12 predictions for the average particle 
multiplicity in the UE for pp collisions at the LHC (charged particles with pt > 0.5 GeV and |r/| < 1). 
The distributions generated by the two models are fundamentally different. Excepting the events with 
Ptije, ^ 3 GeV, PYTHIA6.214 - tuned generates greater activity than PH0JET1.12 in the UE. The average 
multiplicity in the UE for P^^^ > 10 GeV reaches a plateau at 6.5 charged particles according to 
PYTHIA6.214 - tuned and ~ 3.0 according to PH0JET1.12. Compared to the UE distributions measured 
by CDF at 1.8 TeV (figure 7(d) ), PYTHIA6.214 - tuned indicates a plateau rise of ~ 200% at the LHC 
while PHO JET 1.12 suggests a much smaller rise of ~ 40%. 

At the LHC, the minimum bias predictions generated by PYTHIA6.214 - tuned and PH0JET1.12 
for the central plateau of dNc/i/dry, indicate a rise of ~ 70% and ~ 35%, respectively. These are smaller 
than the predicted increase for the UE suggested by both models. As discussed previously, at the Teva- 
tron, for events with P^^^ > 10 GeV the particle density in the underlying event is at least a factor of 
two larger than the equivalent minimum bias prediction. Using similar assumptions as those adopted in 
the analysis for the CDF data, LHC events with P^^^ > 10 GeV are predicted to have a charged particle 
density dNc/j/dr/ of ~ 29 charged particles per pseudorapidity unit according to PYTHIA6.214 - tuned 
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and ~ 13 according to PH0JET1.12. In other words, for F^^^ > 10 GeV the UE at the LHC is predicted 
to have a particle density ~ 4 times larger than its equivalent minimum bias prediction according to 
PYTHIA6.214 - tuned, and ~ 2 times larger according to PH0JET1.12. 

Therefore PYTHIA6.214 - tuned predicts not only that the UE particle density will increase at the 
LHC, but it will also increase its activity compared to the equivalent minimum bias distribution. On the 
other hand, PH0JET1.12 estimates that the increase in charged particle density in the UE at the LHC 
will follow the same rate to the minimum bias density measured at the Tevatron. In both cases however, 
the underlying event density is greater than its equivalent minimum bias counterpart. 

5. USING CORRELATIONS IN THE TRANSVERSE REGION TO STUDY THE UNDERLY- 
ING EVENT IN RUN 2 AT THE TEVATRON ^ 

5.1 Introduction 

Fig. 9 illustrates the way QCD Monte-Carlo models simulate a proton-antiproton collision in which a 
"hard" 2-to-2 parton scattering with transverse momentum, pr(hard), has occurred. The resulting event 
contains particles that originate from the two outgoing partons (plus initial and final-state radiation) and 
particles that come from the breakup of the proton and antiproton (i.e., "beam-beam remnants"). The 
"underlying event" is everything except the two outgoing hard scattered "jets" and receives contributions 
from the "beam-beam remnants" plus initial and final-state radiation. The "hard scattering" component 
consists of the outgoing two jets plus initial and final-state radiation. 




Fig. 9: Illustration of the way QCD Monte-Carlo models simulate a proton-antiproton collision in which a "hard" 2-to-2 parton 
scattering with transverse momentum, pT(hard), has occurred. The resulting event contains particles that originate from the 
two outgoing partons (plus initial and final-state radiation) and particles that come from the breakup of the proton and antiproton 
(i.e., "beam-beam remnants"). The "underlying event" is everything except the two outgoing hard scattered "jets" and consists 
of the "beam-beam remnants" plus initial and final-state radiation. The "hard scattering" component consists of the outgoing 
two jets plus initial and final-state radiation. 

The "beam-beam remnants" are what is left over after a parton is knocked out of each of the initial 
two beam hadrons. It is the reason hadron-hadron collisions are more "messy" than electron-positron 
annihilations and no one really knows how it should be modeled. For the QCD Monte-Carlo models 
the "beam-beam remnants" are an important component of the "underlying event". Also, it is possible 
that multiple parton scattering contributes to the "underlying event". Fig. 10 shows the way PYTHIA 
[1, 23,24] models the "underlying event" in proton-antiproton collisions by including multiple parton 
interactions. In addition to the hard 2-to-2 parton-parton scattering and the "beam-beam remnants", 
sometimes there is a second "semi-hard" 2-to-2 parton-parton scattering that contributes particles to the 
"underlying event". 

Of course, from a certain point of view there is no such thing as an "underlying event" in a 
proton-antiproton collision. There is only an "event" and one cannot say where a given particle in 
the event originated. On the other hand, hard scattering collider "jet" events have a distinct topology. 

■^Contributed by: A. Cruz and R. Field 
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Fig. 10: Illustration of the way PYTHIA models the "underlying event" in proton-antiproton collisions by including multiple 
parton interactions. In addition to the hard 2-to-2 parton-parton scattering with transverse momentum, pT(hard), there is a 
second "semi-hard" 2-to-2 parton-parton scattering that contributes particles to the "underlying event". 

On the average, the outgoing hadrons "remember" the underlying the 2-to-2 hard scattering subprocess. 
An average hard scattering event consists of a collection (or burst) of hadrons traveling roughly in the 
direction of the initial beam particles and two collections of hadrons (i.e., "jets") with large transverse 
momentum. The two large transverse momentum "jets" are roughly "back-to-back" in azimuthal angle. 
One can use the topological structure of hadron-hadron collisions to study the "underlying event" [22, 
25-29]. Here we study the "underlying event" in the Run 2 using the direction of the leading calorimeter 
jet (JetClu, R = 0.7) to isolate regions of rj-cl) space that are sensitive to the "underlying event". 




Fig. 11: Illustration of correlations in azimuthal angle relative to the direction of the leading jet (JetClu, R — 0.7) in the 
event, jet#l. The angle = (p — f^jcti is the relative azimuthal angle between charged particles and the direction of jet# 
1. The "toward" region is defined by |A0| < 60° and \r]\ < 1, while the "away" region is |A<^| > 120° and I?;! < 1. The 
"transverse" region is defined by 60° < |A0| < 120° and I?;] < 1. Each of the three regions "toward", "transverse", and 
"away" have an area in 7;-(j!>space of ArjA<j) = 4n/3. We examine charged particles in the range pr > 0.5 GeV/c and |r/| < 1, 
but allow the leading jet to be in the region j7;(jet#l)| < 2. 

As illustrated in Fig. 11, the direction of the leading jet, jet# 1, is used to define correlations in 
the azimuthal angle, 0. The angle A(j) = 4> — 0jcti is the relative azimuthal angle between a charged 
particle and the direction of jet# 1. The "toward" region is defined by |A(/)| < 60° and \ri\ < 1, while the 
"away" region is |A0| > 120° and \ri\ < 1. The "transverse" region is defined by 60° < |A0| < 120° 
and \7]\ < 1. The three regions "toward", "transverse", and "away" are shown in Fig. 11. Each region 
has an area in r/-0space of ArjAcj) = Air/ 3. The "transverse" region is perpendicular to the plane of the 
hard 2-to-2 scattering and is therefore very sensitive to the "underlying event". We restrict ourselves to 
charged particles in the range pT > 0.5 GeV/c and \r]\ < 1, but allow the leading jet that is used to define 
the "transverse" region to have |r/(jet#l)| < 2. 

In this analysis we look in more detail at the two "transverse" regions defined in Fig. 12. The 
overall "transverse" region in Fig. 11 corresponds to combining the "transverse 1" and "transverse 2" 
regions. Comparing these two "transverse" regions on an event-by-event basis provides a closer look 
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Jet #2 Direction 



Fig. 12: Illustration of correlations in azimuthal angle relative to the direction of the leading jet (highest ET jet) in the event, 
jet#l. The angle A<^ = — 0jcti is the relative azimuthal angle between charged particles and the direction of jet#l. The two 
"transverse" regions 60° < < 120° and 60° < —Acp < 120° are referred to as "transverse 1" and "transverse 2". Each of 
the two "transverse" regions have an area in 77-0 space of At^A^ = 47r/6. The overall "transverse" region defined in Fig. 11 
corresponds to combining the "transverse 1" and "transverse 2" regions. Events in which there are no restrictions placed on the 
on the second highest Et jet, jet#2, are referred to as "leading jet" events (left). Events with at least two jets with |f7(jet)| < 2, 
where the leading two jets are nearly "back-to-back" (\A(j)\ > 150°) with ET{jet^2) / Erijet^l) > 0.8 are referred to as 
"back-to-back" events, (right). 



at the "underlying event". Here we refer to events in which there are no restrictions placed on the 
second highest Et jet, jet#2, as "leading jet" events. Our previous analysis of the "underlying event"we 
only considered "leading jet" events [22, 27-29]. In this analysis we define a second class of events. 
Events with at least two jets with |ry(jet)| < 2, where the leading two jets are nearly "back-to-back" 
{\A4>\ > 150°) with ETijet#2)/ETijet#l) > 0.8 are referred to as "back-to-back" events. "Back-to- 
back" events are a subset of the "leading jet" events. The idea here is to suppress hard initial and final- 
state radiation thus increasing the sensitivity of the "transverse" region to the "beam-beam remnants" 
and the multiple parton scattering component of the "underlying event". 

As in our published Run 1 analysis [22] we consider charged particles only in the region pT > 
0.5 GeV/c and |r/| < 1 where the COT efficiency is high and compare uncorrected data with PYTHIA 
Tune A [27,28] and HERWIG [30-32] after detector corrections {i.e., CDFSIM). Systematic errors are 
calculated in the same way as in our Run 1 analysis. We generate every plot twice, once with our chosen 
track selection cuts and again with the very tight track cuts. The change in each point in every plot due to 
this tighter cut is used as a measure of the systematic error and is added in quadrature with the statistical 
error to form the overall error. 

5.2 Transverse Average Pt vs Nchg 
5.21 Definition 

We study the average transverse momentum of charged particles in the "transverse" region as a function 
of the number of charged particles in the "transverse" region for pT > 0.5GeV/c and \r]\ < 1. The 
average transverse momentum, (pT) , is formed, on an event-by-event basis, and then plotted as a function 
of the charged multiplicity. The idea here is to look for correlations between multiplicity and (px). If, 
for example, there is a mixture of "hard" and "soft" events then one expects that (pr) will increase with 
multiplicity because demanding a large multiplicity will preferentially select the "hard" process that also 
has a lai-ger (pT). On the other hand, it may be possible to get a high multiplicity in a "soft" collision so 
the rate that (pT) rises with multiplicity is a rough measure of the "hard" and "soft" mixture. The steeper 
the slope the larger the "hard" component. There is a very nice published CDF Run 1 analysis that looks 
at this in "min-bias" collisions [33], but it has never previously been studied in the "transverse" region 
of a "hard" scattering process. 
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Number of Charged Particles 



Fig. 13: Run 2 data on the average transverse momentum as a function number of particles for charged particles with pT > 
0.5 GeV/c and < 1 in the "transverse" region for "leading jet" events defined in Fig. 12 with 30 < i5T(jet#l) < 70GeV 
and 130 < BrQet^l) < 250 GeV. Also shown are the data on the average transverse momentum as a function of the number 
particles for charged particles with pr > 0.5 GeV/c and < 1 for "min-bias" collisions at 1.96 TeV. The theory curves 
correspond to PYTHIA Tune A at 1.96 TeV (after CDFSIM). 




Fig. 14; Run 2 data on the average transverse momentum as a function of the number of particles for particles for charged 
particles with pr > 0.5 GeV/c and |r;| < 1 in the "transverse" region for "leading jet" events and for "back-to-back" events 
defined in Fig. 12 with 30 < ii'T(jet#l) < 70 GeV. Also shown are the data on the average transverse momentum as a 
function of the number particles for charged particles with pr > 0.5 GeV/c and < 1 for "min-bias" collisions at 1.96 TeV. 
The theory curves correspond to PYTHIA Tune A at 1.96 TeV (after CDFSIM). 
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5.22 Overall Transverse Region 

Fig. 5.21 shows uncorrected Run 2 data on the ^t) of charged particles versus the number of charged 
particles in min-bias collisions and in the "transverse" region for "leading jet" events with 30 < 
£;^(jet#l) < 70GeV and 130 < ^T(jet#l) < 250 GeV compared with PYTHIA Tune A (after 
CDFSIM). The data suggest that there is more "hard" scattering in the "transverse" region (i.e., initial 
and final-state radiation) than there is in an average "min-bias" collision. 

Fig. 14 shows the (px) of charged particles versus the number of charged particles in "min- 
bias" collisions and in the "transverse" region for "leading jet" and "back-to-back" events with 30 < 
^y(jet#l) < 70 GeV compared with PYTHIA Tune A (after CDFSIM). The "transverse" region for 
the "back-to-back" events looks more like an average "min-bias" collision, which is exactly what one 
expects since the "back-to-back" requirement suppress hard initial and final-state radiation. 
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Fig. 15: Run 2 data on the average transverse momentum as a function of the number of particles for charged particles with 
PT > 0.5 GeV/c and \r]\ < 1 in the "transverse" region for "leading jet" events defined in Fig. 12 with 30 < -ErQet^l) < 
70 GeV and 130 < £'T(jet#l) < 250 GeV compared to HERWIG at 1.96 TeV (after CDFSIM). 
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Fig. 16: Run 2 data on the average transverse momentum as a function of the number of particles for charged particles with 
Pt > 0.5 GeV/c and \ri\ < 1 in the "transverse" region for "leading jet" and "back-to-back" events defined in Fig. 12 with 
30 < £T(jet#l) < 70 GeV compared to HERWIG at 1.96 TeV (after CDFSIM). 



Fig. 15 compares HERWIG (after CDFSIM) with the data on the (px) of charged particles in the 
"transverse" region versus the number of charged particles in the "transverse" region for "leading jet" 
events with 30 < £:T(jet#l) < 70 GeV and 130 < ^T(jet#l) < 250 GeV. 
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Fig. 16 compares HERWIG (after CDFSIM) with the data on the l^pT) of charged particles in 
the "transverse" region versus the number of charged particles in the "transverse" region for "leading 
jet" and "back-to-back" events with 30 < £'T(jet#l) < 70GeV. HERWIG (without multiple parton 
interactions) does not describe the data as well as PYTHIA Tune A (with multiple parton interactions). 
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Fig. 17: Run 2 data on the average number of particles in the "transverse 2" region defined in Fig. 12 as a function of the number 
of particles in the "transverse 1" region for charged particles with pr > 0.5 GeV/c and < 1 for "leading jet" events with 
30 < £T(jet#l) < 70 GeV (top) and 130 < £;T(jet#l) < 250 GeV (bottom). The theory curves correspond to PYTHIA 
Tune A and HERWIG at 1.96 TeV after CDFSIM. 



5.23 Transverse 1 versus transverse 2 

Fig. 17 shows the number of charged particles in the "transverse 2" region versus the number of charged 
particles in the "transverse 1" region for "leading jet" events with 30 < E'rQetT^l) < 70 GeV and 
130 < £:T(jet#l) < 250 GeV compared with PYTHIA Tune A and HERWIG after CDFSIM. 

Fig. 18 shows the (pT) of charged particles in the "transverse 2" region versus the number of 
charged particles in the "transverse 1" region for "leading jet" events with 30 < E'rOet^l) < 70 GeV 
and 130 < £:T(jet#l) < 250 GeV compared with PYTHIA Tune A and HERWIG after CDFSIM. 

Fig. 19 shows the number and (pT) of charged particles in the "transverse 2" region versus the 
number of charged particles in the "transverse 1" region for "leading jet" and "back-to-back" events with 
30 < £:T(jet#l) < 70 GeV compared with PYTHIA Tune A (after CDFSIM). 

Fig. 20 shows the number and (pT) of charged particles in the "transverse 2" region versus the 
number of charged particles in the "transverse 1" region for "leading jet" and "back-to-back" events 
with 30 < ^T(jet#l) < 70 GeV compared with HERWIG (after CDFSIM). HERWIG (without mul- 
tiple parton interactions) does not do nearly as well describing the "transverse 2" versus "transverse 1" 
correlations seen in the data as does PYTHIA Tune A (with multiple parton interactions). 

5.3 Summary 

This analysis takes a closer look at the "underlying event" in hard scattering proton-antiproton collisions 
at 1.96 TeV. We look only at the charged particle component of the "underlying event" and restrict the 
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Fig. 18: Run 2 data on the average transverse momentum of particles in the "transverse 2" region defined in Fig. 12 as a function 
of the number of particles in the "transverse 1" region for charged particles with pt > 0.5 GeV/c and < 1 for "leading jet" 
events with 30 < £T(jet#l) < 70GeV (top) and 130 < £;T(jet#l) < 250 GeV (bottom). The theory curves correspond to 
PYTHIA Tune A and HERWIG at 1.96 TeV after CDFSIM. 





Fig. 19: Run 2 data on the average number of particles (top) and the average transverse momentum of particles (bottom) in the 
"transverse 2" region defined in Fig. 12 as a function of the number of particles in the "transverse 1" region for charged particles 
with pt > 0.5 GeV/c and |j7| < 1 for "leading jet" events and "back-to-back" events with 30 < i5T(jet#l) < 70 GeV. The 
theory curves correspond to PYTHIA Tune A at 1.96 TeV (after CDFSIM). 
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charged particles to be in the range px > 0.5 GeV/c and \rj\ < 1. We use the direction of the leading 
calorimeter jet in each event to define two "transverse" regions of r]-(j) space that are very sensitive to 
the "underlying event". Comparing these two "transverse" regions on an event-by-event basis provides 
more details about the "underlying event". In addition, by selecting events with at least two jets that 
are nearly "back-to-back" (|A0| > 150°) we are able to look closer at the "beam-beam remnants" and 
multiple parton interaction components of the "underlying event". PYTHIA Tune A (with multiple 
parton interactions) does a good job in describing the "underlying event" (i.e., "transverse" regions) 
for both "leading jet" and "back-to-back" events. HERWIG (without multiple parton interactions) does 
not have enough activity in the "underlying event", which was also observed in our published Run 1 
analysis [22]. 
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Fig. 20: Run 2 data on the average number of particles (top) and the average transverse momentum of particles (bottom) in the 
"transverse 2" region defined in Fig. 12 as a function of the number of particles in the "transverse 1" region for charged particles 
with pt > 0.5 GeV/c and |r;| < 1 for "leading jet" events and "back-to-back" events with 30 < ET{}et#l) < 70 GeV. The 
theory curves correspond to HERWIG at 1.96 TeV (after CDFSIM). 



The data presented here show interesting correlations between the two "transverse" regions defined 
in Fig. 12. The charged multiplicity and the ^t) in the "transverse 2" region increases with the charged 
multiplicity in the "transverse 1" region. This is a new type of correlation. It might simply be due to a 
high multiplicity in "transverse 1" biasing in favor of a harder 2-to-2 scattering (i.e., higher pj'(hard)) 
which would result in a higher multiplicity and larger (pr) in "transverse 2". However, we have seen 
in previous studies [22, 25-29] that the average charged particle density in the "transverse" region does 
not change much as one increases ET{iet#l)- It is possible that the "transverse 1" versus "transverse 
2" correlations arise from multiple parton interactions. A large multiplicity in the "transverse 1" region 
would indicate a small impact parameter collision has occurred with several multiple parton scatterings 
which would then cause an increased multiplicity and (px) in the "transverse 2" region. The fact that 
PYTHIA Tune A (with multiple parton interactions) agrees with the data better than HERWIG (without 
multiple parton interactions) is very interesting. However, much more work is necessary to actually 
pinpoint the source of the "transverse 1" versus "transverse 2" correlations. 
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6. SIMULATION OF THE QCD BACKGROUND FOR tt ANALYSES AT THE TEVATRON 
WITH A Z±+ JETS FINAL STATE ^ 

6.1 Introduction 

The top quark mass and the cross section of top quark pair production are important measurements at 
the Tevatron. Compared to Run I, the precision of these measurements is expected to be much higher 
in Run II Ref. [34]. This requires an accurate understanding of all the important backgrounds. In the 
semileptonic top decay channel with one lepton and four jets plus missing energy in the final state, 
the QCD background is one of the main backgrounds. The cross section of this background is large 
compared to the cross section of tt production, but can be reduced strongly with appropriate selection 
criteria. In the past, this background was modeled using + jets events with non-isolated leptons 
coming from experimental data. It is however important that in ti analyses with + jets final states 
the lepton is isolated in order to reduce backrounds with non-isolated leptons. It is not possible to get 
a clean sample of QCD events with isolated leptons directly from the data since these would be mixed 
with ti and W + jets events. Therefore, in the following, we instead model the QCD background via 
Monte Carlo with hhqq and hhgg production and subsequent leptonic i?-decay. 

6.2 Simulation Procedure 

o Hard Process: the hard 2 ^ A process is generated with ALPGEN Ref. [35]. The final state 
consists of a bh pair together with two other partons (light quarks or gluons). In the following, these 
events are referred to as "bb2j" events. They are generated with ^fs^ = 1960 GeV, CTEQ5L, 
and rrib = 4.75 GeV. In order to produce the events in kinematically interesting regions, we 
apply the following phase space cuts for all four partons, including the b and b: prij) > 10 GeV, 
\vU)\ < 3 and AR{j,j) > 0.4. In addition, we require at least one 6-quark with a transverse 
momentum bigger than 30 GeV in order to provide phase space for the leptonic i?-decay. The 
bb2j cross section including these cuts is about 30 nb. 

O Fragmentation: during the fragmentation with PYTHIA Ref. [10] a high pT electron or muon is 
produced in a leptonic S-decay. In order to increase the number of events with energetic leptons, 
the fragmentation is repeated exactly forty times. For this number of repetitions, we obtain about 
one event with a high pT lepton. A constant number of repetitions should not bias the pT spectrum 
of the lepton. Events preselected with a lepton with px of at least 15 GeV are written to a HEPEVT 
file. The high pT lepton requirement reduces the cross section to 14.5 pb (288864 events, the first 
number in Table 2 below), presumably with a large theoretical uncertainty^. 

o CDF detector simulation: The standard CDF software package reads the HEPEVT file, and per- 
forms a detailed simulation of the CDF detector response. After the simulation, the events are 
reconstructed with the same algorithms as used to reconstruct the Run n data. Finally, the top 
working group's standard ntuple is written which we examine in the following. 

6.3 Comparison with Data 

O Preselection: we select events with exactly one lepton' with Et > 20 GeV and with missing 
transverse energy of at least 20 GeV. Furthermore, we ask for at least four jets with |r/(j)| < 2.4 
and Exij) > 8 GeV. Three jets per event are required to be high quality jets with \'r]{j)\ < 2.0 
and E-r(j) > 15 GrV. In addition, we remove events with more than two 6-quarks in case of 

'Contributed by: V. DroUinger 

*The cross section of bb2j is calculated at leading order. Additional background contributions from similar processes like 
cc2j are not included. 

^We use the CDF baseline selection criteria for high quality lepton identification without the isolation cut. Events with 
cosmic muons or more than one lepton are rejected as well. 
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simulated events. This mitigates a potential bias from the repetition of the fragmentation. In case 
of the data, we use good runs only. 

o QCD region: in the plane of the lepton isolation variable iso ^ and the distance of the lepton 
to the closest jet AR„iin{l^ , j) we select a region, where we expect almost purely QCD events 
in the data. In this region the lepton is non-isolated and is well inside a jet: iso > 0.15 and 
ARminil"^ , j) < 0.15 rad. Isolated leptons coming from real W"^ — > l^^u decays, are found 
typically at low iso and at any value of ARmin{l^ , j) which means, those leptons are usually 
isolated and also separated from jets. 

O Comparison: in order to evaluate the simulated bb2j event simulation, we compare the events in 
the QCD region with events from the data which have undergone the same preselection. The data 
sets are called btopOg, btopig, btopOj, and btoplj with the number of events listed in Table 2. All 
events correspond to an integrated luminosity of 126 pb~^. In case of the simulated bb2j events 
"all events" is the number of events after the full simulation and reconstruction. We compare the 
most relevant kinematic variables in Figs. 22 and 21. 



event type 


bb2j 


data 


all events 


288864 


955870 


preselection 


846 


1669 


QCD region 


666 


1043 



Table 2: Number of events after each major selection step. 

In order to validate the simulation of bb2j events, we have studied the shapes of kinematic distributions. 
Whereas the jets (Fig. 22) are more related to the hard QCD process, the lepton and the neutrino, seen 
as missing Et, (Fig. 21) depend rather on the B -decay and fragmentation, respectively. However, the 
kinematic properties of the jets, the lepton, and the neutrino are correlated. There are no major differences 
between simulated bb2j events and the CDF data. 

6.4 Summary 

Top physics at the Tevatron has entered the Run II phase with high luminosity and upgraded detectors, 
and therefore the measurements obtained in the top sector will reach much higher precision. This requires 
a good understanding of all relevant background processes. The use of new Monte Carlo tools enable us 
to simulate the QCD background. 

We have described the bb2j event generation procedure. After the hard process with four partons 
in the final state is generated with ALPGEN, the high pT lepton is obtained from a leptonic 5-decay 
during the fragmentation with PYTHIA. 

We define a region, where we expect to have a clean QCD sample and compare the fully simulated 
events with the CDF Run n data. All of the kinematic variables studied, compare well. From this 
comparison we conclude, that bb2j events can be used to model the QCD background with a + jets 
final state. 

*The repetition of the fragmentation would cause an artificial enhancement of additional bb pairs coming from gluon splitting 
with subsequent leptonic -B-decay. 

'The lepton isolation variable iso is defined as the fractional calorimeter isolation Et in a cone of -R = 0.4 around the 
lepton. 
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Fig. 21: Comparison of bb2j events with 126 pb~^ of CDF Run II data: normalized distributions of Et of the lepton in respect 
to the closest jet axis, Et and 77 of the lepton, missing Et, Acp between missing Et and closest jet, and multiplicities for all 
selected jets. 

7. MONTE-CARLO DATABASE " 
7.1 Problem description 

One of the most general problems for the experimental high energy physics community is Monte-Carlo 
(MC) simulation of physics processes. There are numerous publicly available MC generators. However, 
the correct MC simulation of complicated processes requires in general rather sophisticated expertise on 
the part of their users. Often, a physics group in an experimental collaboration requests experts and/or 
authors of MC generators to create MC samples for a particular process. Furthermore, it is common that 
the same physics process is investigated by various physics groups in need of the the same MC event 
samples. The main motiviation behind the Monte-Carlo Database (MCDB) project is to make MC event 
samples, as prepared by experts, available for various physics groups. In this contribution we present a 
version of the MCDB that is already operative in the CMS collaboration, and discuss future plans. 

There are a number of useful aspects to a central MC Database that motivate its establishment. 

1. Correct and reliable MC event generation of most processes of interest requires considerable ex- 
pertise. Moreover, most MC generators, in particular those calculating higher order perturbative 
corrections, require significant amounts of computer resources. By means of the MCDB, samples 
prepared by experts can be distributed easily and used as many times as needed. 

'"Contributed by: L. Dudko and A. Sherstnev 
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Fig. 22: Comparison of bb2j events with 126 pb ^ of CDF Run II data: normalized Et and i] distributions of the four leading 
jets, and the AR distribution of the pair of closest jets to each other. All jet energies shown, are corrected. 
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2. Public availability of common event files helps speed up the validation procedure of the events. 

3. A central and public location where well-documented MC events and MC generators can be found 
would be very useful. It would also allow rapid communication between authors of MC events and 
their users. 

4. The same MC samples for SM processes can be used for multiple purposes, e.g. to study back- 
grounds to various new physics processes. 

5. Files containing detector and beam-related backgrounds can also be kept in a common location. 

7.2 History 

The first MC Database (version 0) - named PEVLIB [36] - was established at CERN, on the AFS file- 
system. This database provided CompHEP [37] parton level events for CMS users, but lacked a special 
interface. Rather, it was built as a set of directories where event samples were stored. The sample 
documentation consisted of ASCII files (README) located in the same directories as the event files. 

The next version (version 0.5) of the MC Database was estabUshed at Fermilab. This database was 
spUt in two independent parts: 

- MC event files, stored via the FNAL tape system ENSTORE [38] 

- the events' documentation, publicly available via the web [39] 

The latest version of the MCDB (version 1.0), described in this contribution, CMS MCDB [40], was 
created at first for use by the CMS collaboration. This database includes web interfaces both for event 
files (enabling download and upload) and documentation. Its main goal is to store events, only at the 
parton level, generated by MC experts, for use by the LHC community. Note that all files from PEVLIB 
have now been moved to the CMS MCDB. 

7.3 General concepts and practical realization of CIMS IVICDB. 

The MCDB must provide persistent storage of event samples, with pubUc and convenient interfaces for 
users - mainly consisting of LHC experimentalists - and authors of MC generators, or other experts. The 
main features present in the CMS MCDB are: 
o the MCDB is based on web technologies. 

o the MCDB stores parton level events with a standard interface to the next level of simulation (based 

on Les Houches Accord I), 
o the MCDB stores detailed documentation for each set of event samples. 
o the MCDB make rapid communication between users and experts possible via its web pages. 
O the MCDB is divided in two zones: 

> a public area, for all users interested in using MC events. Users can find all necessary infor- 
mation about available event samples and their generators, and download the corresponding 
event files. 

> a restricted area for authors. In this area, authors of MC generators or MC experts can 
change the content of MCDB dynamically - e.g. upload events for new processes, create 
and edit documents, reply to user's comments and questions and upload new event files and 
generators. 

o the MCDB requires users to reference the authors of the event sample in case the events are used 
in a physics analysis. 

The CMS MCDB at CERN is realized as a dedicated web site [40]. The structure of the site 
corresponds to the stored physics processes. For example, all event files relevant for studies of top quark 
production are collected in the category "TOP". This category contains files with events involving top 
quark production via QCD and the electroweak interaction. When a user clicks on a reference in a certain 
category, he/she will see all so-called "articles" related to the selected physics process. 
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The main units of CMS MCDB are HTML documents called articles that describe the event sam- 
ples. The articles are created by authors; these are authorized users who can upload new event samples 
to the CMS MCDB. There are only a few restrictions on the articles in the CMS MCDB. An author can 
create a new article on the basis of a template - this is a web form, where the author fills out various 
fields: author name, abstract, category name, article body, etc. The body and the abstract of the article 
may include HTML tags for more flexible visualization. 

It is easy for users to understand which process is present in the CMS MCDB, via a click on a 
reference to the article that describes the corresponding event samples. Users can download these event 
files directly from the article web page, together with all other files which authors have uploaded to the 
article. 

The web intereface of the present version of the MCDB has the following features, in brief: 

A. Authorized authors can 

o upload files (with events or a MC generator code). There are two different methods to upload file 

(through a web browser and/or directly from the AFS file system). 
o document each set of files in a new article. 

O reply to comments from users (MCDB allows one to organize a special forum for each article). 

O change any of the parameters and content of MCDB, according to permissions (including the web 
design). 

B. Users can 

o read documentation (MCDB articles) for files uploaded by authors, 
o download files. 

o search the CMS MCDB web site (enabled for articles only, not for the event samples themselves), 
o ask questions about a particular document, 
o send new articles to moderators. 

Event files in the CMS MCDB are stored on the AFS filesystem [41]. If a user has access to AFS, 
he/she may download files from AFS directly. To become an authorized author one should send a request 
(by e-mail) to the administrators " of the CMS MCDB. 

7.4 Future plans 

The CMS MCDB is designed to store parton level events. This implies that the size of event files should 
not be too large (typically smaller then 100Mb). The expected number of physics processes in the CMS 
MCDB is several hundred. Note that CMS MCDB is not a SQL database, so that its search engine 
cannot serve complex queries, only keyword phonetic searches. These aspects of the CMS MCDB are 
not a limitation at present. However, we expect that in a few years users will request a more powerful 
MCDB where these restrictions will be removed. The next version of the MCDB now under development 
for use, within the CERN LCG framework, by all CERN collaborations. It will be described in the near 
future [42]. 
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8. RESUMMATION AND SHOWER STUDIES 



8.1 Introduction 

The transverse momentum of a colour-singlet massive particle produced in a hadronic collision provides 
important information on perturbative and nonperturbative effects. A process like qq corresponds 
to p±z = 0, while higher-order processes provide p± kicks as the Z'^ recoils against quarks and gluons. 
At large p±z values the bulk of the p± comes from one hard emission, and perturbation theory is a 
reasonable approach. In the small-pj_z region, on the other hand, many emissions can contribute with p± 
kicks of comparable size, and so the order-by-order approach is rather poorly convergent. Furthermore, 
in this region nonperturbative effects may start to become non-negligible relative to the perturbative ones. 

The traditional solution has been to apply either an analytical resummation approach or a numeri- 
cal parton-shower one. These methods to some extent are complementary. The norm today is for showers 
to be based on an improved leading-log picture, while resummation is carried out to next-to-leading logs. 
However, resummation gives no information on the partonic system recoiling against the Z°, while show- 
ers do, and therefore can be integrated into full-fledged event generators, allowing accurate experimental 
studies. In both approaches the high-p^ tail is constrained by fixed-order perturbation theory, so the 
interesting and nontrivial region is the low-to-medium-p^ one. Both also require nonperturbative input 
to handle the low-p_L region, e.g. in the form of a primordial k± carried by the initiator of a shower. 

One of the disconcerting aspects of the game is that a large primordial k± seems to be required 
and that the required value of this primordial k± can be dependent on the kinematics of the process be- 
ing considered. Confinement of partons inside the proton implies a {k±) « 0.3 GeV, while fits to Z° 
data at the Tevatron favour ^ 2 GeV [43] (actually as a root-mean-square value, assuming a Gaussian 
distribution). Also resummation approaches tend to require a non-negligible nonperturbative contribu- 
tion, but that contribution can be determined from fixed-target data and then automatically evolved to the 
kinematical region of interest. In this note we present updated comparisons and study possible shower 
modifications that might alleviate the problem. We will use the two cases of qq Z° and gg H° (in 
the infinite-top-mass Umit) to illustrate differences in quark and gluon evolution, and the Tevatron and 
the LHC to quantify an energy dependence. 

8.2 Comparison Status 

A detailed comparison of analytic resummation and parton showers was presented in [43]. For many 
physical quantities, the predictions from parton shower Monte Carlo programs should be nearly as precise 
as those from analytical theoretical calculations. In particular, both analytic and parton shower Monte 
Carlos should accurately describe the effects of the emission of multiple soft gluons from the incoming 
partons. 

Parton showers resum primarily the leading logs, which are universal, i.e. process-independent, 
depending only on the initial state. An analytic resummmation calculation, in principle, can resum all 
logs, but in practice the number of towers of logarithms included in the analytic Sudakov exponent 
depends on the level to which a fixed-order calculation was performed for a given process. Generally, 
if a NNLO calculation is available, then the B^'^^ coefficient (using the CSS formalism [44]) can be 
extracted and incorporated. If we try to interpret parton showering in the same language then we say 
that the Monte Carlo Sudakov exponent always contains a term analogous to A^^^ and B^^^ and that an 
approximation to A^'^^ is also present in some kinematic regions. 

In Ref. [43], predictions were made for both and Higgs production at the Tevatron and the 
LHC, using both resummation and parton shower Monte Carlo programs. In general, the shapes for the 
p± distributions agreed well, although the Pythia showering algorithm typically caused the Higgs p± 
distribution to peak at somewhat lower values of transverse momentum. 

'^Contributed by: J. Huston, I. Puljak, T. Sjostrand, E. Thome 
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8.3 Shower Algorithm Constraints 

While customarily classified as leading log, shower algorithms tend to contain elements that go beyond 
the conventional leading-log definition. Specifically, some emissions allowed by leading log are forbid- 
den in the shower description. Taking the Pythia [7, 10] initial-state shower algorithm [23,24,45] as an 
example, the following aspects may be noted (see [46] for further details): 

( i) Emissions are required to be angularly ordered, such that opening angles increase on the way in to the 
hard scattering subprocess. That is, non-angularly-ordered emissions are vetoed. 

(ii) The z and of a branching a ^ be are required to fulfill the condition u = Q'^ — — z) < 0. Here 
s = {Pa +PdY — {Pb + Pd)'^ / z, for d the incoming parton on the other side of the event. In the case that 
b and d form a Z°, say, and c is the recoiling parton, u coincides with the standard Mandelstam variable 
for the a -I- d (Z'^ = b + d) + c process. In general, it may be viewed as a kinematics consistency 
constraint. 

(Hi) The evolution rate is proportional to as((l — z)Q'^) f» as(pi) rather than as(Q^)- Since < 
this implies by itself a larger as and thus an increased rate of evolution. However, one function of the 
Qo ~ 1 GeV nonperturbative cutoff parameter is to avoid the divergent-ag region, so now one must 
require (1 — > Qq rather than > Qq. The net result again is a reduced emission rate. 

(iv) One of the partons of a branching may develop a timelike parton shower. The more off-shell this 
parton, the less the p± of the branching. The evolution rate in x is unaffected, however. 

(v) There are some further corrections, that in practice appear to have negligible influence: the non- 
generation of very soft gluons to avoid the divergence of the splitting kernel, the possibility of photon 
emission off quarks, and extra kinematical constraints when heavy quarks are produced. 

(vi) The emission rate is smoothly merged with the first-order matrix elements at large p±. This is some- 
what separate from the other issues studied, and the resulting change only appreciably affects a small 
fraction of the total cross section, so it will not be considered further here. 

The main consequence of the first three points is a lower rate of x evolution. That is, starting 
from a set of parton densities fiix, Qq) at some low Qq scale, and a matching A, tuned such that stan- 
dard DGLAP evolution provides a reasonable fit to data at > Qq, the constraints above lead to 
X distributions less evolved and thus harder than data. If we e.g. take the CTEQ5L tune [47] with 
A*^^) = 0.192 GeV, the A'^^^ would need to be raised to about 0.23 GeV in the shower to give the same fit 
to data as CTEQ5L when the angular-ordering cut in ( i) is imposed. Unfortunately effects from points 
(ii) and (Hi) turn out to be process-dependent, presumably reflecting kinematical differences between 
q — >^ qg and g gg. There is also some energy dependence. The net result of the first three points 
suggests that Pythia should be run with a A^^) of about 0.3 GeV (0.5 GeV) for Z° (H°) production in 
order to compensate for the restrictions on allowed branchings. 

One would expect the increased perturbative evolution to allow the primordial k± to be reduced. 
Unfortunately, while the total radiated transverse energy, J2i \P±i\' comes up by about 10% at the Teva- 
tron, this partly cancels in the vector sum, p±z = — J2i P±i- For a 2 GeV primordial k± the shift of the 
peak position of the p±z spectrum is negligible. Results are more visible for p^n at the LHC. 

Note that a primordial k± assigned to the initial parton at the low Qq scale is shared between 
the partons at each shower branching, in proportion to the longitudinal momentum fractions a daughter 
takes. Only a fraction Xhard/a^initiai of the initial k± thus survives to the hard-scattering parton. Since the 
typical X evolution range is much larger at the LHC than at the Tevatron, a tuning of the primordial k± 
is hardly an option for H*^ at the LHC, while it is relevant for Z° at the Tevatron. Therefore an increased 
A value is an interesting option. 

We now turn to the point f/vj above. By coherence arguments, the main chain of spacelike branch- 
ings sets the maximum virtuality for the emitted timeUke partons, i.e. the timelike branchings occur at 
longer timescales than the related spacelike ones. In a dipole-motivated language, one could therefore 
imagine that the recoil, when a parton acquires a timelike mass, is not taken by a spacelike parton but 
by other final-state colour-connected partons. A colour-singlet particle, like the Z° or H°, would then be 
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unaffected by the timelike showers. 

The consequences for p^z and p_lh of such a point of view can be studied by switching off timelike 
showers in Pythia, but there is then no possibility to fully simulate the recoiUng event. A new set of 
shower routines is in preparation [48], however, based on jo^-ordered emissions in a hybrid between 
conventional showers and the dipole approach. It is well suited for allowing final-state radiation at later 
times, leaving p^z and p±ii unaffected at that stage. Actually, without final-state radiation, the two 
approaches give surprisingly similar results overall. Both are lower in the peak region than the algorithm 
with final-state radiation, in better agreement with CDF data [49]. The new one is slightly lower, i.e. 
better relative to data, at small p±z values. 

A combined study [46], leaving both the primordial k± and the A value free, still gives some 
preference to {k±) = 2 GeV and the standard A(^) = 0.192 GeV, but differences relative to an alternative 
with {k±) = 0.6 GeV and A^^) = 0.22 GeV are not particularly large, Fig. 23. 
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Fig. 23: Comparison of the CDF p±z spectrum with the new shower algorithm for two parameter sets. 



8.4 Further Comparisons 

Returning to Higgs production at the LHC, in Fig. 24 are shown a number of predictions for the current 
standard PYTHIA shower routines. Using CTEQ5M rather than CTEQ5L results in more gluon radiation 
and a broader p± distribution due to the large value of A. Likewise turning off timelike showers for 
gluons radiated from the initial state also results in the peak of the p± distribution moving outwards. 

We can now compare the results with resummation descriptions and other generators, Fig. 25 [50]. 
As we see, the new Pythia routines agree better with resummation descriptions than in the past [43], 

attesting to the importance of various minor technical details of the Monte Carlo approach. One must 
note, however, that some spread remains, and that it is not currently possible to give an unambiguous 
prediction. 

8.5 Conclusions 

We have studied p±z and p±u spectra, as a way of exploring perturbative and nonperturbative effects in 
hadronic physics. Specifically, we have pointed out a number of ambiguities that can exist in a shower 
approach, e.g. that the shower goes beyond the simpleminded leading-log evolution and kinematics. 
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Fig. 24: Comparison of tiie Pythia p± distributions for Higgs production at tiie LHC using LO and NLO pdf 's, with and 
without (no FSR) timelike parton showering. 





Fig. 25: Comparison of various p± distributions for Higgs production at the LHC. The curves denoted Grazzini [51], ResBos 
[52], Kulesza [53] and Berger [54] are resummation descriptions, while MCNLO [55,56], Herwig [57] and PYTHIA are 
generators, PYTHIA 6.3 refers to the new algorithm outlined above. 
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while still making use of leading-log parton densities. Attempts to correct for mismatches in general tend 
to increase the perturbative p±^z- The need for an unseemly large primordial k± in the shower approach 
is thus reduced, but not eliminated. There is still room for, possibly even a need of, perturbative evolution 
beyond standard DGLAP at small virtuality scales. 

9. NEW SHOWERS WITH TRANSVERSE-MOMENTUM-ORDERING 

9.1 Introduction 

The initial- [23,24,45] and final-state [58,59] showers in the Pythia event generator [7, 10] are based 
on virtuahty-ordering, i.e. uses spacehke and timehke M^, respectively, as evolution variables. Other 
algorithms in common use are the angular-ordered ones in Herwig [57,60] and the -ordered dipole- 
based ones in Ariadne/Ldc [61, 62]. All three have been comparably successful, in terms of abil- 
ity to predict or describe data, and therefore have offered useful cross-checks. Some shortcomings of 
the virtuality-ordering approach, with respect to coherence conditions, have been compensated (espe- 
cially relative to HERWIG) by a better coverage of phase space and more efficient possibiUties to merge 
smoothly with first-order matrix elements. 

Recently, the possibility to combine matrix elements of several orders consistently with showers 

has been raised [63, 64], e.g. W + n jets, n = 0, 1, 2, 3, In such cases, a -ordering presumably 

offers the best chance to provide a sensible definition of hardness. It may also tie in better e.g. with the 
-ordered approach to multiple interactions [65]. This note therefore is a study of how the existing 
Pythia algorithms can be reformulated in pj^-ordered terms, while retaining their strong points. 

The main trick that will be employed is to pick formal definitions of pi^, that simply and unam- 
biguously can be translated into the older virtuality variables, e.g. for standard matrix-element merging. 
These definitions are based on lightcone kinematics, wherein a timelike branching into two massless 
daughters corresponds to p\ = z{l — z)NP and the branching of a massless mother into a spacelike and 
a massless daughter to p\ = {1 — z)Q'^. The actual p±_ of a branching will be different, and e.g. depend 
on the subsequent shower history, but should normally not deviate by much. 

9.2 Timelike Showers 

The new timelike algorithm is a hybrid between the traditional parton-shower and dipole-emission ap- 
proaches, in the sense that the branching process is associated with the evolution of a single parton, 
like in a shower, but recoil effects occur inside dipoles. That is, a dipole partner is assigned for each 
branching, and energy and momentum is 'borrowed' from this partner to give mass to the parton about to 
branch, while preserving the invariant mass of the dipole. (Thus four-momentum is not preserved locally 
for each parton branching a —>■ be. It was in the old algorithm, where the kinematics of a branching was 
not constructed before the off- or on-shell daughter masses had been found.) Often the two partners are 
colour-connected, i.e. the colour of one matches the anticolour of the other, as defined by the preceding 
showering history, but this need not be the case. In particular, intermediate resonances normally have 
masses that should be preserved by the shower, e.g., in t — > bW"*" the W"*" takes the recoil when the b 
radiates a gluon. 

The evolution variable is approximately the p\ of a branching, where p± is the transverse mo- 
mentum for each of the two daughters with respect to the direction of the mother, in the rest frame of 
the dipole. (The recoiling dipole partner does not obtain any p± kick in this frame; only its longitudinal 
momentum is affected.) For the simple case of massless radiating partons and small virtualities rela- 
tive to the kinematically possible ones, and in the limit that recoil effects from further emissions can be 
neglected, it agrees with the dij p_L-clustering distance defined in the PYCLUS algorithm [66]. 

All emissions are ordered in a single sequence p±max > P±i > P±2 > ■ ■ ■ > p±min- That is, 

"Contributed by: T. Sjostrand 
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each initial parton is evolved from the input p±max scale downwards, and a hypothetical branching p± 
is thereby found for it. The one with the largest p± is chosen to undergo the first actual branching. 
Thereafter, all partons now existing are evolved downwards from p±i, and a p±2 is chosen, and so on, 
until pi^min is reached. (Technically, the pj_ values for partons not directly or indirectly affected by a 
branching need not be reselected.) The evolution of a gluon is split in evolution on two separate sides, 
with half the branching kernel each, but with different kinematical constraints since the two dipoles have 
different masses. The evolution of a quark is also spUt, into one p± scale for gluon emission and one for 
photon one, in general corresponding to different dipoles. 

With the choices above, the evolution factorizes. That is, a set of successive calls, where the pxmin 
of one call becomes the p_Lmax of the next, gives the same result (on the average) as one single call for 
the full p± range. This is the key element to allow Sudakovs to be conveniently obtained from trial 
showers [64], and to veto emissions above some p± scale, as required to combine different n-parton 
configurations efficiently. 

The formal p± definition is p5_evoi ~ ^i^ " ^){^^ ~ ^-o)' where pxevoi is the evolution variable, 
z gives the energy sharing in the branching, as selected from the branching kernels, M is the off-shell 
mass of the branching parton and mo its on-shell value. This p±evoi is also used as ccg scale. 

When a p_Lcvoi has been selected, this is translated to a = ttiq + p^evoi/ (-^(^ ~ ^))- N*^*^^ ^^'^^ 
the Jacobian factor is trivial: dM^/ (M^ — rn^) dz = <ip'±cvoi/ P^cvoi Froni there on, the three-body 
kinematics of a branching is constructed as in the old routine. This includes the detailed interpretation of 
z and the related handling of nonzero on-shell masses for branching and recoiling partons, which leads 
to the physical p± not agreeing with the p^cvoi defined here. In this sense, p_Levoi becomes a formal 
variable, while M really is a well-defined mass of a parton. 

Also the corrections to b ^ bg branchings (6 being a generic coloured particle) by merging with 
first-order a beg mattix elements closely follows the existing machinery [59], once the p±evoi iias 
been converted to a mass of the branching parton. In general, the other parton c used to define the matrix 
element need not be the same as the recoiUng partner. To illustrate, consider a Z° — > qq decay. Say the 
q branches first, q —>■ qgi . Obviously the q then takes the recoil, and the new q, gi and q momenta are 
used to match to the — > qqg matrix element. The next time q branches, q qg2, the recoil is taken 
by the colour-coimected gi gluon, but the matrix element corrections are based on the newly created q 
and g2 momenta together with the q (not the gi !) momentum. That way one may expect to achieve the 
most realistic description of mass effects in the coUinear and soft regions. 

The shower inherits some further elements from the old algorithm, such as azimuthal anisotropics 
in gluon branchings from polarization effects. 

The relevant parameters will have to be retuned, since the shower is quite different from the old 
mass-ordered one. In particular, it appears that the five-flavour Aqcd value has to be reduced relative to 
the current default, roughly by a factor of two (from 0.29 to 0.14 GeV). 

9.3 Spacelike Showers 

Initial-state showers are constructed by backwards evolution [23], starting at the hard interaction and 
successively reconstructing preceding branchings. To simphfy the merging with first-order matrix ele- 
ments, z is defined by the ratio of s before and after an emission. For a massless parton branching into 
one spacelike with virtuality and one with mass m, this gives p^ = Q'^ — z{s + Q^){Q'^ + m?)/s, 
or p\ = {1 — z)Q'^ — zQ^/s for m = 0. Here s is the squared invariant mass after the emission, i.e. 
excluding the emitted on-mass-shell parton. 

The last term, zQ^/ s, while normally expected to be small, gives a nontrivial relationship between 
p\ and Q^, e.g. with two possible Q"^ solutions for a given p\. To avoid the resulting technical problems, 
the evolution variable is picked to be p^^cvoi ~ {1 — z)Q'^. Also here pxcvoi sets the scale for the running 
as- Once selected, the p\^^gi is translated into an actual by the inverse relation Q"^ = p\^-^^^i/ {I — z). 
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with trivial Jacobian: dQ"^ /Q"^ dz = dpievoi/^'levoi ^^'^^ the correct p^, including the zQ'^/s 
term, can be constructed. 

Emissions on the two incoming sides are interspersed to form a single falUng p± sequence, 

P±max > P±i > p±2 > • • • > p±miu- That is, the p± of the latest branching considered sets the 
starting scale of the downwards evolution on both sides, with the next branching occurring at the side 
that gives the largest such evolved p±. 

In a branching a be, the newly reconstructed mother a is assumed to have vanishing mass — 
a heavy quark would have to be virtual to exist inside a proton, so it makes no sense to put it on mass 
shell. The previous mother b, which used to be massless, now acquires the spacelike virtuality and 
the correct p± previously mentioned, and kinematics has to be adjusted accordingly. 

In the old algorithm, the 5 kinematics was not constructed until its spacelike virtuality had been 
set, and so four-momentum was explicitly conserved at each shower branching. In the new algorithm, 
this is no longer the case. (A corresponding change occurs between the old and new timelike showers, as 
noted above.) Instead it is the set of partons produced by this mother b and the current mother d on the 
other side of the event that collectively acquire the p± of the new a ^ be branching. Explicitly, when 
the b is pushed off-shell, the d four-momentum is modified accordingly, such that their invariant mass is 
retained. Thereafter a set of rotations and boosts of the whole b + d-produced system bring them to the 
frame where b has the desired p± and d is restored to its correct four-momentum. 

Matrix-element corrections can be applied to the first, i.e. hardest in p±, branching on both sides 
of the event, to improve the accuracy of the high-p^ description. Also several other aspects are directly 
inherited from the old algorithm. 

Work on the algorithm is ongoing. In particular, an optimal description of kinematics for massive 
quarks in the shower, i.e. c and b quarks, remains to be worked out. 

Some first tests of the algorithm are reported elsewhere [67]. In general, its behaviour appears 
rather similar to that of the old algorithm. 

9.4 Outlook 

The algorithms introduced above are still in a development stage. In particular, it remains to combine 
the two. One possibility would be to construct the spacehke shower first, thereby providing a hst of 
emitted partons with their respective emission p± scales. This list would then be used as input for the 
timelike shower, where each emission p± sets the upper evolution scale of the respective parton. This is 
straightforward, but does not allow a fully factorized evolution, i.e. it is not feasible to stop the evolution 
at some p± value and continue downwards from there in a subsequent call. The alternative would be to 
intersperse spacelike and timelike branchings, in one common p^-ordered sequence. 

Obviously the finished algorithms have to be compared with data, to understand how well they 
do. One should not expect any major upheavals, since checks show that they perform similarly to the old 
ones at current energies, but the hope is for a somewhat improved and more consistent description. The 
step thereafter would be to study specific processes, such as W + n jets, to find how good a matching can 
be obtained between the different n-jet multiplicities, when initial parton configurations are classified by 
their p^-clustering properties. The PYCLUS algorithm here needs to be extended to cluster also beam 
jets. Since one cannot expect a perfect match between generated and clustering-reconstructed shower 
histories, it may become necessary to allow trial showers and vetoed showers over some p± matching 
range, but hopefully then a rather small one. If successful, one may expect these new algorithms to 
become standard tools for LHC physics studies in the years to come. 
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10. MATCHING MATRIX ELEMENTS AND PARTON SHOWERS WITH HERWIG AND 
PYTHIA 



10.1 Introduction 

Parton-shower (PS) Monte Carlo event generators are an important tool in the experimental analyses 
of collider data. These computational programs are based on the differential cross sections for simple 
scattering processes (usually 2^2 particle scatterings) together with a PS simulation of additional 
QCD radiation that naturally connects to a model of hadronization. As the PS algorithms are based on 
resummation of the leading soft and coUinear logarithms, these programs may not reliably estimate the 
radiation of hard jets, which, in turn, may bias experimental analyses. 

Improvements have been developed to correct the emission of the hard partons in the PS. In 
PYTHIA [7, 10, 11], corrections were added for e+e" annihilation [68], deep inelastic scattering [69], 
heavy particle decays [59] and vector boson production in hadron-hadron colhsions [45]. Similarly, 
corrections were added to HERWIG [57,70] in Refs. [71-74] following the method described in [75]. 

The Catani-Krauss-Kuhn- Webber (CKKW) algorithm is a method for generalizing such correc- 
tions [63,76]. Along with this development, computer programs have become available [35,77] which 
are capable of efficiently generating multi-parton events in a format (the Les Houches format [78]) that 
can be readily interfaced with HERWIG and PYTHIA. Here, we report on how to use these programs 
combined with the HERWIG and PYTHIA Monte Carlo event generators to implement hard corrections 
to PS predictions. Several approaches are explored. One adheres closely to the CKKW algorithm, but 
uses HERWIG for adding an additional PS. The second is more closely tuned to the specific PS gener- 
ators themselves and calculates branching probabiUties numerically (using exact conservation of energy 
and momentum) instead of analytically. This is accompUshed by generating pseudo-showers starting 
from the various stages of a PS history. A comparison is also made with a much simpler method. 

10.2 Overview of the Correction Procedure 

PS's are used to relate the partons produced in a simple, hard interaction characterized by a large energy 
scale (large means » ^qcd) to the partons at an energy scale near Aqcd- At this lower scale, a tran- 
sition is made to a non-perturbative description of hadronic physics, with physical, long-lived particles 
as the final products. This is possible, because the fragmentation functions for the highly-virtual partons 
obey an evolution equation that can be solved analytically or numerically. This solution can be cast in 
the form of a Sudakov form factor, which can be suitably normaUzed as a probabiUty distribution for no 
parton emission between two scales. Using the Monte Carlo method, the evolution of a parton can be 
determined probabilistically, consisting of steps when the parton's scale decreases with no emission, fol- 
lowed by a branching into sub-partons, which themselves undergo the same evolution, but with a smaller 
starting point for the scale. The evolution is ended when the energy scale of parton reaches the hadroniza- 
tion scale ~ Aqcd- Starting from the initial (simple) hard process, a sampling of PS's generates many 
topologies of many-parton final states, subject to certain phase space and kinematic restrictions. How- 
ever, the evolution equation (as commonly used) only includes the soft and coUinear fragmentation that 
is logarithmically enhanced, so that non-singular contributions (in the limit of vanishing cut-offs) are 
ignored. This means that not enough gluons are emitted that are energetic and at a large angle from the 
shower initiator, since there is no associated soft or colhnear singularity. 

In contrast, matrix element (ME) calculations give a description of a specific parton topology, 
which is valid when the partons are energetic and well separated. Furthermore, it includes interference 
between amplitudes with the same external partons but different internal structure. However, for soft and 
coUinear kinematics, the description in terms of a fixed order of emissions is not valid, because it does 
not include the interference between multiple gluon emissions which cannot be resolved. 

The PS description of hard scattering would be improved if information from the ME were 
"^Contributed by: S. Mrenna 
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included when calculating emission probabilities. A systematic method for this can be developed 
by comparing the PS and ME predictions for a given fixed topology. Consider a hard scattering 
(e+e" — >^ 7/Z ^ qq) followed by a PS off the outgoing qq pair, with each branching i character- 
ized by a variable di The variables di represent some virtuality or energy scales that are evolved down to 
a cut-off dini. The PS rate for this given topology is a product of many factors: (1) the Bom level cross 
section for e+e~ qq, (2) Sudakov form factors representing the probability of no emission on each 
quark and gluon line, and (3) the branching factors at each vertex (or splitting). The ME prediction for 
this topology is somewhat more complicated. First, one needs to calculate the cross section for the full 
initial- and final-state (here e+e~ qqggq'q'). Then, one needs to specify a particular topology. There 
is no unique way to do this, but a sensible method is to choose a clustering scheme and construct a PS 
history. Ideally, the clustering variable would be the same as the virtuahty di used to generate the PS in 
the usual way. Having performed the clustering, one can then make a quantitative comparison of the two 
predictions. 

To facilitate the comparison, we first expand the PS prediction to the same fixed order in ag. This 
is equivalent to setting all the Sudakov form factors to unity. In this limit, we see that the PS product 
of the Bom level cross section and the vertex factors is an approximation to the exact ME prediction. 
As long as the values di are all large, the ME description is preferred theoretically, and the Sudakov 
form factors are indeed near unity. Therefore, the PS description can be improved by using the exact 
clustered ME prediction. When the values di are not all large, and there is a strong ordering of the value 
(di ^ d2 - ■ ■ ^ dini) then the PS description is preferred theoretically. In this limit, the ME prediction 
reduces to the product of Born level and vertex factors, provided that the argument of is chosen to 
match that used in the PS (this should be related to di). Therefore, the ME prediction can be used to 
improve the PS description in all kinematic regions provided that: (1) the correct argument for is 
used, and (2) the Sudakov form factors are inserted on all of the quark and gluon lines. This provides 
then an interpolation scheme between the PS and the ME prediction. As usual, there is a systematic 
uncertainty associated with how one chooses to perform the interpolation. 

This corrects the specific topology considered, but what of the rest of the topologies? ME calcula- 
tions can be performed for those that are simple enough, but technically there is a limitation. Presently, 
e+e~ 6 parton calculations can be performed using computational farms with appropriate cuts. A 
practical solution is to choose the cut-off dini large enough that the ME calculations in hand saturate the 
dominant part of the cross section. Then, an ordinary PS can be used to evolve the parton virtuaUties 
from dini down to the hadronization scale. It has been shown that the correct method for doing this 
consists of starting the PS 's at the scale where a parton was created in the clustered topology, and then 
vetoing branchings with virtualities larger than dini [63]. 

10.3 Results 

10.4 Discussion and Conclusions 

We have compared three different procedures of matching ME predictions with PS's using a methodology 
close to the CKKW algorithm suggested in [63,76]: (1) a slightly expanded version of the CKKW 
procedure using HERWIG as the PS generator (but not limited in principal to HERWIG) and exploiting 
the freedom to choose scales and cut-offs; (2) a version of the CKKW procedure relying on pseudo- 
showers and matched closely to the scales and cut-offs of PYTHIA and HERWIG; and (3) a much 
simpler procedure based on the approach of M. Mangano. Results are summarized in Figs. 1 (a) and (b). 

The HERWIG-CKKW procedure uses all of the elements of the original CKKW procedure, but 
expands upon them. Several choices of scale were investigated as starting points for the vetoed PS, and 
a wide range of prefactors were explored as arguments to the analytic NLL Sudakov form factor and as- 
The variation of the results with these choices is shown in the figures. Optimized choices were settled 
upon based on the smoothness of distributions, the agreement with HERWIG where expected, and the 
apparent improvement over the default HERWIG predictions. While this appears to be a tuning, the final 
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Fig. 26: (a) Differential fcri -cluster distributions da/dkn at the hadron level generated with the pseudo-shower procedure for 
pp +X collisions at = 1.96 TeV, for i = 1 — 5 and also showing the W"'' boson transverse momentum. The default 

result of PYTHIA is shown as a dashed line, while the result of the pseudo-shower algorithm is shown as a solid black line. The 
contribution to the pseudo-shower result from the two (red), three (green), four (blue), five (yellow) and six (magenta) parton 
components is also shown. The matching scale 10 GeV is shown as a vertical arrow; (b) Comparison of the ratio of various 
fcT-cluster distributions from HERWIG and PY using the pseudo-shower procedure, HERWIG using the MLM procedure, and 
HERWIG using the CKKW procedure for a matching scale of 15 GeV. 
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choices are easily justifiable. Since HERWIG is an angular-ordered shower, a variable such as fcy-cluster 
values is well suited as a starting point for the HERWIG shower. Because of the details of the HERWIG 
shower, a prefactor of ^ for the scale used in the Sudakov form factor is understandable, as well as a 
pref actor of | for the scale used in evaluating as- The results presented are better at the parton level than 
at the hadron level, which may require some tuning of the HERWIG hadronization model, and are less 
important at higher energies or when the cut-offs are larger. 

The pseudo-shower procedure uses the Sudakov form factors of HERWIG and PYTHIA to nu- 
merically calculate the Sudakov suppression. Since the Sudakov form factor is a probability distribution 
for no parton emissions, the suppression factor can be determined by starting showers from different 
stages of the PS history and discarding those events with emissions above a given cut-off. Because of 
the nature of this approach, there is less tuning of parameters. To match the argument used in as by 
default in HERWIG and PYTHIA, a different clustering scheme was used: pt clustering or LUCLUS- 
clustering. Final results at the hadron level are shown in the figures. In general, the hadron-level results 
are better than the parton-level ones. The use of LUCLUS over KTCLUS was driven by the kinematics of 
the PYTHIA shower. We have not checked whether KTCLUS works as well or better for the HERWIG 
results, and we leave this for future investigation. We should also investigate the advantages of using the 
exact clustering scheme of the individual generators: invariant mass and angular ordering for PYTHIA 
or just angular ordering for HERWIG. Also, since this work began, a new model of final-state showering 
was developed for PYTHIA which is exactly of the LUCLUS type. This should also be tested. 

The MLM procedure is a logical extension of the procedure developed by M. Mangano for adding 

PS's to VF-i-multijet events. It entails fcr-clustering the parton-level events, adding a PS (with HERWIG 
in practice, but not limited to it), and rejecting those events where the PS generates a harder emission 
(in the A;T-nieasure) than the original events. This approach yields a matching which is almost as good 
as the more complicated procedures based on the CKKW procedures explored in this work. The reason 
is not a pure numerical accident. The MLM procedure rejects events (equivalently, reweights them to 
zero weight) when the PS generates an emission harder than the lowest kx value of the given kinematic 
configuration. This is equivalent to the first step of the pseudo-shower procedure in the calculation of 
the Sudakov suppression when applied to the highest multiplicity ME. The remaining difference is in 
the treatment of the internal Sudakov form factors and the argument of a^. The agreement between the 
pseudo-shower and MLM procedures implies that the product of Sudakov form factors on internal lines 
with the factors of as evaluated at the clustering scale is numerically equivalent to the product of as 
factors evaluated at the hard scale. It is worth noting that, for the process at hand, qq' ^ W + X, only 
two of the cluster values can be very close to the cut-off, and thus only two of the as(kT) values can be 
very large. Also, at the matching scales considered in this study, 10 — 20 GeV, with a factorization scale 

on the order of Mw, Qf = \J Myy + P^^r, a fixed order expansion is of similar numerical reliability 
as the "all-orders" expansion of a resummation calculation. In fact, the resummation (PS) expansion is 
ideally suited for Q ^ Mw, whereas the fixed order expansion is best applied for Q ~ Mw- 

Based on the study of these three procedures, we can make several statements on the reliability of 
predicting the shapes and rates of multijet processes at collider energies. 

1 . The three matching procedures studied here can be recommended. They are robust to variation of 
the cut-off scale. 

2. The relative distributions in kr, for example, are reliably predicted. 

3. The variation in the relative distributions from the three procedures depends on the variable. For 
variables within the range of the ME's calculated, the variation is 20%. For variables outside this 
range, which depend on the truncation of the ME calculation, the variation is larger 50%. Of 
course, it is important to study the experimental observables to correctly judge the senstivity to the 
cut-off and methodology of matching. 

4. More study is needed to determine the best method for treating the highest multiplicity ME con- 
tributions. 
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5. The subject of matching is far from exhausted. The procedures presented here yield an improve- 
ment over previous matching prescriptions. However, these methodologies are an interpolation 
procedure. 
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11. W BOSON, DIRECT PHOTON AND TOP QUARK PRODUCTION: SOFT-GLUON COR- 
RECTIONS 15 

11.1 Introduction 

VF-boson, direct photon, and top quark production are all processes of considerable interest, useful in 
testing the Standard Model and searching for new physics. The hadroproduction cross sections for these 
processes have been calculated fully through next-to-leading order (NLO). Threshold corrections are 
known to be important in current hadron colliders and attempts have been made to calculate these soft- 
gluon corrections at next-to-next-to-leading order (NNLO) and beyond. Here I present results from the 
latest calculations. 

In general, at each order in perturbation theory the partonic cross section a for a hard-scattering 
process includes "plus" distributions with respect to akinematical variable, denoted say as S2, that mea- 
sures distances from partonic threshold. At ra-th order in the strong coupling ag (beyond the leading 
order) these distributions are of the type 



W(s2/M2 
S2 



Z < 2n - 1 , (8) 



where M is a hard scale, such as a mass or transverse momentum, relevant to the process at hand. 
These logarithmic terms are the soft-gluon corrections and they arise from incomplete cancellations near 
partonic threshold between graphs with real emission and virtual graphs. This is due to the limited phase 
space available for real gluon emission near partonic threshold. These threshold corrections, calculated in 
the eikonal approximation, can be formally shown to exponentiate [79-82] as a result of the factorization 
properties of the cross section. The logarithms with / = 2n — 1 are denoted as leading (LL), with 
Z = 2ra — 2 as next-to-leading (NLL), with Z = 2n — 3 as next-to-next-to-leading (NNLL), and with 
Z = 2n — 4 as next-to-next-to-next-to-leading (NNNLL). We note that the virtual corrections appear in 
S{s2) terms. A unified approach and master formulas for the calculation of these corrections for any 
process at NNLO have recently been presented in Ref. [83]. For the processes discussed in the next three 
sections, the LL, NLL, and NNLL terms have been calculated fully. In the NNNLL terms we have not 
included some process-dependent two-loop contributions [83] which, however, we expect to be small. 



11.2 W Boson Production 

The production of W bosons in hadron colliders is a process of relevance in testing the Standard Model, 
calculating backgrounds to new physics such as associated Higgs boson production, and luminosity 
monitoring. 

The calculation of the complete NLO cross section for W hadroproduction at large transverse mo- 
mentum was presented in Refs. [84—86]. The NLO results displayed an enhancement of the differential 

'^Contributed by: N. Kidonakis 
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Fig. 27: The differential cross section, da/dQ^, for W hadroproduction in pp collisions at the Tevatron at (left) y/S — 1.8 
TeV and (right) 1.96 TeV, with ^j.f = fiR ^ Qt- Shown are the Born, NLO, NNLO-NNLL, and NNLO-NNNLL results. 



distributions in transverse momentum Qt of the W boson. The Qt distribution falls rapidly with in- 
creasing Qt, spanning five orders of magnitude in the 30 GeV < Qt < 190 GeV region at the Tevatron. 

l^-boson production at high transverse momentum receives corrections from the emission of soft 
gluons from the partons in the process. The resummation and NNLO-NNLL corrections were studied 
in Ref. [87]. More recently the NNLO-NNNLL corrections were studied in Ref. [88]. These threshold 
corrections further enhance the cross section and reduce the scale dependence [88]. 

For the hadronic production of a high-Qj- W boson, with mass mw, the lowest-order partonic 
subprocesses are q{pa) + g{pb) — * W{Q) + q{pc) and q{pa) + q{pb) — > W{Q) + g{pc)- The partonic 
kinematical invariants in the process are s = {pa + Pb)^, t = {Pa — QY, u = {Pb — Q)^, which satisfy 
S2 = s + t + u — Q'^ = Oat partonic threshold. Here S2 = {pa + Pb — Q)"^ is the invariant mass of the 
system recoiling against the W boson and it parametrizes the inelasticity of the parton scattering. The 
partonic cross section a includes distributions with respect to S2 of the type ^ti {s2 / Q\) / 82]+- 

In Fig. 27 we plot the transverse momentum distribution, da jdQ\, for W hadroproduction at the 
Tevatron Run I with ^fS = 1.8 TeV and Run II with ^/S = 1.96 TeV. We use the MRST2002 NNLO 
parton densities [89]. We set the factorization scale ^xp and the renormalization scale equal to Qt- 
We focus on the high-Q^ region where the soft-gluon approximation holds well and the corrections 
are important. We see that the NLO corrections provide a significant enhancement of the Bom cross 
section. The NNLO-NNLL corrections provide a further modest enhancement of the Qt distribution. If 
we increase the accuracy by including the NNNLL contributions, which are negative, then we find that 
the NNLO-NNNLL cross section hes between the NLO and NNLO-NNLL results. 

The i^-factors are shown at ^/S = 1.96 TeV in the left frame of Fig. 28. We see that the i^-factors 
are moderate, and nearly constant over the Qt range shown even though the distributions themselves 
span two orders of magnitude in this range. It is also easy to see from the NLO/NLO-NLL curve that 
in the high Qt region the soft-gluon approximation holds very well, as the NLO-NLL cross section is 
almost identical to the full NLO result. 

The scale dependence of the differential cross section is shown on the right frame of Fig. 28 for 
Qrp = 80 GeV and = 1.8 TeV. We plot da/dQ^ versus h/Qt over two orders of magnitude: 
0.1 < /x/Qt < 10. We note the good stabilization of the cross section when the NLO corrections are 



43 




Fig. 28: Left: The iC'-factors for the differential cross section, da/dQ^, for W hadroproduction at the Tevatron with x/S = 
1.96 TeV and /if ~ fJ-R ~ Qt- Right: The fi dependence of the differential cross section, da/dQ^, for W hadroproduction 
at the Tevatron with ^fS = 1.8 TeV and Qt = 80 GeV. Here /i = = ^.r- 



included, and the further improvement when the NNLO-NNNLL corrections (which include all the soft 
and virtual NNLO scale-dependent terms) are added. 

11.3 Direct Photon Production 

Direct photon production is an important process for determinations of the gluon distribution function. 
The NLO cross section for direct photon production has been given in Refs. [90-92]. The role of higher- 
order soft-gluon corrections has also been addressed more recently. Threshold resummation and NNLO 
corrections for direct photon production have been presented in Refs. [93,94]. 

At lowest order, the parton-parton scattering subprocesses are q{pa) +g{Ph) lip-f) + q{pj) 
Q{Pa) + qiPb) liP-r) + g{pj)- We define the Mandelstam invariants s = ipa+ Phf, t = {pa - P-y)'^, 
and u = {ph — P-yY, which satisfy S4 = s + t + n = 0at threshold. Note that the photon transverse 
momentum is pT = (tu/s)^^"^. Here we calculate the cross section Ey d?a/d?py in the MS scheme. 
The soft corrections to the cross section appear in the form of plus distributions [W(s4/p|,)/s4]+. 

In order to show the effect of including the NNLO threshold terms, we start with a complete NLO 
calculation of the appropriate cross section using a program [95] which employs the phase-space slicing 
technique [96]. The original NLO calculation has been extended to include a complete NLO treatment 
of the bremsstrahlung contribution. The Set 2 fragmentation functions of [97] are used along with the 
CTEQ6M parton distribution functions [98]. In all cases the factorization and renormalization scales 
have been set equal to a common scale ^. Once the NLO results are obtained, the approximate NNLO 
contributions are then added to them. 

In the left frame of Fig. 29 a comparison is made to data from the E-706 Collaboration [99]. The 
NNLO-NNNLL curve at ^ = pj'/2 is practically indistinguishable from the NLO except at high pT- 
However, the ^ = 2pT NNLO result is much higher than NLO, and as a result the scale dependence 
at NNLO is considerably reduced. The theoretical band lies below the data at the lower end of the 
range covered by the data. In the right frame of Fig. 29 the rapidity dependence is shown for the UA-6 
proton proton data [100]. Again, the NNLO-NNNLL terms give a negligible contribution for the choice 
= /2 and the overall scale dependence is greatly reduced when the NNLO terms are added. 
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p p — > Y + X 




Fig. 29: NLO and NNLO-NNNLL results for direct photon production in hadronic collisions. Left: theory compared to data 
from the E-706 Collaboration [99] at pbeam = 530 GeV/c. Right: theory compared to pp data for the rapidity distribution 
from the UA-6 Collaboration [100] aiVS = 24.3 GeV. 



11.4 Top Quark Production 

Recent calculations for top hadroproduction include NNLO soft-gluon corrections to the double dif- 
ferential cross section [101-103] from threshold resummation techniques. The latest calculation [103] 
includes NNNLL and some virtual C terms (defined in Ref. [103]) at NNLO. When not all terms are 
known there is some difference between single-particle-inclusive (IPX) and pair-invariant-mass (PIM) 
kinematics. When NNNLL terms are included, the kinematics dependence of the cross section vanishes 
near threshold and is reduced away from it relative to NNLL accuracy. The factorization and renormal- 
ization scale dependence of the cross section is also greatly reduced. 

We study the partonic process ij ti with ij = qq and gg. In IPI kinematics, a single top quark 
is identified, i{pa) + j{Pb) — > t{pi ) + X\i] (^2) where t is the identified top quark of mass m and X\t] 
is the remaining final state that contains the t. We define the kinematical invariants s = {pa + Pb)^, 
ti = {Pb — PiY — fJ-i = {Pa — Pi)^ — "rn? and S4 = s + ti + ui. At threshold, 54 0, and the 
soft corrections appear as [In'(s4/m^)/s4]+. In PIM kinematics, we have instead i{jPa) + jipb) — > 
tt{p) + X{k). At partonic threshold, s = M^, with the pair mass squared. The soft corrections 
appear as [W(l - z)/{l - z)]+, with z = M^/s ^ 1 at threshold. 

In Fig. 30 we present the NLO and approximate NNLO ti cross sections at = 1.8 TeV 
(left frame) and 1.96 TeV (right frame) for fi = fip = /i^ = m. We use the MRST2002 NNLO 
parton densities [89]. The NNLO results include the soft NNNLL and virtual C terms in IPI and PIM 
kinematics. We also show the average of the two kinematics results which may perhaps be closer to the 
full NNLO result. 

In the left frame of Fig. 31 we show the scale dependence of the cross section, in the region 0.2 < 
^/m < 10, at = 1.96 TeV. The NLO cross section has a milder dependence on scale than the LO 
result. The NNLO cross section exhibits even less dependence on ///m, approaching the independence of 
scale corresponding to a true physical cross section. The change in the NNLO cross section in the range 
m/2 < ji < 2m, normally displayed as a measure of uncertainty from scale variation, is less than 3%. 
In the right frame of Fig. 31 we show the top quark transverse momentum distributions at = 1.96 
TeV. At NNLO we observe an enhancement of the NLO distribution with no significant change in shape. 
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Fig. 30: The ti total cross sections in pp collisions at ^/S =1.8 TeV (left frame) and 1.96 TeV (right frame) are shown as 
functions of m for n = m. The NLO (solid), and approximate NNLO IPI (dashed), PIM (dot-dashed) and average (dotted) 
results are plotted. 

Finally we note that recently [104] NNLO threshold corrections were calculated for top quark 
production via flavor-changing neutral-current (FCNC) processes at the Tevatron and HERA colliders. 
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12. EXTENDING THRESHOLD EXPONENTIATION BEYOND LOGARITHMS FOR DIS 
AND DRELL-YAN 

12.1 Introduction 

Threshold resummation [105, 106] sums terms in cross sections that grow as a production threshold is 
approached. There is however empirical [107], theoretical [105, 108-112], as well as speculative [113] 
evidence of the fact that the formalism enabling threshold resummation could be extended to include 
classes of terms that are either constant or decrease upon the approach of threshold. The resumma- 
tions of such terms, once put on a theoretically sound basis, would have significant phenomenological 
consequences, as shown in Ref. [107] for the case of Higgs production at hadron colliders. 

The first evidence that the dominant non-logarithmic perturbative contributions could be expo- 
nentiated goes back to [108], where it was shown that the partonic Drell-Yan cross section in the DIS 
factorization scheme contains the ratio of the timelike to the spacelike Sudakov form factor: large per- 
turbative contributions are left over in the exponentiated form of this ratio after the cancellation of IR 
divergences. This observation was made more precise in Ref. [105]. There, the resummation of thresh- 
old logarithms for the Drell-Yan process was proven to all logarithmic orders, making use of a procedure 

'^Contributed by: T.O. Eynck, E. Laenen, L. Magnea 
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Fig. 31: production in pp collisions at ^ = 1.96 TeV with m = 175 GeV and MRST20002 NNLO pdf 's [89]. Left: The 
scale dependence of the total cross section as a function of fj./m. The LO (dot-dot-dot-dashed), NLO (solid), and approximate 
NNLO IPI (dashed), PIM (dot-dashed) and average (dotted) results are shown. Right: The NLO and approximate NNLO top 
quark pr distributions are shown. 



of refactorization: the Mellin transform of the cross section is expressed near threshold, approached by 
letting the Mellin variable grow very large, as a product of functions, each organizing a class of in- 
frared and collinear enhancements; the refactorization is valid up to corrections which are suppressed by 
powers of at large N, so that terms independent of N can be treated by the same methods used to 
resum logarithms of A^. 

It has also been clear for some time [1 10] that at least a subset of the terms enhanced by logarithms 
but suppressed by a power of N can be resummed: in the MS scheme, for example, they arise from 
the exponentiation of the MS quark distribution. More recently, a factorization analysis of these terms 
was performed [111, 112], for the case of the longitudinal structure function in DIS, where however 
(log A^)/A^ terms are the leading ones. It would be of considerable interest to extend this analysis to 
other processes. Since these logarithms arise at one loop from finite remainders of collinear singularities, 
which are suppressed by infrared power counting, it is conceivable that joint resummation [114] might 
provide a framework for an all-order treatment. 

Here we report on recent work [1 15], where the results of Refs. [105, 109] were exploited to show 
that for processes which are electroweak at tree level (such as DIS and vector boson production through 
electroweak annihilation at colliders) the exponentation of A^-independent terms is in fact complete. 
The generalization of this result to processes with nontrivial color exchange is not straightforward, but 
it would of phenomenological interest for many processes to be studied at the LHC, including prompt 
photon, heavy quark and jet production. 



12.2 Extended exponentiation for electroweak annihilation 

Consider the Drell-Yan partonic cross section in the DIS scheme. The refactorization procedure de- 
scribed in Ref. [105] leads to the expression 

n ,.n_ (^(Q^6))^ f MN,e) y Ur{N,e) 1 
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The presence of the final-state jet function Jr{N, e) e.g. is direct consequence of dividing the unfactor- 
ized Drell-Yan cross section by the deep-inelastic structure function. A similar expression emerges in 
the MS scheme, 



(iV) 



( \nQ',e) 



+ 0{1/N) . 



(10) 



Notice the absence of Jr, which does not occur in the MS density. Here F is the quark form factor; 

(j), tp and X are different quark distributions, with the same collinear singularities but different finite 
contributions: (p, in particular, is the MS distribution, consisting only of collinear poles; U and V are 
eikonal functions describing soft radiation at wide angles; finally J is a jet function responsible for soft 
and colhnear contributions from the final state DIS current jet. All functions exponentiate, and in each 
case real (r) and virtual (v) contributions have been separated. Finite virtual contributions in both cases 
can be expressed just in terms of the quark form factor. Note that although each function is divergent and 
thus depends on e = 2 — d/2, the partonic cross sections are finite. 

All functions involved in Eqs. (9) and (10) can be precisely defined in terms of quark fields and 
eikonal lines. Renormalization group analysis and explicit evaluation lead to expressions that have the 
familiar exponential form [106], with corrections involving AT-independent terms. For the DIS scheme 
one finds 
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Similarly, for the MS scheme one has the expressioon 
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The functions ^dis and J^y^ are given to order below in Eqs. (15) and (16). As is well-known [105, 
106], the function A is responsible for leading logarithms of N to all orders. To achieve NLL accuracy 
one must compute it to two loops, obtaining 
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To NLL accuracy, one also needs the functions B and D to one loop, which are given by 



(13) 



(14) 



The remaining ingredients, collecting A?^-independent terms, are easily computed at one loop. For the 

DIS scheme one needs 



^(Q^e) 


2 r 


r(-g2,e) 


= exp 



as{Q) 



Cf(3C(2)) 



a 



J'Dis{as) = ^CF{- + a2) 



For the MS scheme, one further needs 
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Notice finally that by taking the ratio of Eq. (12) and Eq. (11) one finds directly the (square of) the DIS 
structure function (ms)(-^)' factorized in the MS scheme, which can then be computed to the same 
accuracy without introducing other information. The ratio of form factors drops out from F^ (ms)(-^)' 
as does the function D [116], except for a contribution to the running between the physical scales 
and (1 — which can be systematically reexpressed as a modification of the function B. 

All the functions appearing in Eq. (11) and Eq. (12) can be exphcitly evaluated at two loops by 
matching with the complete two-loop calculation of Ref. [1 17]. An alternative, and often simpler method 
to determine the two-loop coefficients in these expressions uses equations derived from the fact that real 
and virtual contributions in the factorized expressions (16) and (10) are separately finite [115]. 

It is also important to keep in mind that the exponentiation of A^-independent terms does not have 
the predictive power of the standard resummation of threshold logarithms. In that case, typically, an 
entire tower of logarithms can be exactly predicted to all orders by performing just a low order calcu- 
lation. Here, on the other hand, functions such as Tdis ^^jg receive new nontrivial contribution 
at each perturbative order. The exponentiation pattern is nonetheless nontrivial, and higher order terms 
predicted by the exponentiation can be considered representative of the size of the complete higher order 
correction [115]. 

13. JOINT RESUMMATION FOR TOP QUARK PRODUCTION 

The formalism of joint resummation [114, 118] for hadronic cross sections of distributions singular at 
partonic threshold and at zero recoil has so far been applied to only a few processes. Recent studies 
involve processes that proceed at lowest order through a 2 ^ 1 electroweak {Z/W production [1 19]) or 
Yukawa interaction (Higgs production [53]). For these cases, the observables are the production cross 
sections at fixed mass M and measured Qt- Partonic threshold is then defined hy z = M'^ /s = 1, where 
s is the partonic center of mass (cms) energy squared, and zero recoil by Qt = 0. At any finite order, the 
disfi-ibutions take the form of plus-disfi-ibutions [ln''(l - z)/{l - z)]^ and \hi^ {M / Qt) / Qt\ j^- Note 
that in these observables the latter distributions enter the physical cross sections, whereas the former are 
defined, after factorization, in the context of a perturbative analysis of the hard scattering. 

In Ref. [120] the case of the prompt photon hadroproduction cross section at measured pt was 
analyzed, and a preliminary numerical study performed. In this case, a single-particle inclusive process 
proceeding through a 2 ^ 2 reaction at lowest order, the identification of the recoil variable Qt can 
only be made in the context of a refactorization analysis, just like the threshold variable z. Through 
such an analysis, it is possible to identify a reduced hard scattering with cms energy squared s and at 
transverse momentum Qt with respect to the hadronic cms system. Note that this transverse momentum 
is invariant w.r.t. longitudinal boosts. In this way, the hard scattering need produce a photon with 
transverse momentum pr' = pr — Qt/"^- It still remains to implement a procedure that consistently 
matches the joint resummed cross section to finite order calculations. In Refs. [114, 120] a simple cut-off 
p, was used in the integration over Qt, with p, smaller that 2pT. For the prompt photon case, the hard 
scattering is singular at Qt = 2pT. 

In this brief report we present a preUminary study of the application of joint resummation to 
another prominent single-particle inclusive cross section, the pT distribution of top quarks produced in 
hadronic collisions. Our motivation is to see what effect joint resummation has on a distribution in a TeV 
collider process that is nevertheless near threshold. Two key differences with the prompt photon case 
are ( i) the heavy quark mass m, preventing a singularity in the reduced hard scattering function when 
Qt = 2pT, and (ii) the possibility of multiple colored states for the produced top quark pair. 

"Contributed by: A. Banfl, E. Laenen 
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The jointly resummed expression for this observable can be written as 



[d^QreiH-QT)^^^^^^, (17) 
dpT J dprd^Qr 



where 



X 



g-2CFt(iV) (RcL^+1) 



7V+1 (18) 



where the impact vector b is Fourier conjugate to Qt, and the variable N is Laplace conjugate to the 
variable 

1 _ 4 = 1 _ _T (19) 
where x\ = 1 defines hadronic threshold. Furthermore 



(20) 



where /3 = ^/Y^^4rn?Js. The exponentional function is given in Ref. [119] and Ml{N),M^{N) 
are the Laplace moments of the lowest order heavy quark production matrix elements for either the qq or 
gg channel, the index labelUng the color-state of the heavy quark pair. 

In two figures we illustrate the effect of joint resummation over threshold resummation for the 

top quark (m = 175 GeV) px spectrum for at Run 11 Tevatron, and only for the dominant qq channel. 
It produces the top quark pair at lowest order in an octet state. We use ag = 0.1 and toy densities 
'Pa{x) = x{l — x)^. In Fig. 32 we exhibit the Qt profile of Eq. (18) for two rather large pt values, 




Fig. 32: The Qt profile for two pr values. 



in analogy to Fig. 1 in Ref. [120]. Note that these profiles are only of theoretical relevance, only their 
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integral over Qt can be measured. Recoil effects can be removed by neglecting Qt in the last factor 
in Eq. (18). We observe a small enhancement over the threshold-resummed result, in particular at very 
large pt- Although not shown in the figure, there is in our case no divergence when Qt approaches 2pT. 

To keep the recoil soft with respect to the hard scale, we choose the cut-off fi = 0.2mT- We can 
now show, in Fig. 33 the effect of joint resummation on the pT distribution, via Eq. (17). We observe 
only a small enhancement at very large pt values. The suppression at lower pt values is a consequence 
of the Qt cut fi. A proper matching procedure should resolve this issue. Actually we note that this 
suppression is absent if we choose /Z around 200 GeV. The enhancement at large pt is a consequence of 




Pj [GeV] [GeV] 



Fig. 33: The pr spectra for top quark producation at Tevatron run II, using toy PDF's a;(l — a;)^, in the qq channel, and their 
ratio. 

what is effectively a smearing of the recoilless pT spectrum induced by perturbative radiation. 

A more extenstive study including the gg channel and use of realistic parton distribution functions 
is forthcoming [121]. 
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14. A COMPARISON OF PREDICTIONS FOR SM HIGGS BOSON PRODUCTION AT THE 

LHC i« 

14.1 Introduction 

The dominant mechanism for the production of a SM Higgs boson at the LHC is gluon-gluon fusion 
through a heavy (top) quark loop. For this reason this channel has attracted a large amount of theoretical 
attention [122]. Recently, the total cross section has been calculated to NNLO in the strong coupling con- 
stant as (i.e. at order ©(af )) [123-127] and also contributions from multiple soft gluon emission have 
been consistently included to NNLL accuracy [128]. In addition to the size of the total rate, a knowledge 
of the shape of the Higgs boson pT distribution is essential for any search and analysis strategies at the 
LHC. In particular, the pt distribution for the Higgs boson is expected to be harder than the one of its 

'^Contributed by: C. Balazs, M. Grazzini, J. Huston, A. Kulesza, I. Puljak 
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corresponding backgrounds. The Higgs boson pT distribution has been computed with LL parton shower 
Monte Carlos (HERWIG [70] and PYTHIA [7]), and through various resummed calculations. The latter 
techniques are the more powerful ones, but it is primarily the former that experimentahsts at the LHC 
have to rely upon, because of their flexibility in allowing to test the effects of the various kinematic cuts 
which may optimize search strategies. 

In the kinematic region p\ <C m\, where most of the events are expected, large logarithmic 
corrections appear of the form a" In"* m?^ Ip\ that spoil the vaUdity of the fixed order perturbative 
expansion. The pT distribution can be written as 



du da^^- da 



■fin. 



dp^ dp\ dp\ ^ ^ 

The first term contains all logarithmically-enhanced contributions and requires their resummation to all 
orders. The second term is free from logarithmically-enhanced contributions and can be evaluated at 
fixed order in perturbation theory. The method to perform the all-order resummation is well known: to 
correctly take into account momentum conservation the resummation must be performed in the impact 
parameter (6) space [129, 130]. The large logarithmic contributions are exponentiated in the Sudakov 
form factor, which in the CSS [131] approach takes the form 



bi/b'^ 



A, {as{ix)) In 1^^^ + Be 



(22) 



where = 2e and c = q,g. The Ac and Be functions are free of large logarithmic corrections and 
can be computed as expansions in the strong coupling constant a^: 



oo 



Aeias)= ^^^^ 

n=l 
oo 

Be{as)= J2{^yBjrl (24) 

n=l 

The functions Ac and Be control soft and flavour-conserving colUnear radiation at scales mn. 
Purely soft radiation at a very low scales fi^l/b cancels out because the cross section is infrared safe 
and only purely coUinear radiation up a scale /x ~ 1/6 remains, which is taken into account by the 
coefficients 

oo 

Cab{as,z)=^{^yC^^{as,z). (25) 

n=l 

Beyond NLL accuracy, to preserve the process independence of the resummation formula, an additional 
(process dependent) coefficient H is needed [132], which accounts for hard virtual corrections and has 
an expansion 

oo 

Hc{as) = l + Y.{^yH(^l (26) 

n=l 

In the case of Higgs boson production through gg fusion, the relevant coefficients A^p, A^p and Bg^^ 
are known [133] and control the resummation up to NLL accuracy The NNLL coefficients C^^ and 
Hg^^ are also known [132, 135]. The NNLL coefficient Bg^^ has been computed in Refs. [136, 137], 
whereas A^g^ is not yet known exactly. In the following we assume that its value is the same that appears 
in threshold resummation [138]. 

"There are two different classification sciiemes of the LL, NLL, NNLL, etc terms and their corresponding B contents. Here 
we use the most popular scheme. Another is discussed in Ref. [134]. 
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14.2 Predictions for pt spectra and comparisons 

In the 1999 Les Houches workshop, a comparison [43, 139] of the HERWIG and PYTHIA (2 versions) 
predictions for the Higgs boson px distribution with those of a pr resummation program (ResBos [52, 
140]) was carried out. This comparison was continued in the 2001 workshop and examined the impact 
of the B^'^^ coefficient [122]. In the meantime, a number of new theoretical predictions have become 
available, both from resummation and from the interface of NLO calculations with parton shower Monte 
Carlos. For these proceedings, we have carried out a comparison of most of the available predictions for 
the Higgs boson pT distribution at the LHC. We have used a Higgs boson mass of 125 GeV and either 
the MRST2001 or the CTEQ5M pdf 's. The difference between the two pdf 's for the production of a 125 
GeV mass Higgs boson is of the order of a few percent. Before comparing the different predictions, we 
comment on the various approaches in turn. 

Parton shower MC programs such as HERWIG, which implements angular ordering exactly, im- 
plicitly include the A^'^\ A^"^^ and B^^^ coefficients and thus correctly sum the LL and part of the N'^LL 
contributions. However, in the most straightforward implementations, MC cannot correctly treat hard 
radiation. By contrast, the PYTHIA MC, which does not provide an exact implementation of angular 
ordering, has a hard matrix element correction Recently, an approach to match NLO calculations to 
parton showers generators, MC@NLO [55,56], has been proposed, and applied, amongst the other, to 
Higgs production. This method joins the virtues of NLO parton level generators (correct treatment of 
hard radiation, exact NLO normalization) to the ones of MC. It thus can be compared to a resummed 
calculation at NLL+NLO accuracy. 

As far as resummed calculations are concerned, we first consider two implementation of the CSS 
approach. The ResBos code includes the Ac , Be and C^iJ coefficients in the low-p^ region and 
matches this to the NLO distribution at high pT- NNLO effects at high pT are approximately taken into 
account by scaling the second term in Eq. (21) with a K-f actor. The matching is performed through a 
switching procedure whose uncertainty will be considered in the following. The calculation of Berger 
and Qiu [54] also performs a pT resummation in b space and is accurate to NLL. The coefficient B^^^ 
is included but the matching is still to NLO. Note that in both these approaches the integral of the 
spectrum is affected by higher-order contributions included in a non-systematic manner whose effect is 
not negligible for Higgs production. 

The prediction by Bozzi, Catani, de Florian and Grazzini [51] (labeled Grazzini et al. in the 

following) is based on an implementation of the 6-space formalism described in [51,132]. The calculation 
has the highest nominal accuracy since it matches NNLL resummation at small px to the NNLO result 

(2) (2) 

at high Pt [142]. This approach includes the coefficients C^^^ and Hg in approximated form. The 
main differences with respect to the standard CSS approach are the following. A unitarity constraint is 
imposed, such that the total cross section at the nominal (NNLO) accuracy is exactly recovered upon 
integration. A study of uncertainties from missing higher order contributions can be performed as it is 
normally done in fixed order calculations, that is, by varying renormalization and factorization scales 
around the central value, that is chosen to be mn- 

Finally, we discuss the px distribution of Ref. [53] (Kulesza et al.). This is obtained using a joint 
resummation formalism, by which both threshold and low-p^ logarithmic contributions are resummed 
to all orders. This approach has been formally developed to NLL accuracy, but the NNLL coefficients 

can also be incorporated. The matching is still performed to NLO. Even 
though a low mass Higgs boson at the LHC is produced with relatively low x partons, threshold effects 
can still be significant due to the large color charge in the gg initial state as well as steep x dependence 
of the gluon distribution functions at low x. This leads to an increased sensitivity to Sudakov logarithms 
associated with partonic threshold for gluon-induced processes, as shown in Ref. [128]. 

It is known that the low-pT region is sensitive to non-perturbative effects. These are expected 
^"Very recently hard matrix element corrections for Higgs productions have been implemented in HERWIG as well [141]. 
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to be less important in the gluon channel due to the larger colour charge of the gg initial state [43]. 
Different treatments of non-perturbative effects are included in the ResBos, Berger et al. and Kulesza et 
al calculation, wherei 
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Fig. 34: The absolute predictions for the production of a 125 GeV mass Higgs boson at the LHC. 

The absolute predictions for the cross sections are shown in Figure 34. All curves are obtained 
in the m-top oo limit. HERWIG and PYTHIA cross sections are significantly smaller than the other 
predictions, their normalization being fixed to LO. In the high-p^ region, the HERWIG prediction drops 
quickly due to the lack of hard matrix element corrections. PYTHIA, in contrast, features the hard matrix 
element corrections. We also note that PYTHIA prediction is significantly softer than all the other curves, 
and thus its overall shape is fairly different from all the other predictions. 

The MC@NLO cross section, about 32.4 pb, is roughly twice that of the HERWIG and PYTHIA 
predictions, being fixed to the NLO total cross section. 

Two predictions {step, smooth) are shown for ResBos which differ in the manner in which the 
matching at high pT is performed. Their difference can be considered as an estimate of the ambiguity 
in the switching procedure. The two curves correspond to the same total cross section of about 36.2 
pb, which is about 8 % higher than the NLO cross section. This is the effect of the higher-order terms 
that enter the prediction for the total rate in the context of the CSS approach. A slightly softer curve is 
obtained by Berger and Qiu. The predicted cross section (37 pb) is close to that of ResBos. 

The Grazzini et al. prediction has an integral of about 39.4 pb, which corresponds to the total 
cross section at NNLO. Contrary to what is done in Ref. [51], here the curve is obtained with MRST2002 
NNLO partons and three-loop a^- The difference with the result obtained with MRST2001 NNLO PDFs 
is completely negligible. 

Concerning the Kulesza et al curve, the subleading terms associated with low x emission (i.e. in 
the limit opposite to partonic threshold) and of which only a subset is included in the joint resummation 
formalism, play an important role numerically. As a result, the total cross section turns out to be 35 pb, 
about 10% lower than the pure threshold result, which is 39.4 pb [53]. 
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We now want to examine in more detail the relative shaoes of the oredictions clotted in Figure. 34. 
In Figure. 35 all the p 9.4 pb. 
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Fig. 35: The predictions for the production of a 125 GeV mass Higgs boson at the LHC, all normalized to the same cross 
section for better shape comparison. 



In the region of small and moderate pT (say, smaller than 100 GeV) all of the predictions are 
basically consistent with each other, with the notable exception of PYTHIA, which predicts a much 
softer spectrum. The curve of Kulesza et al. is also softer than the others. 

For larger px, HERWIG gives unreliable predictions, since the transverse momentum is generated 
solely by means of the parton shower, and therefore it lacks hard matrix element effects. The Grazzini et 
al. and ResBos curves are harder than MC@NLO for large px- There are two reasons for this. Grazzini 
et al. implement the NNLO matrix elements exactly, corresponding to the emission of two real partons 
accompanying the Higgs in the final state [142]; ResBos mimics these contributions, by multiplying the 
NLO matrix elements by the K factor. MC@NLO, on the other hand, contains only NLO matrix elements 
(one real parton in the final state). Secondly, Grazzini et al. and ResBos choose the renormalization and 
factorization scales equal to mn, whereas in MC@NLO these scales are set equal to the transverse mass 

of the Higgs, ^ m?^ + p^. The difference is small at the level of total rates, but it is not negligible in the 
tail of the pT distribution. 

14.3 Conclusions 

Up to now, the ATLAS and CMS experiments have relied primarily on the predictions from HERWIG 
and PYTHIA in designing both their experiments as well as defining their search and analysis strate- 
gies. In the last few years, a number of tools for and predictions of the Higgs boson cross section at the 
LHC have become available, with the inclusion of beyond-the-leading-order effects at different level of 
accuracy. In the case of total rates, NNLO results have recently become available; their consistent inclu- 
sion in experimental analysis will allow to further decrease the estimated lower bound on the integrated 
luminosity to be collected for discovery. 
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In this contribution, we primarily focused on the predictions for the pT spectrum, comparing the 
results of Monte Carlos with those obtained with analytically-resummed calculations. In contrast to 
the situation in 1999, all of the predictions, with the exception of PYTHIA, result in the same general 
features, most notably in the position of the peak. However, differences do arise, because of different 
treatments of the higher orders. It is an interesting question beyond the scope of this review that of 
whether these differences are resolvable at the experimental level, which may lead to modify the strategy 
for searches. In order to answer this, studies including realistic experimental cuts must be performed 
with the newly available tools. 

15. MATRIX-ELEMENT CORRECTIONS TO gg/qq HIGGS IN HERWIG 

15.1 The Higgs transverse momentum 

In order to investigate Higgs boson production via gg Higgs (see Ref. [143]), one needs to account for 
multi-parton radiation for the sake of performing trustworthy phenomenological analyses [43, 122, 144]. 
Standard Monte Carlo (MC) algorithms [7,57,70] describe parton radiation in the soft and/or collinear 
approximation of the parton shower (PS), but can have regions of phase space, so-called 'dead zones', 
where no radiation is allowed. Here, one can however rely on higher-order tree-level results, as in 
this region the radiation is neither softly nor colUnearly enhanced. Several methods have been recently 
suggested in order to match PS and fixed-order matrix elements (MEs) [59,75], also including the virtual 
one-loop terms [55,56, 145]. 

15.2 The HERWIG implementation 

In this note, we briefly mention that the same strategy which has already been used to implement real ME 
corrections to e+e" annihilation into quark pairs [71], Deep Inelastic Scattering (DIS) [72], top quark 
decay [73] and vector boson hadro-production [74] has now also been adopted for the case of Higgs 
hadro-production via gluon-gluon fusion, in the context of the HERWIG event generator [57,70]. That 
is, the dead zone is here populated by using the exact next-to-leading order (NLO) tree-level ME result 
and the PS in the already-populated region is corrected using the exact amplitude any time an emission 
is capable of being the hardest so far. 

15.3 Numerical results and comparisons 

The MEs squared for the real corrections to gg H that we have used can be found in [146], where top 
mass effects are fully included. The real NLO corrections to qq ^ H are instead rather straightforward: 
the formulae we used can be read from Eq. (3.62) of [147] with appropriate Yukawa couplings and 
crossing. In the new HERWIG default version, in line with [74], ME corrections use the Higgs transverse 
mass my = q^ + as the scale for as and for the Parton Distribution Functions (PDFs) while the 
99 J QQ ^ H contributions use mjj. We shall also assume that the intrinsic transverse momentum of the 
initial-state partons is equal to qT,mt = 0> the HERWIG default value. 

By adopting the HERWIG defaults, we first consider Higgs production at the Tevatron and the 
LHC within the MC itself, by plotting the qx distribution with (solid histogram) and without (dotted) 
ME corrections: see Fig. 36. Beyond gr ~ ran/ 2 the ME-corrected version allows for many more 
events. In fact, one can prove that, within the standard algorithm, qt is constrained to be qx < tuh- 
At small qT the prediction which includes ME corrections displays a suppression. By default, after the 
latter are put in place, the total normalization still equals the LO rates. Hence, it is obvious that the 
enhancement at large qT implies a reduction of the number of events which are generated at small qT 
values. 

^'Contributed by: G. Corcella, S. Moretti 
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Fig. 36: Higgs transverse momentum distribution according to HERWIG with (solid) and without (dotted) ME 
corrections, at Tevatron (left, ^fspp = 2 TeV) and LHC (right, Vs^p = 14 TeV). We have set the Higgs mass to 
niH = 115 GeV. 
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Fig. 37: Left: comparison of ME-corrected HERWIG predictions (solid) to the + jets' result from [146] (dotted). 
Centre: comparison of ME-corrected HERWIG predictions (solid) to the NLO and resummed calculation of [51] 
(dotted). Right: comparison of ME-corrected HERWIG predictions (solid) to the MC@NLO results from the code 
described in Ref. [148] (dotted). Here, qq ^ H processes have been turned off. 

In Fig. 37 (left plot) we present the improved HERWIG spectrum (solid) for the LHC, along with 
the result obtained running the so-called -i- jets' process (dotted), where the hard process is always one 
of the corrections to gg H. In order to perform such a comparison, we have turned the qq ^ H hard 
process off, as + jets' in HERWIG does not currently implement the corrections to quark-antiquark 
annihilation. Furthermore, we have chosen gTmin = 30 GeV for the 'H -i- jets' generation. As expected, 
at small qr the two predictions are fairly different but at large transverse momentum they agree well. 

In Fig. 37 (centre plot) we compare the new HERWIG version with the resummed calculation 
of Ref. [51]. For the sake of comparison with HERWIG, which includes leading logarithms and only 
some subleading terms, we use the results of [51] in the NLL approximation (rather than the default 
NNLL one), matched to the NLO prediction. In order for such a comparison to be trustworthy, we 
have to make parameter choices similar to [51]: namely, we adopt a top quark with infinite mass in 
the loop and tuh = 125 GeV, with as and PDFs (both from HERWIG defaults) evaluated at m|^. 
While the normalization (LO in HERWIG, NLO in Ref. [51]) and the small-qr behaviour of the two 
curves are clearly different, the large transverse momentum predictions are in good agreement, as in both 
approaches it is the real NLO ME that dominates the event generation at large qr- 

Finally, in Fig. 37 (right plot), we compare the results of standard HERWIG after ME corrections 
with the MC@NLO event generator (version 2.2) of Ref. [148], the latter implementing both real and 
virtual corrections to the hard-scattering process, in such a way that predicted observables (including nor- 
malization) are correct to NLO accuracy. As version 2.2 of the MC@NLO includes only the corrections 
to Higgs production in the gluon-fusion channel, we again have turned the quark-annihilation process off 
in our routines. As observed in the comparison with the resummed calculation, the two spectra differ in 
normalization and at small qr, but agree in the large-transverse-momentum region. 
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15.4 Conclusions 



Between the described implementation and the one available within the MC@NLO option, we believe 
that HERWIG is presently a reliable event generator for (direct) Higgs production from parton fusion 
at hadron colliders both at small and large transverse momentum. In fact, all currently available ME 
corrections will play an important role to perform any analysis on Higgs searches at present and future 
colliders. In particular, the option described here may be the most convenient choice for when the phase 
space is limited to transverse momentum values such that qr ^ mn- 
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16. CAESAR: AUTOMATING FINAL-STATE RESUMMATIONS ^2 

Event shapes and jet resolution parameters (final-state 'observables') measure the extent to which the 
energy-flow of the final state departs from that of a Born event. Their study has been fundamental for 
measurements of the strong coupling [149, 150] as well as the QCD colour factors [151]; final states 
also provide valuable information on the yet poorly understood transition from parton to hadron level 
(see [152] for a recent review). In the region where an observable's value v is small, one should resum 
logarithmically enhanced contributions that arise at all orders in the perturbative series. For a number 
of observables such a resummation has been carried out manually at next-to-leading logarithmic (NLL) 
accuracy [153]. But achieving NLL accuracy requires a detailed analysis of the observable's properties, 
and is often technically involved. We have instead recently proposed [154] a new approach based on a 
general NLL resummed master formula valid for a large class of final-state observables (which includes 
many of the final-state observables for which a NLL resummation already exists). We have also formu- 
lated the formal requirements that an observable should satisfy so as to be within the scope of the master 
formula. 

To illustrate these requirements we consider a Bom event consisting of n hard partons or 'legs'. 
For an observable (a function V of all final-state momenta) to be resummed in the n-jet limit it should: 

1 . vanish smoothly after addition of an extra soft particle colUnear to any leg i, with the following 



Here Q is a hard scale of the problem; {p} represents the Bom (hard) momenta after recoil from 
the emission, which is defined in terms of its transverse momentum kt and rapidity rj with respect 
to leg £, and where relevant, by an azimuthal angle relative to a Born event plane. 

2. be recursively infrared and coUinear (rlRC) safe: meaning roughly that, given an ensemble of 
emissions, the addition of a relatively much softer or more coUinear emission does not significantly 
alter the value of the observable, no matter how soft and/or coUinear the other emissions are. 
This condition is necessary to ensure exponentiation of leading logarithms, and is not satisfied for 
example by the JADE algorithm three-jet resolution parameter [155]. 

3. be continuously global, meaning that the observable is sensitive to emissions in the whole of the 
phase space ('global') in a continuous way, the formal requirement being d^jdinkt 1^ ^{{p}i k) =0 
and d(j)dinktl^V({p}, k) = (kt being defined with respect to the nearest leg). For non- 
continuously-global observables one must account for non-global logarithms (known only for large 
Nc) [156]. Their evaluation for a general observable would necessitate a (quite non-trivial) deter- 
mination of the phase-space boundaries associated with the observable. 

^Contributed by: A. Banfl, G.P. Salam, G. Zanderighi 
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Given the above conditions, the NLL resummation for the observable's distribution (the probabil- 
ity S(f ) that the observable's value is less than v) for a fixed Bom configuration is given by the 'master' 
formula [154]: 

S(^) = e-^('')^(i?'(7;)) , R'{v) = -v^^ ■ (28) 

The function R{v) is a Sudakov exponent that contains all leading (double) logarithms and all NLL 
(single-log) terms that can be taken into account by exponentiating the contribution to T,{v) from a 
single emission. This function depends parametrically on a^, 6^, and on the azimuthal average (In 5^); 
its full expression is reported in [154]. Multiple emission effects, for example the fact that even if all 
V{{p}, ki) < V, one might nevertheless have V{{p}, fci, . . . , kn) > v, are accounted for by the NLL 
function T, which can be computed via a Monte Carlo procedure [157]. 

The advantage of having introduced a master formula is that the resummation of the observable 
can be performed entirely automatically. The master formula and applicability conditions are encoded in 
a computer program (CAESAR, Computer Automated Expert Semi-Analytical Resummer), which given 
only the observable's definition in the form of a computer routine, returns the observable's distribution 
S(v) at NLL accuracy (where possible). 

As an example we present explicit results for the specific case of the (global) transverse thrust in 
hadronic dijet production, defined as 

r,^l-ma. ^f^'-"^^l , (29) 

where the sums run over all particles in the final state, the p±i are the particle transverse momenta (with 
respect to the beam direction) and n± is a unit transverse vector. The program, probing the observable 
with randomly chosen soft and collinear emissions, is able to verify that the applicability conditions hold 
and to determine the parameters ae, be, di, as well as the function gi{(l)). It then apphes the Monte Carlo 
procedure introduced in [157] to compute the function T. The results from this analysis are then plugged 
into the master formula (28) to compute ^{t±) at NLL accuracy. The resulting differential distribution 
D{t±) = d'E{Tj_)/dlnT± (integrated over a range of Born configurations, with the cuts given below) 
is shown in figure 38 for the most relevant partonic subprocesses at the Tevatron run II c.o.m. energy 
= 1.96TeV. We select events with two outgoing jets with E± > 50GeV and < 1.0, use the 
CTEQ6M parton density set [98], corresponding to as{Mz) = 0.118, and set both the renormalisation 
and factorisation scale at the Born partonic c.o.m. energy. The curves in figure 38 show a degree of 
separation between the various partonic channels — this information could perhaps be exploited in fits 
of parton distributions. 

To conclude we remark that since the only input to CAESAR is a computer subroutine for an 
observable, it offers for the first time the possibility for non-experts to easily obtain rigorous NLL re- 
summed predictions. Work remains to be done both to release the first version of CEASAR and to 
implement automated matching of NLL resummation with fixed order results. In particular, addressing 
this last issue will open the way for a vast amount of phenomenological analyses. 



17. COMBINED EFFECT OF QCD RESUMMATION AND QED RADIATIVE CORRECTION 
TO W BOSON MASS MEASUREMENT AT THE LHC ^3 

17.1 Introduction 

As a fundamental parameter of the Standard Model (SM), the mass of the H^-boson (Mjy) is of particular 
importance. Aside from being an important test of the SM itself, a precision measurement of Mw, 
together with an improved measurement of top quark mass (M^), provides severe indirect bounds on 

"Contributed by: Q.-H. Cao, C.-P. Yuan 
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Fig. 38: The resummed differential distribution at NLL accuracy for the global transverse thrust. 

the mass of Higgs boson (Mh)- With a precision of 15 MeV for Mw [158] and 2 GeV for Mt at 
the LHC [159], Mh in the SM can be predicted with an uncertainty of about 30% [158]. Comparison 
of these indirect constraints on Mh with the results from direct Higgs boson searches, at the LEP2, 
the Tevatron and the CERN Large Hadron Collider (LHC), will be an important test of the SM. In 
order to have a precision measurement of M]v, the theoretical uncertainties, dominantly coming from 
the transverse momentum of the VF-boson (Pj^), the uncertainty in parton distribution function (PDF) 
and the electroweak (EW) radiative corrections to the W boson decay, must be controlled to a better 
accuracy [158, 160]. 

At the LHC, most W bosons are produced in the small transverse momentum region. When 
is much smaller than Mw, every soft-gluon emission will induce a large logarithmic contribution to 
the Pj^ distribution so that an order-by-order perturbative calculation in the theory of Quantum chro- 
modynamics (QCD) cannot accurately describe the P^ spectrum and the contribution from multiple 
soft-gluon emission, which contributes to all orders in the expansion of the strong coupling constant a^, 
needs to be summed to all orders. It has been shown that by applying a renormalization group analysis, 
the multiple soft-gluon radiation effects can be resummed to all orders to predict the P^ distribution 
which agrees with experimental data [140, 161]. RESBOS, a Monte Carlo (MC) program [140] resum- 
ming the initial-state soft-gluon radiations of the hadronically produced lepton pairs through EW vector 
boson production and decay at hadron colliders pp/pp —>■ £il2)X, has been used by the CDF and 
D0 Collaborations at the Tevatron to compare with their data in order to determine M^- However, 
RESBOS does not include any higher order EW corrections to describe the vector boson decay. The 
EW radiative correction, in particular the final-state QED correction, is crucial for precision measure- 
ment of W boson mass at the Tevatron, because photon emission from the final-state charged lepton can 
significantly modify the lepton momentum which is used in the determination of Mw- In the CDF Run 
lb W mass measurement, the mass shifts due to radiative effects were estimated to be —65 ± 20 MeV 
and —168 ± 10 MeV for the electron and muon channels, respectively [162]. The full next-to-leading 
order (NLO) 0(a) EW corrections have been calculated [163, 164] and resulted in WGRAD [164], a MC 
program for calculating 0(a) EW radiative corrections to the process pp — > U£l{'^). However, WGRAD 
does not include the dominant correction originated from the initial-state multiple soft-gluon emission. 
The inclusion of both the initial-state QCD and final-state QED corrections into a parton level MC pro- 
gram is urgently required in order to reduce the theoretical uncertainties in interpreting the experimental 
data at the Tevatron. It was shown in Refs. [163, 164] that at the NLO, the EW radiative correction in 
pp — > lvi{'^) is dominated by the final-state QED (FQED) correction. Hence, in this paper we present 
a consistent calculation which includes both the initial-state multiple soft-gluon QCD resummation and 
the final-state NLO QED corrections, and develop an upgraded version of the RESBOS program, called 
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Fig. 39: Transverse mass distribution of boson 

RESBOS-A to simulate the signal events. Here, we only present the phenomenological impacts on 
a few experimental observables, the transverse mass of W boson (M^f ) and the transverse momentum 
of charged lepton (p^), that are most sensitive to the measurement of Mw- We focus our attention on 
the electron only, though our analysis procedure also applies to the muon. The detailed formula, the SM 
inputting parameters and the kinematics cuts are given in Ref. [165]. 

17.2 Precision Measurement of W Mass 

In Fig. 39, we show various theory predictions on the Mjf distribution. The legend of the figure is 
defined as follows: 

O LO : including only the Bom level initial-state contribution, 

O RES : including the initial-state multiple soft-gluon corrections via QCD resummation, 
O LO QED : including only the Born level final-state contribution, 
O NLO QED : including the final-state NLO QED corrections. 

For example, the soUd curve (labelled as RES-i-NLO QED) in Fig. 39(a) is the prediction from our 
combined calculation, given by Eqs. (1) and (2) of Ref [165]. 

As shown in Fig. 39(a), compared to the lowest order cross section (dotted curve), the initial state 
QCD resummation effects (dashed curve) increase the cross section at the peak of the distribution 
by about 6%, and the final state NLO QED corrections (dot-dashed curve) decrease it by about —12%, 
while the combined contributions (solid curve) of the QCD resummation and FQED corrections reduce 
it by 6%. In addition to the change in magnitude, the Une-shape of the distribution is significantly 
modified by the effects of QCD resummation and FQED corrections. To illustrate this point, we plot the 
ratio of the (RES-f-NLO QED) differential cross sections to the LO ones as the solid curve in Fig. 39(b). 
The dashed curve is for the ratio of (LO-i-NLO QED) to LO. As shown in the figure, the QCD resum- 
mation effect dominates the shape of Mjf distribution for 65 GeV < Mw < 95 GeV, while the FQED 
correction reaches its maximal effect around the Jacobian peak (Mjf ~ Mw). Hence, both corrections 
must be included to accurately predict the distribution of Mjf around the Jacobian region to determine 
M\Y- We note that after including the effect due to the finite resolution of the detector (for identifying an 
isolated electron or muon), the size of the FQED correction is largely reduced [163, 164]. 

^■^A Fortran code that implements the theoretical calculation is presented in this work. 
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Fig. 40: Transverse momentum distributions of e 



Although the Mjf distribution is an optimal observable for determining Mw at the LHC with 
a low luminosity, it requires an accurate measurement of the missing transverse momentum direction 
which becomes more difficult to control with a high luminosity option (when multiple scattering be- 
comes important). On the other hand, the transverse momentum of the decay charged lepton (p^) is less 
sensitive to the detector resolution, so that it can be used to measure Mw and provide important cross- 
check on the result derived from the Mjf distribution, for they have different systematic uncertainties. 
Another important feature of this observable is that distribution is more sensitive to the transverse 
momentum of W boson. Hence, the QCD soft-gluon resummation effects, the major source of pjf, 
must be included to reduce the theoretical uncertainty of this method. In Fig. 40(a), we show the pj^ 
distributions predicted by various theory calculations, and in Fig. 40(b), the ratios of the higher order to 
lowest order cross sections as a function of p^. The lowest order distribution (dotted curve) shows a clear 
and sharp Jacobian peak atp^ ~ Mvi//2, and the distribution with the NLO final-state QED correction 
(dot-dashed curve) also exhibits the similar Jacobian peak with the peak magnitude reduced by about 
15%. But the clear and sharp Jacobian peak of the lowest order and NLO FQED distributions (in which 
Pj^ = 0) are strongly smeared by the finite transverse momentum of the W boson induced by multi- 
ple soft-gluon radiation, as clearly demonstrated by the QCD resummation distribution (dashed curve) 
and the combined contributions of the QCD resummation and FQED corrections (solid curve). Similar 
to the distribution, the QCD resummation effect dominates the whole p^ range, while the FQED 
correction reaches it maximum around the Jacobian peak (half of M\y). The combined contribution of 
the QCD resummation and FQED corrections reaches the order of 45% near the Jacobian peak. Hence, 
these lead us to conclude that the QCD resummation effects are crucial in the measurement of Mw from 
fitting the Jacobian kinematical edge of the pj, distribution. 

As shown in Ref. [140], the effect from the initial state QCD gluon resummation to the W trans- 
verse mass distribution is dominated by the perturbative Sudkov contribution and is not very sensitive to 
the non-perturbative parameters {gi,g2 and 53) of the CSS resummation formalism. On the other hand, 
as shown in Fig. 39, the final state QED correction can largely modify the W transverse mass distribution, 
though a definite conclusion can only be drawn after including the effect of detector resolution. 

In our calculation we have included the contributions from the final state QED correction together 
with the initial state QCD resummation to predict the production and decay of W bosons produced at the 
LHC. Since the exact matrix elements have been used in the calculation, the spin correlations among the 
initial state partons and final state leptons are correctly implemented. Hence, the kinematic distributions 
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of the final state leptons, and the corresponding experimental observables, can be reliably predicted. 

In order to study the impact of the presented calculation to the determination of the W boson 
mass, the effect due to the finite resolution of the detector should be included, which will be presented 
elsewhere. 

We thank P. Nadolsky and J.W. Qiu for helpful discussions. This work was supported in part by 
NSF under grand No. PHY-0244919 and PHY-0100677. 

18. RESUMMATION FOR THE TEVATRON AND LHC ELECTROWEAK BOSON PRODUC- 
TION AT SMALL X ^5 

In the production of electroweak bosons, precise knowledge of the transverse mass Mt and transverse 
momentum qx provides detailed information about the production process, including the mass of the 
boson and associated radiative corrections. At the Tevatron, qx distributions of bosons offer insight 
into soft gluon radiation, and this information is then used for precision extraction of the W boson 
mass. At the LHC, good knowledge of the transverse distribution of Higgs bosons will be needed to 
efficiently separate Higgs boson candidates from the large QCD background. Accurate predictions for the 
small-gr region are obtained via resummation of large logarithms In^ [qT/Q] arising from unsuppressed 
soft and colUnear radiation in higher orders of perturbation theory. 

As we move from the 2 TeV Tevatron to the 14 TeV LHC, typical values of partonic momentum 

fractions x for producing W, and bosons become smaller, thus enhancing ln(l/x) terms in 
higher orders of a^. It is not entirely known how these terms (not included in a fixed-order cross section 
or conventional qx resummation) will affect W, Z^, and production at the LHC energies, in part 
because no Drell-Yan qx data is available yet in the relevant region of a; of a few 10~^ or less. 

Studies [166, 167] in the crossed channel of semi-inclusive deep-inelastic scattering (SIDIS) sug- 
gest that hadronic qx distributions at small x cannot be straightforwardly described within the Collins- 
Soper-Sterman (CSS) resummation framework [131], if the nonperturbative Sudakov function behaves 
like its large-x counterpart from the Drell-Yan process. A qx distribution in SIDIS at x < 10^^ is 
substantially broader than the conventional CSS prediction. The broadening effect can be modeled by 
including an extra x-dependent term in the Sudakov exponent. To describe the data, the extra term must 
grow quickly as x ^ 0. It noticeably contributes to the resummed form factor at intermediate impact 
parameters (6 ~ l/gy < 1 GeV~^), which hints at its origin from perturbative physics. A possible inter- 
pretation of this term is that it mimics higher-order contributions of the form ln"(l/x), which are not 
included in the resummed cross section. Due to the two-scale nature of the qx resummation problem, the 
non-resummed ln{l/x) terms may affect the qx distribution even when they leave no discernible trace 
in inclusive DIS structure functions. The DIS structure functions depend on one hard scale (of order Q), 
while the CSS resummation formula (cf. Eq. (30)) also includes contributions from large impact param- 
eters b (small momentum scales). As b becomes large, the series Q!^(1/6) ln"(l/x) in the CSS formula 
may begin to diverge at a larger value of x than the series a™((5) ln"(l/x) in the inclusive structure 
functions. For this reason, transition to fer-unordered (BFKL-like [168, 169]) physics may happen at 
larger x in qx distributions than in inclusive (one-scale) observables. 

The qx broadening discussed above was observed in semi-inclusive DIS processes. In this study, 

we explore its possible implications for the (crossed) Drell-Yan process. We begin by examining the 
resummed transverse momentum distribution for the Drell-Yan process [131], following notations from 
Ref. [170]: 

d^"? /^e-^^^'^^^(&,Q,a;A,xs)-Fy(gT,Q,XA,XB). (30) 

Here xa,b = Qe"^^ / VS, the integral is the Fourier transform of a resummed form factor W given in 
^'Contributed by: S. Berge, P. Nadolsky, F. Olness, C.-P. Yuan 
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impact parameter (6) space, and Yj^s a regular (finite at qt 0) part of tlie fixed-order cross section. In 
the small-6 limit, the form factor W is given by a product of a perturbative Sudakov exponent e""^^ and 
generaUzed parton distributions V{x,h): 



W{b,Q,XA,XB) , , =e-^^(^''3)p(x^,6)P(xB,6). (31) 



At moderately small x, where the representation (31) for W holds, we write these generalized parton 
distributions in the form 

^(^'^)L2«A-, ^ {C®f){xM/b) e-''^^)^', (32) 

where C{x,bo/b) are coefficient functions, f{x,fi) are conventional parton distributions, and bo = 
2g-7B = 1.12... is a commonly appearing constant factor. 

The term e~f'^^'>^^ in V{x,b) will provide an additional qr broadening, with an x dependence 
specified by p{x). For example, it may approximate x-dependent higher-order contributions that are not 
included in the finite-order expression for (C® /). We parametrize p{x) in the following functional form: 

(33) 

such that p{x) ~ cq/x for x <C xq, and p{x) ~ for x » xq. This parameterization ensures that the 
formaUsm reduces to the usual CSS form for large x (x » xq) and introduces an additional source of 
broadening (growing as 1/x) at small x (x <C xq). The parameter cq determines the magnitude of the 
broadening for a given x, while xq specifies the value of x below which the broadening effects become 
important. In principle, cq and xq may depend on the hard scale Q; in this first study, we neglect this 
dependence. Based on the observed dependence p(x) ~ 0.013/x at x < 10~^ in SIDIS energy flow 
data [167], we choose cq = 0.013 and xq = 0.005 as a representative choice for our plots. 

As X — >^ 0, the additional broadening term in Eq. (32) affects the form factor W both at perturba- 
tive (6 < 1 GeV~^) and nonperturbative (6 > 1 GeV~^) impact parameters. In addition, the resummed 
cross section contains conventional non-perturbative contributions from power corrections, which be- 
come important at large impact parameters {b>l GeV^^). We introduce these corrections by replacing 




the impact parameter in functions 5p and with a variable 6* = 6/y 1 + 6^/(0.25 GeV ^) [131] 

and including a nonperturbative Sudakov exponent exp {— 5Arp(6, Q)}. The function SNp{b,Q) is 
parametrized by a 3-parameter Gaussian form from a recent global fit to low-energy Drell-Yan and Teva- 
tron Run-1 data [170]. Combining all the terms, we have: 

^ = 9 J-0j2^'''^'hc<^f){xA,bo/b.){C0f){xB,bo/b.) 



Figs. 41 and 42 show the comparison of the resummed cross section (34) with the additional 
broadening term (p(x) ^ 0) to the resummed cross section without such a term (p(x) = 0). We con- 
sider cross sections for the production of and bosons, calculated according to the procedures in 
Refs. [140] and [52], respectively. The numerical calculation was realized using the programs Legacy 
and ResBos [140, 170], and with the CTEQ6M1 parton distribution functions [171]. The perturbative Su- 
dakov factor was included up to 0{a^), and the functions (C (g) /) up to 0{as). The relevant perturbative 
coefficients can be found in Refs. [131-133, 136, 172-175]. 

Fig. 41(a) shows the differential distribution da/dqx for Z boson production in the Tevatron Run- 
2, integrated over the rapidity y of the Z bosons. We observe that the distribution (34) with the additional 
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Fig. 41: qr distributions of Z'^ bosons in the Tevatron Run-2; (a) integrated over the full range of Z boson rapidities; (b) in- 
tegrated over the forward regions \y\ > 2. The solid curve is a standard CSS cross section, calculated using the 3-parameter 
Gaussian parametrization [170] of the nonperturbative Sudakov factor. The dashed curve includes additional terms responsible 
for the qr broadening in the small-a; region (cf. Eq. (34)). 



small-x term (the dashed curve) essentially coincides with the standard CSS distribution (the solid curve). 
When y is integrated over the full range, both resummed cross sections are dominated by contributions 
from X ^ 0.05 » xq, where the additional broadening (given by the function p{x)) is negligible. For 
this reason, the Tevatron distributions that are inclusive in y (e.g., the Run-1 boson data) will not 
be able to distinguish the small- j; broadening effects from uncertainties in the nonperturbative Sudakov 
function Snp- 

In contrast, the small-x broadening does lead to observable differences in the qx distributions in 
the forward rapidity region, where one of the initial-state partons carries a smaller momentum fraction 
than in the central region. Fig. 41(b) shows the cross section da/dgx for Z bosons satisfying \y\ > 2. 
The peak of the curve with p{x) 7^ is lower and shifted toward higher qx- While this difference was 
not large enough to be observed in the Tevatron Run- 1 , it seems to be measurable in the Run-2 given the 
improved acceptance and higher luminosity of the upgraded Tevatron collider. The small-x broadening 
is more pronounced in W boson production due to the smaller mass of the W boson. 

We now turn to the LHC, where the small-x broadening may be observed in the whole rapidity 
range due to the increased center-of-mass energy. Fig. 42(a) displays the distribution da/dqT for Z^ 
production with and without the small-x terms. Here, the difference is striking even if y is integrated out. 
Effects of a similar magnitude are present in W boson production, and they are further enhanced in the 
forward regions. 

The small-x broadening is less spectacular, but visible, in the production of light Standard Model 
Higgs bosons via the effective ggH vertex in the limit of a heavy top quark mass. Fig. 42(b) displays the 
resummed cross sections for production of Higgs bosons with a mass Mh = 120 GeV for several choices 
of Snp and the broadening term. We first compare cross sections for p{x) = and p{x) / (thick 
lines), where the functions p{x) and SNp{b, Q) are taken to be the same as in Z^ boson production. The 
difference between the two cross sections is not large, due to a harder qx spectrum in the Higgs boson 
case. The peaking in the g<^-dominated distribution occurs at = 10 — 20 GeV, i.e., beyond the 
region where the function p{x)lP' play its dominant role. This is different from the gg-dominated Z^ 
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distribution, where the peak is located at (?t ~ 5 — 10 GeV and is strongly affected by ex.p {— p{x)b'^} . 
Hence, for the same function p{x) as in the boson case, the difference between the curves with and 
without p{x) is minimal. 

The harder qx spectrum in the Higgs boson case is induced by a larger leading-logarithm coef- 
ficient (Ca) in gg channels, as compared to the leading-logarithm coefficient Cp in qq channels. This 
suggests that the Q-dependent part (and possibly other terms) of the nonperturbative Sudakov function 
^Arp in Higgs boson production is also multiplied by a larger color factor than in the Drell-Yan process. 
We estimate this effect by multiplying Sj\[ p by the ratio of the leading color factors in Higgs and bo- 
son production processes, Ca/Cf = 9/4 (the thin solid line). The resulting change turns out to be small 
because of the reduced sensitivity of the Higgs boson cross section to nonperturbative contributions. 

The ln(l/2;) terms may be enhanced in the case of the Higgs bosons as well, due to the direct 
coupling of the Higgs bosons to gluon ladders. At present, we do not have a reliable estimate of the 
small-x broadening in gluon-dominated channels. However, this broadening would have to be quite 
large to affect qp of 10-20 GeV or more, i.e., in the region where selection cuts on qp of the Higgs boson 
candidates will be imposed. For example, increasing the function p{x) by a factor of two as compared 
to the Z^ boson case would lead to a distribution shown by the thin dashed Une. While at qp > 20 GeV 
this effect is relatively small as compared to other theory uncertainties (e.g., higher-order corrections), it 
may affect precision calculations of qp distributions needed to separate the Higgs boson signal from the 
background in the 77 mode. 

Additional constraints on the small-x behavior of the resummed cross sections in the gg channel 
could be obtained from examination of photon pair production away from the Higgs signal region. As 
the mass of the photon pair decreases, 77 production in the gluon fusion channel via a quark box diagram 
becomes increasingly important. For instance, the subprocess gg 77 contributes up to 40% of the 
total cross section at Q = 80 GeV [176]. By comparing qp distributions in pp 77 and — > Z in the 
same region of Q, one may be able to separate the qq and gg components of the resummed cross section 
and learn about the x dependence in the gg channel. 

To summarize, we argue that a measurement of transverse momentum distributions of forward Z 
bosons at the Tevatron will provide important clues about the physics of QCD factorization and possibly 
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discover broadening of qt distributions associated with the transition to small-x hadronic dynamics. 
Based upon the analysis of qt broadening effects observed in semi-inclusive DIS, we have estimated 
similar effects in the (crossed) processes of electroweak boson production at hadron-hadron colliders. 
While the estimated impact on the Higgs boson cross section da / dqx at high qt was found to be minimal, 
much larger effects may occur in W and Z boson production in the forward region at the Tevatron Run-2, 
and at the LHC throughout the full rapidity range. If present, the small-x broadening will have to be taken 
into consideration in precision studies of electroweak boson production. Additionally, its observation will 
provide insights about the transition to fcr-unordered (BFKL-Uke) dynamics in multi-scale distributions 
at hadron-hadron colliders. 
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19. THE HIGH ENERGY LIMIT OF QCD AND THE BFKL EQUATION ^6 

The Balitsky-Fadin-Kuraev-Lipatov (BFKL) [168, 169, 177-180] formalism resums a class of large 
logarithms dominant in the Regge limit of scattering amplitudes, where the center of mass energy is 
large and the momentum transfer fixed. The cross-section for a general process A + B ^ A' + B' 
within this approach in the high energy limit can be written in the factorised form 

= / SI / S| *^(''«' ' ("»• ^ = • 

where are the impact factors characteristic of the particular scattering process and / (ka, k^, A) is 
the gluon Green's function describing the interaction between two Reggeised gluons exchanged in the 
t-channel with transverse momenta k^ b. When these two transverse momenta are large and of similar 
magnitude the Regge scale sq = |k(i| |k;,| is chosen as the scaling factor in the large logarithms. The 
implicit SQ-dependence of the NLL impact factors cancels that of the BFKL gluon Green's function 
so as to render the cross section independent of sq to this accuracy. In the leading logarithmic (LL) 
approximation terms of the form (a^ A)" are resummed, while in the next-to-leading logarithmic (NLL) 
approximation [181, 182] contributions of the type Ug (a^A)" are also taken into account. 

The gluon Green's function is obtained as the solution of an integral equation, the so called BFKL 
equation, where radiative corrections enter through its kernel. At LL it is possible to construct the 
complete eigenfunctions of this kernel and, consequently, calculate the solution analytically. At NLL 
this is only possible up to terms directly related to the running of the strong coupling; this means that 
solving the BFKL equation at NLL with the full kernel is a very challenging problem. Good progress has 
been made in the last few years in this field: Studies of the stability of the perturbative expansion were 
performed in [183-193] and of running coupling effects in [194—203]. Among the most recent work 
studying the gluon Green's function an analysis based on a new renormalisation group improved small 
X resummation scheme was proposed in Ref. [204]; in Ref. [205] an anomalous dimension including 
running couplings effects was constructed; and in Ref. [206,207] the Green's function and gluon splitting 
function were studied including a particular resummation scheme. 

In this contribution we report on the progress made in the last year to solve the NLL BFKL equa- 
tion exactly in a novel way: Using the numerical implementation of an iterative solution presented in 
Ref. [208, 209] found by explicitly separating the virtual contributions to the kernel from the real emis- 
sions in transverse momentum space. This is achieved by introducing a phase space slicing parameter in 

^^Contributed by: J. R. Andersen, V. Del Duca,A. De Roeck, A. Sabio Vera 
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dimensional regularisation. How to obtain this solution, which includes the angular correlations present 
in the kernel and the running coupling effects, is presented in the next section. 



19.1 A novel solution to the NLL BFKL equation -^^ 

It is convenient to introduce the MelUn transform in A space 

f-a+ioo 



/(k„,k,,A) 



2 ra-tt 
27r? Ja-io 



in order to write the NLL BFKL equation in dimensional regularisation as 

CO fa, {K, kb) = (5(2+2^) (k„ -H)+J d'+'^k JC {K, k + k„) (k + k„, kb) , 



(36) 



(37) 



with the kernel /C (ka, k) = 2tj(^) (ka) 5(^+2^) (k^ - k) + /Cr (ka, k) depending on the gluon Regge 
trajectory, which includes the virtual contributions, and a real emission component [181]. The delta 
function in the driving term of the integral equation corresponds to the limit of two gluon exchange. 

The phase space slicing parameter. A, is introduced through the approximation 



/^(k + k„,kb) = /^(k + k„kb)(e(k2-A2)+e(A2-k2)) 

~ (k + KM) 9 (k2 - a2) + (KM) e (A^ - k^) , 



(38) 



which is a valid one for small values of the infrared parameter A. With this separation it is possible 
to show that the e poles cancel when the real emission below the cut off is combined with the virtual 
contributions; i.e. 

cjo (q, A) = lim !^2J'^ (q) + J (f+^'klC^'^ (q, q + k) (A^ - k^) 
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+ ( o-- + oTr)ln^-6C(3) 



(39) 



Using the notation 



and 



a;o(q,A) = -^(|q|A)ln3! + «2 3^(3) 



e(X) = a,+ 



a. 



3 3 3Nc Nc "^M^ 



(40) 



(41) 



the NLL BFKL equation takes the simple form 

(a;-a;o(k„,A))^(k„,kb) = S^'^Hka-^) 



+ J d-'k (^-^^ (k2) 9 (k2 - A2) + ICr iK, K + k)) {ka + k, k^) , (42) 

where /C^ (q, q') can be found in Ref. [208]. 

Eq. (42) can now be solved using an iterative procedure in the to plane similar to the one in 
[210-212] for the LL approximation. The final solution is obtained after Mellin transform back to energy 



^^Contributing authors: J. R. Andersen, A. Sabio Vera 
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space, A. The expression for the gluon Green's function then reads (using the notation yo = 

/(k„, kb, A) = exp {uo {K, A) A) { S^^') (k„ - k^) (43) 

oo ra „ 



n=l i=l 
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(k2) +ICrlK + J2^l,K + Y,^l 
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dyi exp 
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Wo ka + ^ k;, A -Wo ka + y]k;,A 



1=1 



1=1 



(5(2) K^k, + k„-k, 



-.1=1 




where n corresponds to the number of emissions, or, alternatively, to the number of iterations of the 
kernel. 

This solution has been implemented in a Monte Carlo integration routine to study the behaviour 
of the gluon Green's function. In Ref. [208, 209] it was shown how, for a fixed value of A, only a finite 
number of iterations contributes to the final value of the solution. As the available energy in the scattering 
process increases more terms in the expansion in Eq. (43) are needed. Independence on the A scale is 
achieved when its value is small compared to the initial transverse momenta \^a,b- 

As an example of the potential of this approach we reproduce here some results. In Fig. 43 the 
value of the modulus of k^ is fixed and the dependence on the modulus of ka is studied. At LL there 
is complete agreement with the analytic solution, while the NLL result is always lower. This plot is 
calculated for a particular low value of the energy scale. The discontinuity present in this figure has its 
origin in the initial condition of the integral equation and its effect diminishes as the available energy 
in the scattering is larger. In Fig. 44 the rise of the gluon Green's function with energy is calculated. 
The slower rise at NLL compared to LL is a well known feature of the NLL corrections. The central 
lines for both the LL and NLL results are obtained by choosing the renormahsation scale fi = kb- The 
coloured bands correspond to a variation of the renormahsation scale from kb/2 to 2 ki,. An advantage 
of this numerical method of solution is that it is possible to study the angular dependences in the BFKL 
ladder. As an example we show Fig. 45 where the contribution to the NLL solution from different angles 
between the two dimensional vectors k^ is plotted. This analysis shows how the emissions are less 
correlated when the energy is larger. 

Future work using this approach will include the study of the gluon Green's function in the non 

forward case; the calculation of the solution to the = 4 Supersymmetric equation; an analysis of the 
effect of resumming the strong coupling to all orders; and the investigation of the gluon distribution at 
small X including all the scale invariant and running couphng NLL effects. 



19.2 BFKL Phenomenology at CoUiders 

When confronting BFKL predictions with data, several points are worth observing. First of all, present 
day colliders do not operate at "asymptotic energies" where the high energy exponent dominates the 
BFKL prediction, leading to a prediction of an exponential rise in cross section with an intercept of 
as4: In 2. The logarithms resummed are kinematically generated, and in the derivation of the standard 
analytic solution to the BFKL equation, the transverse momentum of the gluons emitted from the BFKL 
evolution has been integrated to infinity. It is therefore apparent that any limits on the phase space probed 
in an experiment can have a crucial impact on the theoretical prediction. Such limits can either be the cuts 
implemented in the measurement or the overall limit on the available energy at a collider. Taking hadronic 
dijet production as an example, the energy constraint will obviously not just Umit the possible rapidity 
separation of the leading dijets, but also the amount of possible radiation from the BFKL evolution, 
especially so when the leading dijets are close to the kinematical boundary. For a multi-particle final 
state described by two leading dijets with transverse momentum and rapidity ipa/bA-iUa/b) ^tid gluons 

Contributing autliors: J. R. Andersen, V. Del Duca 
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Fig. 43: ka dependence of the LL and NLL gluon Green's function at p = fc;, = 30 GeV for two values of A. 

described by {ki±,yi), the square of the total energy of the event is given by s = XaXjjS where ^/s is the 
energy of the hadron colUder and 

xai,) = ^e(-)^<' + ^e(-)^^ + J2 (44) 

^' Vs Vs 7^ ^ 

where the minus sign in the exponentials of the right-hand side applies to the subscript b on the left- 
hand side. Since in the analytic BFKL approach the contribution to x^i^b) from the BFKL radiation is 
inaccessible, this approach systematically underestimate the exact value of the x's, and can thus grossly 
overestimate the parton luminosities. 

The iterative approach of Refs. [211,212] to solving the LL BFKL equation solves this problem. 
The results of this approach coincides with the LL limit of the solution of the NLL BFKL equation of 
Sec. 19.1. At LL the change in the t-channel momentum at each step in the iteration corresponds to 
the momentum carried by one emission. The method therefore allows for the reconstruction of the full 
final state configurations contributing to the BFKL evolution, and it is possible to study quantities such 
as multiplicities and distribution in transverse momentum of the emitted gluons [213]. Only this recon- 
struction of the full final state allows for the observation of energy and momentum conservation. The 
effects of energy and momentum conservation have been studied in several processes [214—216]. When 
no phase space constraints are imposed, the iterative solution reproduces the known analytic solution to 
the LL BFKL equation. 

In Ref. [215] we have, in hght of the recent DO measurement, re-analysed the Mueller-Navelet 
proposal of the study of dijets at hadron colliders in search of BFKL signatures. The main result of 

this study is that the difference between DO and the Mueller-Navelet analysis in the reconstruction of 
the parton momentum fractions, the presence of an upper bound on the momentum transfer, and the 
contribution of the BFKL gluon radiation to the parton momentum fractions (at Tevatron energies) lower 
the parton flux in such a way as to approximately cancel the rise in the subprocess cross section with 
increasing dijet rapidity separation {ajj ~ cxp(AAy)) predicted from the standard BFKL approach. 
This strong pdf suppression is due to the dijet production being driven by the gluon pdf, which is very 
steeply falling in x for the region in x of interest. This means that even the slightest change in x has a 
dramatic impact on the parton flux and thus the prediction for the cross section. Also, the experimental 
cuts implemented have been shown to invalidate the Mueller-Navelet analysis and extraction of a 'BFKL 
intercept' . 
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Fig. 44: A-dependence of the NLL gluon Green's function evaluated for ka ~ 25 GeV and kt — 30 GeV. 

As shown by this analysis, the expected rise in cross section over the LO QCD results does not 
survive once energy and momentum conservation is taken into account at the LHC or the Tevatron. 
However, other BFKL signatures should still be present; among those studied the most is the azimuthal 
decorrelation of dijets [217, 218]. In Fig. 46 we have plotted the data and prediction for the angular 
correlation between the leading dijets as a function of the inter-jet rapidity separation at the Tevatron. 
It should be noted that no detailed jet definition was appUed in the BFKL prediction. It is, however, 
believed that a proper implementation of a jet finding algorithm will not change the partonic LL BFKL 
prediction significantly, since the gluons emitted from the BFKL evolution are typically well separated in 
rapidity. It is seen from Fig. 46 that the BFKL prediction respecting energy and momentum conservation 
is predicting slightly too much azimuthal decorrelation, but far less than an estimate based on a naive LL 
BFKL analysis ignoring the energy taken up by the BFKL radiation. On the other hand, a fixed NLO 
analysis shows too little decorrelation, while the prediction from HERWIG is in agreement with data. 
This shows that the decorrelation is dominated by soft gluon effects. 

19.3 Experimental Opportunities at the LHC 

Studies of low-x and BFKL dynamics at colliders typically require experiments with a large acceptance. 
Presently five experiments are planned at the LHC. Two of these, CMS [219] and ATLAS [220] are gen- 
eral purpose experiments with an acceptance in pseudorapidity rj of roughly \r]\ < 2.5 for tracking based 
measurements and \r]\ < 5 for calorimeteric based measurements. Here r] is defined as — lntan6'/2, 
with 6 the polar angle of the particle. Hadronic jets can be detected and measured up to approximately 
r/ = 4.5 while muons and electrons can be identified up to about r] = 2.5. Extensions of the detector 
range are being investigated, as discussed below. 

The TOTEM experiment [221] will use the same interaction point as the CMS experiment. 
TOTEM is an experiment to measure the elastic and total cross section and will use roman pot detectors 
to measure scattered protons in elastic and diffractive pp interactions, and charged particle detectors for 
tagging inelastic events in the regions 3 < jr^l < 5 (Tl) and 5.3 < |r/| < 6.5 (T2). CMS is studying 
to add a calorimeter, called CASTOR, in the forward region directly behind T2. CASTOR will have 
an electromagnetic and hadronic readout section and an acceptance in the range 5.4 < |?7| < 6.7. The 
TOTEM trigger and data readout will be such that these can be included in the CMS datastream, such 

^'Contributing authors: A. De Roeck 
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Fig. 45: The dependence of the gluon Green's function at NLL on the angle between (ka = 25 GeV) and k;, (ki, = 30 GeV) 
for the choice of A = 3 and A = 5. The renormaUsation point is chosen at /x = kb- 

that both experiments can run as one experiment. That way, a full detector with good coverage for jets 
and electrons/photons with rapidities of up to 7 will be available. The position of the Tl, T2 trackers and 
the forward calorimeter, along the beamline and integrated with CMS, is shown in Fig. 47. 

Roman pot detectors are presently also considered by ATLAS [222]. ATLAS further plans for a 
cerenkov/quartz fiber detector, called LUCID, for luminosity monitoring, which will cover the region 
5.2 < r) < 6.2. This detector could possibly be used for the tagging of rapidity gaps in an event. 

The extended coverage of the detectors will allow to probe processes at low-x via Drell-Yan, jets, 
W and direct photon production. E.g. the Drell-Yan process qq e~^e~ has a simple experimental 
signature. The xi,2 values of the two incoming quarks relate to the invariant mass of the two electron 
system Mge as • a;2 • s — Mgg, hence when one of the xi^2 values is large (say x > 0.1), very low-x 
can be probed with low mass Drell-Yan pairs. From Fig. 48 [223] we observe that in order to reach small 
masses (small scales) and low-x, will require to probe large values of rj. Thus the resulting electrons will 
dominantly go in the very forward direction. 

If low mass Drell-Yan pairs, prompt photons or jets (which have the same kinematics) can be 
measured in the CMS/TOTEM forward detectors, the parton distributions can be probed down to values 
of 10~^ — 10~^, i.e. at lower values than for any other collider Depending on the low-x dynamics, 
the predictions for jet, Drell-Yan and other cross sections can differ by a factor two or more in this 
region [224]. Experimental challenges on extracting the low-x signals with respect to machine and QCD 
backgrounds are still under study. 

The extended detector capabilities also allow for larger rapidity distances (i.e longer gluon lad- 
ders) between two measured hard jets, one of the golden BFKL measurements at a hadron collider (see 
Sec. 19.2). Note that the experiments will also measure all activity in between these two jets. Since 
BFKL does not have strong kr ordering the number and energy spectrum of the mini-jets produced 
along the ladder may be a significant footprint for BFKL. The main challenge will be to define a suitable 
trigger for this two jet channel. The presently foreseen trigger thresholds in ATLAS and CMS are around 
150-200 GeV for one and two jet events, while for BFKL studies particularly jets in the energy range 
of 20-100 GeV are more of interest, since the distributions of these will be less distorted by kinematical 
constraints. 

In principle it is also possible to run the LHC at lower energies, e.g. at 3 and 8 TeV, with significant 
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Fig. 46: Dijet angular correlation at the Tevatron. DO data are marked by the diamonds, BFKL prediction (LL supplemented 
by running coupling and including energy/momentum conservation) given by the stars. The LL BFKL prediction without 
energy/momentum conservation is marked by the solid line. Figure taken from Ref. [212]. 



luminosities (> lO^^cm ^s allowing to make parton density independent ratios of di-jet production 
at different energies, over a large range of rapidity. 

For luminosities up to 10^^cm~^s~^ or below, the number of overlap events in one bunch cross- 
ing is small enough such that rapidity gaps can be used to detect interactions with colour singlet ex- 
change [225]. The extended range covered in ry will allow to study events with multi-gaps. Hard diffrac- 
tive phenomena, often linked to small-x dynamics, can be studied further in detail via the tagging of the 
scattered protons in the roman pot detectors along the beam-line. Together with the measurement of 
the hard probes in the central LHC detectors, these events will allow for a detailed study of diffractive 
phenomena in pp collisions at the highest energies. 



MB/ 1/4 




Fig. 47: Position of the inelastic event tagging detectors of TOTEM, Tl, T2 and CASTOR integrated with CMS 
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LHC parton kinematics 




Fig. 48: The kinematic plane (a, Q^) and the reach of the LHC, together with that of the existing data (HERA, fixed target). 
Lines of constant pseudo-rapidity are shown to indicate the kinematics of the produced objects in the LHC centre of mass 
frame [223]. 
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20. PION PAIR PRODUCTION AT THE LHC : COMPARING QCD @ NLO WITH PYTHIA 
20.1 Introduction 

The indirect mass bounds for the Higgs boson from precision measurements indicate that the favoured 
range is below 200 GeV. A prominent search channel for neutral Higgs bosons in the mass window 
80 GeV < niH < 140 GeV is the rare decay H 77, as electromagnetic signals are preferable 
compared to hadronic ones due to huge backgrounds for the latter. The important backgrounds for this 
channel are prompt di-photon production and processes with one real photon and a jet with an isolated 
neutral meson. In addition, the abundant two jet rate at the LHC leads to a large reducible background 
where a meson in both jets produces an electromagnetic signature in the detector. This background needs 
to be rejected very efficiently. It is thus an important question to predict the respective rates as precisely 
as possible to allow for reliable simulations. 

Experimental studies dominantly use PYTHIA [10] to predict the reaction rates for signal and 
background processes. PYTHIA is based on leading order 2^2 matrix elements combined with initial 
and final state parton showers coupled to a hadronization model. One can expect that calculations in 
next- to-leading order in the QCD coupling lead to an improvement of the predictions in two respects. 
Firstly, distributions which are sensitive to hard radiation of an extra parton will not be described properly 
by a parton shower which is on the other hand expected to describe the colUnear regions in phase space 
more reliably than a fixed order calculation. Secondly, inclusion of loop effects typically reduces the 
dependence on unphysical scales and thus leads to more reliable predictions. 

The experimental studies show that the reducible pion-i-photon and di-pion contribution can be 
reduced by isolation cuts to a level which is acceptable for Higgs searches in the two photon channel 
[226, 227]. It is the aim of this study to cross check if this statement holds also at next-to-leading order 
(NLO) precision. We compare in the following production rates of pion pairs at high pT calculated 
with PYTHIA and the DIPHOX code [228]. The latter contains partonic matrix elements at NLO for 
the production of photon pairs, photon+hadron and hadron pairs and is flexible enough to account for 
diverse experimental cuts which can be modeled on the partonic level. For the description of hadrons 
in DIPHOX the model of collinear fragmentation is used. Partonic matrix elements are folded with 
fragmentation functions, Dp^^^x, Q) which stand for the probabihty to find a hadron h with an energy 
fraction z relative to the jet formed by the initial parton p = g,q. This probability is scale dependent. The 
scale dependence is governed by AltarelU-Parisi evolution. Initial distributions are fitted to experimental 
data in a certain x range and extrapolated to the high (low) x range by some standard parameterizations. 
We note at this place that in the case of severe isolation cuts typically the tails of the fragmentation 
functions become numerically important which are not experimentally well restricted [229]. 

DIPHOX has been confronted with many different data sets from fixed target experiments to 
Tevatron data and showed an excellent agreement between di-pion or di-photon observables, whenever 
the condition for fixed order perturbative calculations were fuUfilled [228, 230, 231]. We are thus in 
a situation to make a NLO prediction for photon pair, photon-i-pion and pion pair production at the 
LHC [229, 232, 233]. We want to mention that there is another calculation for di-hadron production at 
NLO available [234]^\ 

A similar study for photon-i-pion production was presented by us in a previous Les Houches report 
[134]. We found good agreement between PYTHIA and DIPHOX for loose isolation cuts whereas for 
severe isolation cuts the PYTHIA rates were above the NLO predictions. Qualitatively we observed that 
isolation cuts are more efficient for the NLO calculation than in the PYTHIA result. We were lead to the 
conclusion that the experimental studies are conservative for the photon pion background. 

^"Contributed by: T. Binoth, K. Lassila-Perini 

^'in [176] corrections for two photon production via a quark loop at the two-loop level are taken into account. 
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20.2 Comparison: PYTHIA VS. DIPHOX 

The basic underlying reaction for the production of pion pairs at leading order is two jet production with 
numerous different subprocesses like q + q^q' + q',q + q'^q + q',q + q^g + g,q + g^p + g, 
g + g ^ g + g- It is a priori not clear, if there is a dominant process for different kinematical regions 
and experimental cuts. All subprocesses have to be considered. At NLO virtual corrections to the 2^2 
processes have to be combined with numerous 2^3 matrix elements with one unresolved particle. 
Details of the calculation can be found in [228]. 

We compared two distributions relevant for Higgs boson searches, first the invariant mass dis- 
tribution of the pion pair, da/dM^^T^ with Mt^tt = {p^j + Pwi)'^ and second the transverse momentum 
distribution of the pion pair, da/dqx with qx = \pt -ki + Pt ^21- 

PYTHIA [10] version 6.205 was used to generate the events. High px QCD processes were 
simulated and the events with two neutral pions were accepted. The pions were required to have pt > 
25 GeV and to be in the pseudo-rapidity range of \r]\ < 2.5. The isolation was implemented by setting 
a threshold to the sum of Et of all particles (excluding the pion itself) in the isolation cone AR = 
^(Ar/)2 + A0)2. 

The default parameters were chosen, apart from the longitudinal fragmentation function, where a 
hybrid scheme was used and for the multiple parton interactions, where a set of parameters tuned with 
the collider data [235] were used. The pile-up effects from collisions within the same bunch crossing 
were not taken into account, as we are here basically interested in a comparison of event rates of partonic 
interactions. We note that the multiple interactions belong to the regime of soft QCD, which is modeled 
in PYTHIA. The relevance of multiple interactions is restricted to pr values which are much below the 
scales used in our study. We are not in a regime where these effects are expected to become numerically 
sizable. 

As already mentioned DIPHOX contains all matrix elements indicated above together with their 
next-to-leading order virtual and real emission corrections. Using the same experimental cuts defined 
above, the same distributions as with PYTHIA were produced by varying the renormalization (//), fac- 
torization (M) and fragmentation (Mf) scales. We took /i = M = Mj = cM^^/2 with c = 1/2, 1, 2. 
For the parton distribution functions we used [236], for the pion fragmentation function we used a recent 
parametrization of Kramer, Kniehl and Potter (KKP) [237] and an older one from Binnewies, Kramer 
and Kniehl (BKK) [238]. This was done to get an idea about the uncertainties due to fragmentation 
functions, although the KKP set should be more reliable, as it contains newer data. 

In Fig. 49 we show the invariant mass distribution for two different isolation criteria, Exniax = 
20, 100 GeV in a cone AR < 0.4. LO and NLO matrix elements were used for the DIPHOX predictions 
for three different choices of the scales. Higher scales mean lower cross sections in the plot. In the left 
plot KKP fragmentation functions were used, on the right the BKK parametrization. The PYTHIA pre- 
diction is shown as a full line. One observes that for ETmax = 100 PYTHIA undershoots the DIPHOX 
curves considerably. For the KKP fragmentation functions there is difference between a factor 2 to 6 for 
LO and 3 to 6 for NLO matrix elements. By using the BKK fragmentation functions, LO matrix ele- 
ments and high scales (c = 2) the accordance is improved. For the harder isolation cut, Exmax = 20, the 
PYTHIA prediction shows a fair agreement with the DIPHOX result for both fragmentation functions. 
It is inside the scale variation of the NLO result apart for the last bin which is slightly below in the KKP 
case. One observes that the isolation cut acts more efficiently in DIPHOX than in PYTHIA. This was 
observed already in a earlier study [134]. It has to do with the different fragmentation models used. For 
more severe isolation preliminary results indicate that the PYTHIA distributions start to fall below the 
DIPHOX curves 

The inclusion of NLO matrix elements increases the cross section. In Fig.50 the ratio of NLO and 
LO result for KKP (left figure of Fig.49) is shown. Evidently the K-factors defined as this ratio depend 

^^The curve for ETmax = 10 GeV is not included because more statistics is needed to make a definite statement. 
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Fig. 49: The invariant mass distribution of pion pairs at LHC for loose (ETmax = 100 GeV, AR < 0.4) and a more realistic 
isolation {ETmax = 20 GeV, AR < 0.4) criterion. The full line is the PYTHIA prediction. For DIPHOX NLO (full triangle) 
and LO (open triangle) matrix elements were used for 3 different scale choices ii = M = Mf = c M,r7r/2, c = 0.5 (upper), 
c = 1 (middle), c = 2 (lower). The left plot uses the KKP, the right plot the BKK fragmentation functions for the neutral pions. 

on the scale choice and the experimental cuts, they vary between 1 and 2 for ETmax = 100 GeV and 2 
and 3 for Exmax = 20 GeV, see Fig. 50. The inclusion of higher order terms improves the stability of 
the DIPHOX prediction under scale variations. This is much less pronounced for the harder isolation cut 
with ETmax = 20 GcV Here compensations between LO and NLO terms are spoiled by the isolation cut. 
We note that asymmetric pT cuts, not applied here, lead to larger K-factors for the same reason [233]. 

Finally we present the transverse momentum distribution of the pion pair in Fig. 51. Now, only the 
NLO prediction is shown as this distribution is sensitive to 2 — >^ 3 matrix elements not present in the LO 
calculation. For a partonic 2^2 reaction there is a kinematical restriction pT < ETmax- Distributions 
beyond this point are filled by the parton shower in the PYTHIA case and with hard 2^3 matrix 
elements in DIPHOX which are not present in PYTHIA. For ETmax = 100 GeV this does not affect the 
plotted range. Again, whereas PYTHIA undershoots the NLO curves considerably for ETmax = 100 
GeV, the shape is very similar. Fair agreement is found for ETmax = 20 GeV for not too high values of 
qt- The drop in statistics of the PYTHIA curve can be understood by the missing 2—^3 matrix elements 
which are important for high qt values in the case of hard isolation. 

20.3 Conclusion 

We have presented a comparison between PYTHIA version 6.205 and DIPHOX version 1.2 for the pair 
production of high pT pions at the LHC. This is an important reducible background for Higgs search 
in the two photon channel. We observe that for loose isolation cuts the NLO result of DIPHOX is 
significantly higher than from PYTHIA which is due to matrix elements not present in PYTHIA and 
different fragmentation models in both codes. As isolation cuts act harder in DIPHOX than in PYTHIA, 
due to different fragmentation models, there is a fair agreement of both codes for a harder isolation cut. 
The trend is that for more severe isolation cuts PYTHIA will even overshoot the NLO predictions. We 
note that a further suppression of the dipion rate than shown in our plots is possible for example by using 
asymmetric pt cuts and pion/photon identification methods, but this is beyond the scope of this study. 
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Fig. 50: K factors (top) and relative scale variations (bottom) for the two isolation criteria. The scale variations are plotted 
for the LO (dashed) and the NLO prediction (full). The scales have been varied as explained in the text. KKP fragmentation 
functions have been used. 




Fig. 51: The qr distribution of the pion pair. PYTfflA vs. DIPHOX NLO. KKP (full triangle) and BKK (open triangle) 
fragmentation functions are used. 
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We arrive at the tentative conclusion that experimental studies based on PYTHIA seem to be 
slightly below the more complete NLO calculation although the discrepancie is not dramatic. We have to 
stress that this statement depends strongly on the behaviour of the fragmentation model. Hard isolation 
means that fragmentation models are tested at high x (the ratio between the pion energy and the jet 
energy). Currently the fragmentation functions are not well constrained experimentally in the high x 
tails. Hopefully the situation improves soon by analyzing high pT pions at the Tevatron. Our study 
strongly motivates a refitting of the fragmentation functions which eventually leads in turn to a good 
understanding of the two photon rate at the LHC which is, apart from Higgs search, an interesting 
observable concerning our understanding of QCD. 



Acknowledgements 

We would like to thank Sasha Nikitenko for many useful discussions and also the conference organizers 
for the pleasant workshop. 



21. QCD-INDUCED SPIN PHENOMENA IN TOP QUARK PAIR PRODUCTION AT THE 

LHC 33 

21.1 Introduction 

The NLO QCD corrections to hadronic tt production have been known for quite some time [239-242]. 
These results were refined by resummation of soft gluon [243] and threshold [102,244] logarithms. These 
predictions were made for the production of tt pairs averaged over their spins. Because of the extremely 
short top-quark Ufetime the spin properties of top quarks are transferred to its decay products without be- 
ing diluted by hadronization. Thus quantities that involve the t and/or t spin are also "good" observables 
in the sense that they can be rehably calculated perturbatively, in particular within perturbative QCD. It 
is expected that such observables will be very useful in exploring the interactions that are involved in top 
quark production and decay. Besides standard model (SM) studies they allow, for instance, the search 
for non-SM interactions, in particular CP violating interactions of top quarks [245-247], or to pin down 
the nature of heavy resonances that strongly couple to top quarks [247,248], if such objects exist. 

Needless to say, searches for new interactions with top spin observables will require, on the the- 
oretical side, rather precise SM precictions. While SM interactions induce only small polarizations of t 
and t quarks in hadronic tt production, the QCD-induced correlation between the t and t spins is large. 
They can be studied both at the Tevatron and at the LHC by means of double differential angular distri- 
butions of ti decay products. In this contribution we present results [249-252] for these distributions at 
NLO QCD. 



21.2 Theoretical framework 

We consider hadronic ti production and subsequent decay into the following channels: 

i+l'- + X 



hi h2 



tt + X 



(45) 



Jt Jt 



+ X 



Here /ii 2 = p,p, = e,^,T, and jt (jt) denotes a jet originating from nonleptonic t (t) decay. A 
complete next-to-leading order QCD analysis of (45) involves the parton reactions 

tt 



gg, qq — 'bb + 4/, 



(46) 



Contributed by: W. Bernreuther, A. Brandenburg, Z.G. Si, P. Uwer 
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gg,qq^bb + Af + g, (47) 

9 + q{q)^bb + 4f + q{q). (48) 

where f = q,i, V£. In view of the fact that the total width of the t quark is much smaller than its mass, 
Vt/mt ~ 1%, it is adequate to use the leading pole approximation (LPA) in calculating the NLO dif- 
ferential cross sections of the above parton reactions. Within the LPA, the radiative corrections can be 
classified into factorizable and non-factorizable contributions. The factorizable corrections were com- 
puted in [249-252] in the narrow width approximation Tt/nit for which the contributions to the 
squared matrix element M of the respective parton reaction are schematically of the form 

\M\^ OC Tr [pRp\ = pa'aRaa',/3p'Pf3'/3- (49) 

Here R denotes the density matrix which describes the production of on-shell tt pairs in a specific 
spin configuration, and p, p are the density matrices describing the decay of polarised t and t quarks, 
respectively, into specific final states. The subscripts in (49) denote the t, t spin indices. The spin- 
averaged production density matrices yield the NLO cross sections for M being produced by qq, gg, gq, 
and gq fusion, which were first determined by [239-242]. At NLO the qq and gg fusion reactions receive 
also non-factorizable corrections. They were computed by [253] in the semi-soft gluon approximation. 
While these contributions have an impact on, e.g., the t, t, and tt invariant mass distributions, they cancel 
in the integrated NLO cross sections of the above reactions [254,255]. Moreover, it can be shown [252] 
that these non-factorizable corrections do not contribute either to the double differential distributions 
which will be discussed below. Therefore we will not discuss them here any further. 

The factorizable contributions (49) must be consistently evaluated to order a^. This involves 
also the matrix elements to order ag of the main SM decay modes of the (anti)top quark in a given 
spin state, that is, the semileptonic modes t bl^Vf, bi^veg {i = e, p, r), and the non-leptonic decays 
t bqq' , bqq'g where qq' = ud, cs for the dominant channels. For the computation of the double angular 

distributions (53) the matrix elements of the 2-particle inclusive parton reactions i ai -I- 02 -I- X are 
required. Here ai, 02 denotes a lepton or a jet. In the LPA this involves the 1 -particle inclusive t decay 
density matrix = r(^)(]l + ki r • qi)Q'Q, where qi is the direction of flight in the t rest frame 

and r^^) is the partial width of the respective decay channel. An analogous formula holds for i decay. 
The factor ki is the t spin analyzing power of particle/jet ai. It is clear that its value is crucial for the 
experimental determination of top spin effects, in particular of tl spin correlations. For the standard 
V — A charged current interactions these coefficients are known to order for semileptonic [256] and 
non-leptonic [257] modes. The charged lepton is a perfect analyzer of the top quark spin, which is due 
to the fact that 

K£ = 1 - O.OlSa^ . (50) 

In the case of hadronic top quark decays the spin analysing power of jets can be defined. A detailed 
analysis was made in [257]; we give here only two examples: 

Kb = -0.408 X (1 - 0.340as) = -0.393 , (51) 
Kj = +0.510 X (1 - 0.654as) = +0.474 , (52) 

where kj, is the analysing power of the 6 jet and kj refers to the least energetic non-b-quark jet defined 
by the Durham algorithm. Using hadronic final states to analyze the spins of t and/or t results in a loss of 
analyzing power. However, this is (over)compensated by the gain in statistics and by the efficiency with 
which the t (t) rest frames can be reconstructed. 

Let us now discuss double angular distributions da / d cos 6 id cos 62 for the channels (45) 
which are appropriate observables to measure the (QCD induced) ti spin correlations. Here ^1(^2) is 
the angle between the direction of flight of particle/jet 01(02) in the t{t) rest frame with respect to ref- 
erence directions a (b) which will be specified below. We define the t (t) rest frame by a rotation-free 
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boost from the ti zero-momentum frame, for reasons given below. As mentioned above, non-factorizable 
corrections do not contribute, at NLO, to this type of distributions [252]. Thus, to this order in the QCD 
coupling a (b) can be interpreted as t{t) spin quantization axes, see below. Integrating over the full phase 
space and choosing a, b as given below it can be shown [252] that these distributions have the generic 
form^^ 

-1 A - = -il-Ccos9icos92). (53) 

a a cos Old cos 02 4 

The coefficient C which signals the correlation between the ti spins depends, for a given reaction, on the 
choice of a and lo. For the factorizable corrections the exact (to all orders in a^) formula C = K1K2D 
holds [251]. Here D is the double spin asymmetry 

^ ^ Nm)+N[ii)-N{U)-N{i^) 
N{]])+N{ii) + N{]i) + N{i]y 

where A^(TT) denotes the number of ti pairs with t (i) spin parallel to the reference axis a (lo), etc. Thus 
a and b can be identified with the quantization axes of the t and i spins, respectively, and D directly 
reflects the strength of the correlation between the t and i spins for the chosen axes. 

For ti production at the Tevatron it is well known that the so-called off-diagonal basis [260], which 

is defined by the requirement that (5"(tJ,) = (T(J,t) = for the process qq ti at tree level, yields a 
large coefficient D. It has been shown in [251] that the beam basis, where a and b are identified with 
the hadronic beam axis, is practically as good as the off-diagonal basis. A further choice is the helicity 
basis, which is suitable for the LHC. 

An important theoretical issue of top quark spin physics beyond leading-order QCD is the con- 
struction of infrared and collinear safe observables at the parton level. In the case at hand it boils down to 
the question in which frame the reference directions a and b are to be defined. It has been shown [252] 
that, apart from the t and i rest frames, the ti zero momentum frame (ZMF) is the appropriate frame for 
defining collinear safe spin-momentum observables. The off-diagonal, beam, and helicity bases which 
we refer to in the next section are defined in the tiZMF. Details will be given in [252]. 



21.3 Spin correlations at nlo: predictions for the Tevatron and the LHC 

We have computed the 2-particle inclusive differential cross sections for (45) and, in particular, the 
double angular distributions (53) to NLO QCD, with ag and the top mass being defined in the MS and 
in the on-shell scheme, respectively. The mass factorization was performed in the MS scheme.) 

In Table 3 we list our predictions [252] for the spin correlation coefficient C in the double differ- 
ential distribution (53) at the Tevatron for the three different choices of references axes a, lo discussed 
above. We use the CTEQ6L (LO) and CTEQ6.1M (NLO) parton distribution functions (PDF) [171]. 
Numbers are given for the dilepton, lepton-|-jet and all-hadronic decay mode of the ti pair, where in the 
latter two cases the least energetic non-6-quark jet (defined by the Durham cluster algorithm) was used 
as spin analyser. One sees that the spin correlations are largest in the beam and off-diagonal basis. The 
QCD corrections reduce the LO results for the coefficients C by about 10% to 30%. 

For the LHC it turns out that the spin correlations w.r.t. the beam and off-diagonal basis are quite 
small due to a cancellation of contributions from the gg and qq initial states. Here, the helicity basis is 

a good choice, and Table 4 shows our results for the C coefficient in that case. The QCD corrections 
are smaller for the LHC than for the Tevatron; they vary between 1 and 10%. For both colhders the 
relative corrections |(Cnlo — C'lo)/C'lo| are largest for the all-hadronic decay modes and smallest for 
the dilepton decay modes. 

'''QCD (or SM) generated absorptive parts in the parton scattering amplitudes induce a small t and t polarization, which to 
order a'l is normal to the qq, gg — » scattering planes [258,259]. 



81 







dilepton 


lepton+jet 


jet+jet 


Chel 


LO 


-0.471 


-0.240 


-0.123 




NLO 


-0.352 


-0.168 


-0.080 


Cbeam 


LO 


0.928 


0.474 


0.242 




NLO 


0.111 


0.370 


0.176 


Coff 


LO 


0.937 


0.478 


0.244 




NLO 


0.782 


0.372 


0.177 



Table 3: LO and NLO results for the spin correlation coefficient C of the distribution (53) in the case of pp collisions at 
y/s = 1.96 TeV for different decay modes. The PDF CTEQ6L (LO) and CTEQ6.1M (NLO) of [171] were used, and 
H = mt = 175 GeV. 





dilepton 


lepton+jet 


jet+jet 


Chel LO 


0.319 


0.163 


0.083 


NLO 


0.326 


0.158 


0.076 



Table 4: Results for Chei for pp collisions at = 14 TeV using the same PDF and parameters as in Table 3. 

An interesting aspect of these double distributions is their high sensitivity to the quark and gluon 
content of the proton [251]; the reason being that the contributions to C from qq and gg initial states 
have opposite signs. Table 5 shows, for dilepton final states, the dependence of the NLO results on the 
choice of the PDF. While the results for the recent CTEQ6.1M and MRST2003 [261] PDF agree at the 
percent level (this is not the case for previous versions of the CTEQ and MRST PDF), the GRV98 [262] 
PDF give significantly different results at the Tevatron. This shows that measurements of (53) may offer 
the possibiUty to further constrain the quark and gluon content of the proton. 



Tevatron 




CTEQ6.1M MRST2()03 


GRV98 


Chel 


-0.352 -0.351 


-0.313 


Cbeam 


0.777 0.777 


0.732 


Coff 


0.782 0.782 


0.736 


LHC 


Chel 


0.326 0.327 


0.339 



Table 5: Spin correlation coefficients at NLO for different PDF for pp at \/s = 1.96 TeV (upper part) and pp at = 14 TeV 
(lower part) for dilepton final states. 

All the results above were obtained with jj. = /jir = np = rrit = 175 GeV. A variation of the 
scale between mt/2 and 2mt changes the central results for C at = mt by roughly ±5%. Varying 
mt from 170 to 180 GeV changes the results for C at the Tevatron by less than 5%, while for the LHC, 
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Chel changes by less than a percent. 
21.4 Conclusion 

Top-antitop spin correlations, which are predicted to be large within the SM, are expected to become a 
good tool for analyzing in detail top quark pair production and decay dynamics. They can be studied at 
the Tevatron and - in view of the expected large ti data samples - especially at the LHC in the dilepton, 
single lepton and all-hadronic decay chaimels by measuring suitably defined double angular distributions. 
While the NLO QCD corrections to these distributions are of the order of 10 to 30% for the Tevatron, 
they are below 10% for the LHC. Work on soft gluon and threshold resummations will further reduce the 
theoretical uncertainties. 

22. QCD RADIATIVE CORRECTIONS TO PROMPT DIPHOTON PRODUCTION IN ASSO- 
CIATION WITH A JET AT THE LHC 

22.1 Introduction 

Higgs production in association with a jet of high transverse energy with a subsequent decay into two 
isolated photons, pp ^ if + jet ^ 77 + jet, is considered a very promising discovery channel for a 
Higgs boson of intermediate mass (100 GeV < uih < 140 GeV) [263-265]. The main background to 
this signal is the pp — > 77 + jet chaimel, where the photons are isolated. In this channel the signal is a 
small and narrow peak on a flat background [265]. Thus the QCD prediction is not needed for predicting 
the background, which can be well measured from the sidebands, but rather to optimize the selection and 
isolation cuts for the experimental search. For this purpose usually a Monte Carlo program is used [70], 
which however, does not take into account the QCD radiative corrections. These corrections are large 
and strongly dependent on the photon isolation parameters [265] and therefore caimot be ignored in the 
analysis. 

In perturbative QCD, the cross section for the production of a single isolated photon in a colhsion 
of two hadrons A and B of momenta pA and ps, respectively, has the following general form: 

dcrAB{PA,PB;P'r) = X] y^^a fa/AiXa, l^F)h/B{xb, fJ-p) da'^l^{XaPA, Xi,Pb;P^; fiR, HF, fJ-y) 
a,b 

+ XI ^^f> fa/A{Xa, IJ'F)fb/BiXb, I^f) ^Kb^c (xaPA, XbPB', i^R^ I^F' /^t) -^7/0(2, fJ'j) ■ 

a,b,c 

(55) 

The first term is called the direct component, where the subtracted partonic cross sections da^]^ get con- 
tributions from all the diagrams with a photon leg. The second term is called fragmentation component, 
where the subtracted partonic cross sections da^^^], get contributions from diagrams with only coloured 
external partons, with one of the partons eventually fragmenting into a photon, in a way described by 
the (perturbatively uncalculable but universal) parton-to-photon fragmentation function D^/^- If there 
are two isolated photons, the equation contains four terms, a double direct, two single-fragmentation 
and a double-fragmentation component [228]. The direct and fragmentation components are not defined 
unambiguously (except for the case of the 'smooth' photon isolation to be used here), finite terms can be 
shifted among the terms, only the sum is meaningful physically. The precise definition of the direct and 
fragmentation terms, valid to all orders in perturbation theory, can be found for instance in Ref. [266]. 

In perturbation theory beyond leading order, the isolated photon cross section is not infrared safe. 
To define an infrared safe cross section, one has to allow for some hadronic energy inside the photon 
isolation cone. In a parton level calculation it means that soft partons up to a predefined maximum 

^'Contributed by: V. Del Duca, F. Maltoni, Z. Nagy, Z. Trocsanyi 
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energy are allowed inside the cone. This is also natural in the experiment: complete isolation of the 
photon is not possible due to the finite energy resolution of the detector. 

There are two known ways to implement the photon isolation. The standard way of defining an 

isolated prompt photon cross section, that matches the usual experimental definition, is to allow for 
transverse hadronic energy inside the photon isolation cone up to i?_L,max = ^P7±» with typical values 
of e between 0.1 and 0.5^^, and where p^±_ is taken either to be the photon transverse momentum on 
an event-by-event basis or to correspond to the minimum value in the range. Using this isolation 
prescription, both the direct and the fragmentation terms contribute to Eq. (55). In Ref. [266] it was 
shown that a small isolation cone for the photon leads to unphysical results in a fixed order computation. 
For a small cone radius R^, an all-order resummation of ln(l/i2^) terms combined with a careful 
study of the border Une between perturbative and non-perturbative regions has to be undertaken. 

Even with resummation, the very narrow isolation cone is not favoured from the practical point of 
view. The ratio of the fragmentation component to the direct one slowly increases with decreasing i?^. 
Since the fragmentation function is non-perturbative, therefore, it has to be measured and the available 
functions have relatively large errors. The fragmentation component is also strongly dependent on e, 
increasing rapidly with increasing e. Thus, the theoretical uncertainty related to the photon fragmentation 
can be decreased with large cone sizes and small e. 

In order to avoid completely the uncertainties related to the fragmentation component, Frixione 
introduced a 'smooth' photon-isolation when the fragmentation contribution is zero [267]. This isolation 
means that the energy of the soft parton inside the isolation cone has to converge to zero smoothly if the 
distance in the rj — (/) plane between the photon and parton vanishes. Explicitly, the amount of hadronic 
transverse energy E± (which in our NLO partonic computation is equal to the transverse momentum of 
the possible single parton in the isolation cone) in all cones of radius r < must be less than 



The smooth isolation prescription can be viewed as if the singularity in the quark-photon splitting were 
factorized into the fragmentation component using phase-space cuts such that the finite remainder in the 
fragmentation is shifted completely into the direct component in a perturbatively computable way. 

The two prescriptions have their own advantages and less appealing features. The standard isola- 
tion is easier to employ experimentally, but has the following disadvantages from the theoretical point 
of view: (i) it involves a non-perturbative component in its theoretical prediction (fragmentation) with 
large errors due to pure experimental information; (ii) the theoretical prediction of the fragmentation is 
more difficult to compute than that of the direct component. The smooth isolation is easier to implement 
in NLO perturbation theory. However, it is less favoured in the experiment because (i) it is difficult to 
employ it due to the finite granularity of the detector; (ii) it was found to be less efficient for isolat- 
ing photons from jets than the cone isolation, approaching the efficiency of the latter with n tending to 
zero [268]. 

Experimentally the standard cone isolation is recovered from the smooth isolation when n ^ 0. 
The smooth isolation leads to smaller cross sections for the same values of the cone radius and e, because 
it means a more severe cut into the phase space. However, due to the finite granularity of the detector, 
the two prescriptions yield equal cross sections even for small, but non-zero values of n. Small values of 
n can be chosen in perturbation theory, but n = caimot be taken without including the fragmentation 
component. If we are interested only in the Higgs-boson search, we need a perturbative prediction that 
is reliable for those values of the isolation parameters which are found useful to increase the signal 
significance. Studying the dependence of the perturbative prediction with the smooth prescription on the 
isolation parameters, we can estimate the reliability of the theoretical predictions. 

^*In experiments, often a fixed value of -Ex, max in the order of several GeV is used. 




(56) 
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In order to assess the dependence of the radiative corrections on the isolation parameters, we use 
a partonic Monte Carlo program which employs the dipole subtraction method [269], slightly modified 
for better numerical control [270] as implemented in the NLO JET++ package [271], and use the smooth 
isolation prescription at the parton level. We compute cross section values at leading order and at NLO 
for a Large Hadron Collider (LHC) running at 14 TeV. The values shown at leading order were obtained 
using the leading order parton distribution functions (p.d.f.'s) and those at NLO accuracy were obtained 
using the NLO p.d.f.'s of the CTEQ6 package [98] (tables cteq611 and cteq6m, respectively). We used 
the two-loop running of the strong coupling at NLO with Aqcd = 226 MeV and one-loop running with 
Aqcd = 165 MeV at leading order. The renormalization and factorization scales are set to /ir = fip = 
XfM/J'O, where for the reference value /iq we use /Xq = (m^^ + pjetx)/^ ^^^^ '^t^ the invariant mass of the 
photon pair. This definition reduces to the usual scale choice for inclusive photon pair production if Pjet± 
vanishes. Our prediction for the 77 + jet production cross section is intended for use in the detection 
of a Higgs boson lighter than the top quark, therefore, we assume 5 massless flavours and restrict all 
cross sections to the range of 80GeV < m^^ < 160 GeV. The electromagnetic couphng is taken at the 
Thomson limit, Oem = 1/137. We use a jet reconstruction algorithm and a set of event selection cuts, 
expected to be typical in Higgs searches. In particular, in order to find the jet, we use the midpoint 
cone algorithm [272] with a cone size of R = ^yArf + A^ = 1, with Ar] the rapidity interval and 
the azimuthal angle^^. Then, we require that both the jet and the photons have p± > 40 GeV and 
rapidity within \r]\ < 2.5. These are the same selection cuts as used in Ref. [273] for computing the 
gluon initiated 0(a|) corrections. Furthermore, we isolate both photons from the partons in a cone of 
Size I^ji^t 

22.2 Results 

In Fig. 52(a) we plot the invariant mass distribution of the photon pair. Here we see the continuum back- 
ground on which the Higgs signal is expected to manifest itself as a narrow resonance in the intermediate- 
mass range. The dotted (red) line is the leading order prediction and the solid (red) one is the differential 
cross section at NLO accuracy. The striking feature of the plot is the rather large correction. The large 
corrections are partly due to the appearance of new subprocesses at NLO as can also be read off the 
figure. The gluon-gluon scattering subprocess begins to appear only at NLO accuracy, and therefore 
it is effectively leading order. It is shown with a long dashed-dotted (magenta) line: it yields a very 
small contribution. The bulk of the cross section comes from quark-gluon scattering both at leading 
order and at NLO, shown with sparsely-dotted (blue) and long-dashed (blue) fines. The quark-quark 
scattering receives rather large corrections because the leading order subprocess can only be a quark- 
antiquark annihilation process, shown with short-dashed (green) line, while at NLO, shown with short 
dashed-dotted (green) line, there are (anti-)quark-(anti-)quark scattering subprocesses. Thus at NLO the 
parton luminosity is sizeably larger. In addition, more dynamic processes are allowed, which include 
t-channel gluon exchange. These contribute to enlarge the cross section in phase space regions which 
are disfavoured at leading order. 

A part of the large radiative corrections is accounted for by the new subprocesses; another part 
is due simply to the enlarged phase space, as can be seen from of Fig. 52(b), where the differential 
distributions in the distance Rj-y = v^|?77^-f?TJ^~+1^7^-^7p between the jet and the photons in the 
r]-(j) plane are shown, with a selection cut at Rj^, > 0.4. From the distribution for the harder photon 
we see that the radiative corrections recieve contributions from a larger part of the phase space than the 
prediction at leading order. In fact, on one hand a cut on Rj^^ > 1.5 affects the leading order and the 
NLO evaluation in the same way because they have the same shape. Thus in that case the cut does not 
reduce the correction to the m,^^ distribution. On the other hand, a cut on Rj^^.,^ > 1.5 cuts the NLO 
correction, but does not cut the leading order, so the correction is reduced. Nevertheless, the reduction 
is less then 10%, thus cutting on Rjj at Rjj > 1.5 reduces the NLO correction to the invariant mass 

^'in our NLO computation the midpoint and seedless cone algorithms yield identical cross sections. 
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Fig. 52: (a) The invariant mass distribution of tlie piioton pair at LHC energy for smooth isolation with Fty — 0.4 and e = 0.5. 
(b) The distributions in the distances between the jet and the photons in the rj-cf} plane. The photons and the jet are required 
to have transverse momentum |pxl > 40 GeV and lie in the central rapidity region of \ri\ < 2.5. The jet is reconstructed 
according to the midpoint algorithm. 



distribution of the photon pair, but marginally. In the rest of this study we require Rj^ > 1.5 and 
P-y-yi ^ 40 GeV, in addition to the same selection cuts as used in Fig. 52. 

Next we consider the dependence of the isolated cross section on the power n in the definition of 
the smooth isolation Eq. (56). We show the dependence on n for = 0.4, e = 0.5 and R-y = 0.7, 
e = 0.1 in Fig. 53(a). We see that the radiative corrections are much more than 100 % and depend on n 
strongly if the isolation cone is narrow and e is large, therefore, the fixed-order perturbative prediction 
at the NLO accuracy is not reliable in this case. On the other hand, the dependence on n is much milder 
and remains well below 100 % for Rj = 0.7, e = 0.1. Thus, the n = 0.1 line can be considered a good 
approximation to the prediction with standard photon isolation. 

In order to assess the stability of the predictions against scale variations, we show the cross section 
in a 3 GeV bin around m^-j = 120 GeV, that is the background for a hypothetical Higgs signal for a Higgs 
particle of mass 120 GeV. Fig. 53(b) shows the cross section for two sets of photon-isolation parameters. 
We show the scale variations for varying the renormalization and factorization scales separately, keeping 
the other scale fixed, as well as varying them simultaneously. The lower three curves represent the 
leading order predictions. At leading order the dependence on the renormalization and factorization 
scales is rather small, especially when the two scales are set equal (densely dotted line). Observing the 
predictions we conclude that the scale dependence at leading order does not represent the uncertainty 
of the predictions due to the unknown higher orders. The inclusion of the radiative corrections results 
mainly in the substantial increase of the cross section. For the isolation parameters R-y = 0.4 and e = 0.5 
the NLO corrections at = 1 are about 120% of the leading order prediction. In addition, the scale 
dependence is not reduced by the inclusion of the radiative corrections. If we require more stringent 
photon isolation cuts, then we find smaller corrections and a more stable prediction. For instance, in 
Fig. 53(b) we also show the scale dependence of the cross section obtained with i?^ = 1 and e = 0.1. 
We find that the cross section at NLO is about 40 % larger than the leading order prediction, and in this 
case the scale dependence at NLO is reduced as compared to the one at leading order accuracy. The 
reduction of the scale dependence when the stronger cuts (Rj = 1 with e = 0.1) are used indicates that 
the large correction obtained with the looser (default) cuts are due mainly to real emission of soft and 
coUinear partons. 
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Fig. 53: (a) Dependence of the invariant mass spectrum on the isolation parameter n. (b) Dependence on the renormalization 
and factorization scales of the cross section in a bin of 1 18.5 GeV < m^^ < 121.5 GeV at leading-order and at NLO accuracy. 



22.3 Conclusions 

We conclude by summarizing our observations concerning the QCD radiative corrections to the pp — > 
77 + jet process that constitutes part of the irreducible background to the ^ + jet — > 77 + jet 
discovery channel of an intermediate-mass Higgs boson at the LHC. We used a smooth photon isolation 
that is infrared safe to all orders in perturbation theory and independent of the photon fragmentation into 
hadrons. We found large radiative corrections, which are rather sensitive to the selection cuts and photon 
isolation parameters. Choosing a mild photon isolation, i.e. a small isolation cone radius FLy = 0.4 
with relatively large hadronic activity allowed in the cone results in more than 100 % correction with as 
large residual scale dependence at NLO as at leading order. In this case, the dependence on the isolation 
parameter n is also large. Since the radiative correction with standard photon isolation would even be 
larger, we conclude that NLO perturbation theory does not yield a reliable theoretical prediction for such 
a mild photon isolation. Making the photon isolation more stringent, for instance increasing the cone 
radius to Ry > 0.7 and decreasing the hadronic activity in the cone (using for instance e = 0.1) reduces 
both the magnitude of the radiative corrections as well as the dependence on the renormalization and 
factorization scales and that on the isolation parameter n. The more stringent isolation is also useful for 
decreasing the relative weight of the poorly known fragmentation components in the cross section, and 
the prediction obtained with the smooth isolation can be considered a good approximation (valid to the 
extent of scale ambiguities) to those with the standard isolation. Our results show that a constant factor 
is not appropriate for taking into account the radiative corrections to the irreducible background of the 
pp ^ + jet ^ 77 + jet discovery channel at the LHC. 
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23. ELECTROWEAK RADIATIVE CORRECTIONS TO HADRONIC PRECISION OBSERV- 
ABLES AT TEV ENERGIES 



23.1 Weak corrections at TeV scales 

At TeV energy scales, next-to-leading order (NLO) Electro- Weak (EW) effects produce leading cor- 
rections of the type ckew log^(s/M^), where oew = otEusin^ Ow, with oem the Electro-Magnetic 
coupling constant and 9w the Weinberg angle. In fact, for large enough s values, the centre-of-mass 
(CM) energy at parton level, such EW effects may be competitive not only with next-to-NLO (NNLO) 
(as a-Ew ~ ai) but also with NLO QCD corrections (e.g., for = 3 TeV, log^(s/M^) 10). 

These 'double logs' are due to a lack of cancellation between virtual and real W^-emission in higher 
order contributions. This is in turn a consequence of the violation of the Bloch-Nordsieck theorem in 
non-Abelian theories [274]. The problem is in principle present also in QCD. In practice, however, it 
has no observable consequences, because of the final averaging of the colour degrees of freedom of par- 
tons, forced by their confinement into colourless hadrons. This does not occur in the EW case, where 
the initial state has a non-Abelian charge, as in an initial quark doublet in proton-(anti)proton scatter- 
ing. Besides, these logarithmic corrections are finite (unlike in QCD), since Mw provides a physical 
cut-off for W^-emission. Hence, for typical experimental resolutions, softly and coUinearly emitted weak 
bosons need not be included in the production cross section and one can restrict oneself to the calcula- 
tion of weak effects originating from virtual corrections only. By doing so, similar logarithmic effects, 
^ oew log^(s/M|), are generated also by Z-boson corrections. Finally, all these purely weak contribu- 
tions can be isolated in a gauge-invariant manner from EM effects which therefore may not be included 
in the calculation. In fact, we have neglected the latter here. 

In view of all this, it becomes of crucial importance to assess the quantitative relevance of such 

weak corrections affecting, in particular, key QCD processes studied at present and future hadron col- 
liders. We show here results for the case of 6-jet-, prompt-photon and Z-production at Tevatron and 
LHC. 



23.2 Corrections to b-jet-production 

In Fig. 54 (left and right panels) we show the effects of the full ©(aloEw) contributions to the pp — 
bb{g) and pp bb{g) cross sections at Tevatron and LHC, respectively. (For details of the calculation, 
see Ref. [275].) Results are shown for the total inclusive 6-jet production rate as a function of the jet 
transverse energy. (Tree-level EW and NLO QCD effects are also given for comparison.) At Tevatron, 
O(agaEw) terms are typically negligible in the inclusive cross section, as the partonic energy available 
is too small for the mentioned logarithms to be effective. At LHC, the contribution due to such terms 
grows accordingly to the collider energy, reaching the -2% level at transverse momenta of « 800 GeV. 

Next, we study the forward-backward asymmetry at Tevatron, defined as follows: 

^ ^+\PP bb{g)] - a-[pp bbjg)] ^^^^ 
<T+ Ipp bb{g)] +a-[pp^ bb{g)] ' 

where the subscript +(— ) identifies events in which the &-jet is produced with polar angle larger(smal- 
ler) than 90 degrees respect to one of the two beam directions (hereafter, we use the proton beam as 
positive z-axis). The polar angle is defined in the CM frame of the hard partonic scattering. In the center 
plot of Fig. 54, the solid curve represents the sum of the tree-level contributions, that is, those of order 
and whereas the dashed one also includes the higher-order ones a^OEW- The effects of the 

one-loop weak corrections on this observable are rather large, indeed comparable to the effects through 
order ag [276,277]. In absolute terms, the asymmetry is of order —4% at the W, Z resonance (i.e., for 
Pt ~ Mz/2) and fractions of percent elsewhere, hence it should be measurable after the end of 

^^Contributed by: E. Maina, S. Moretti, M.R. Nolten, D.A. Ross 
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Run 2. We expect even larger effects at LHC, however, some care is here necessary in order to define the 
observable, which depends on the configuration and efficiency of the experimental apparata (so we do 
not present the corresponding plot in this instance). The ag results presented here are from Ref. [278]. 




Fig. 54: The corrections (NLO-LO)/LO due to the Qew> qIoew and Qg terms relative to the ag ones vs. the transverse 
momentum of the 6-jet for pp —> hb(g) and pp hb(g) production at Tevatron and LHC, left and right frame, respectively. 
(For LHC, we do not show the corrections due to a% terms as results are perturbatively unreliable.) In the middle frame, the 
forward-backward asymmetry vs. the transverse momentum of the 6-jet for pp — > bh{g) events at Tevatron, as obtained at 
tree-level 0(a|w) and one-loop 0(a|aEw) orders. 

23.3 Corrections to 7- and Z-production 

The neutral-current processes {V = 7, Z) 

qq gV and q{q)g q{q)V (58) 

are two of the cleanest probes of the partonic content of (anti)protons, in particular of antiquark and 
gluon densities. In order to measure the latter it is necessary to study the vector boson pr spectrum. That 
is, to compute V production in association with a jet (originated by either a quark or a gluon). We briefly 
report here on the full one-loop results for processes (58) obtained through ©(asalw)- C^^^ technical 
details of the calculation, see Ref. [279].) 

Fig. 55 shows the effects of the ©(asal-^^) corrections relatively to the 0{asaEw) Born results 
(ctEM replaces ctEw for photons), as well as the absolute magnitude of the latter, as a function of the 
transverse momentum. The corrections are found to be rather large, both at Tevatron and LHC, partic- 
ularly for Z-production. In case of the latter, such effects are of order -10% at Tevatron and -15% at 
LHC for Pt ~ 500 GeV. In general, above px ~ 100 GeV, they tend to (negatively) increase, more or 
less linearly, with px- Such effects are indeed observable at both Tevatron and LHC. For example, at 
FNAL, for EJ^-production and decay into electrons and muons with BR{Z e,fi) 6.5%, assuming 
L = 2 — 20 fb~^ as integrated luminosity, in a window of 10 GeV at px = 100 GeV, one finds 650-6500 
Z + j events through LO, hence a 5a/a ^ —1.5% EW NLO correction corresponds to 10-100 fewer 
events. At CERN, for the same production and decay channel, assuming now L = 30 fb~^, in a window 
of 40 GeV at px = 450 GeV, we expect about 1200 + j events from LO, so that aSa/a - 12% EW 
NLO correction corresponds to 140 fewer events. In Une with the normalisations seen in the top frames 
of Fig. 55 and the size of the corrections in the bottom ones, absolute rates for the photon are similar to 
those for the massive gauge boson while 0{asa'^^j^^) corrections are about a factor of two smaller. 

Altogether, these results point to the relevance of one-loop 0{asa^yif) contributions for precision 
analyses of prompt-photon and neutral Drell-Yan events at both Tevatron and LHC, also recalling that 
the residual scale dependence of the known higher order QCD corrections to processes of the type (58) 
is very small in comparison [84,280-282]. Another relevant aspect is that such higher order weak terms 
introduce parity- violating effects in hadronic observables [283], which can be observed at (polarised) 
RHIC-Spin [284]. 
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p^ (GeV) I Pt (GeV) 

Fig. 55: The LO results through 0(asQ:Ew) for the 7- and Z-production cross sections at Tevatron and LHC, as a function of 
the transverse momentum (top) as well as the size of the NLO corrections through ©(asalvv) relatively to the former. 

24. TOWARDS AUTOMATED ONE-LOOP CALCULATIONS FOR MULTI-PARTICLE PRO- 
CESSES 3^ 

24.1 Introduction 

In this decade experiments at hadron colliders explore the TeV scale. The large center of mass energies 
lead generically to multi-particle final states created by QCD initial states. The application of perturbative 
methods is justified if the scales of the problem are considerably larger than the proton mass. Then a 
systematic separation of long and short distance effects is possible and predictions for cross sections can 
be made at the loop level, which is mandatory for a reliable estimate of production rates especially at 
hadron colliders. 

The calculation of multi-particle production at the one-loop level is very challenging due to the 
combinatorial complexity of the Feynman diagrammatic approach. Although the calculation of partonic 
2^2 amplitudes at one-loop is meanwhile standard, already the number of known 2^3 1-loop 
amplitudes is very restricted. Up to now not a single Standard Model process which has generic 2 — > 4 
kinematics is computed at the one-loop level. Needless to say this is highly relevant to many search 
channels for the Higgs boson at the LHC, like gluon fusion and weak boson fusion, where additional jets 
have to be tagged to improve the signal to background ratio. For signal reactions like PP — ^ H + 0,1,2 
jets, with H — > 77, WW* ,t^t^ which are available at one-loop level, not all amplitude calculations for 
the background exist. This is due to the fact that the signal case typically contains only 5 -point functions 
at 1-loop, whereas the background has generic 2 — > 4 kinematics. As an example for such reactions 
consider PP bbbb + X, PP 77 + 2 jets + X or PP ZZ + 77 + X. These require the 
evaluation of hexagon graphs like 

The aim of our working group is to develop methods for the calculation of such hexagon amph- 

^'Contributed by: T. Binoth, J.Ph. Guillet, G. Heinrich, N. Kauer, F. Mahmoudi 
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tudes. The final goal is a tool to tackle tiie computation of general 2^4 processes at the one-loop 
level in a completely automated way. A basic ingredient of such a tool are algebraic reduction for- 
mulas. Our reduction formalism is described in the next section. As an example for the efficiency of 
our methods we discuss the 5-point 1-loop amplitude gg — > 773 in Section 24.3. The fully analytical 
treatment of hexagon one-loop amplitudes seems to be feasible analogously, if massless particles are 
considered. Examples of hexagon amplitudes in the Yukawa model calculated with our approach can be 
found in [285,286]. Other, phenomenologically relevant examples are presently under study. For the 
massive case we suggest a numerical evaluation described in Section 24.4. 



24.2 Reduction Formalism 

We will very briefly discuss the basic reduction formulas for A^-point scalar and tensor integrals. More 
details can be found in [287,288]. 



24.21 Tensor reduction 

Feynman diagrams correspond to combinations of tensor integrals. The momentum space representation 
of an A^-point tensor integral of rank Rin D = A — 2e dimensions is given by 



Here qj is a linear combination of the loop momentum k and external momenta pi, . . .p^- If N > 5 
standard reduction methods lead to a proliferation of terms with complicated denominators. This has 
to be avoided, as otherwise a stable numerical evaluation of the amplitude is hardly possible. By using 
helicity or projection methods loop momenta can always be combined with external momenta such that 
they are expressed by combinations of inverse propagators, e.g. with sj = p|, Sij = {pi + pj)"^ , 
qi = k-pi and q2 = k - pi - P2: 

2k-p2 = {qf - ml) - {qj - m|) + {si2 - ml) - (si - mj) (60) 

After expressing products of loop momenta with external vectors by propagators and canceling as many 
propagators as possible, one is left with tensor integrals which are irreducible. One finds that at most 
rank 1 A^-point integrals have to be reduced for an A^-point problem. Explicit representations for the 
irreducible tensor integrals can be found in [288]. After complete tensor reduction one is left with a 
linear combination of scalar integrals. 



24.22 Scalar reduction 

To achieve as many analytic cancellations as possible the amplitude has to be expressed by a basis of 
scalar integrals. To this end scalar A^-point integrals have to be reduced further. The scalar A^-point 
function in momentum and Feynman parameter space is given by 

1 



' = (-l)-r(Ar - D/2) I d-.i^i^^fc^ 



ITT 
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where S,, 



^3 



G 



Gii/2 



2n 



Pi 



reduction formula relates an AT-point scalar integral to {N — 1) -point scalar integrals where the 

jth propagator is omitted, and {D + 2)-dimensional remainder terms: 



TV 



-(1 



o^N, dct(G) tD+2 



Oie) 







N 



N>6 



Pj 



The basic 



(62) 



The reduction coefficients are 6j = —"^iLiSji^. One finds by iteration that each A^^-point function 
and therefore each amplitude can be expressed by 2- and 3-point functions in D dimensions and 4-point 

functions in Z?+2 dimensions only. These form the basic building blocks for an irreducible representation 
of the amplitude in terms of scalar functions. The coefficients of these scalar integrals are expected to 
simplify to a large extent. 



24.3 Analytic Calculation of gg jjg 

To give an example for our algebraic approach to A'^-point amplitudes we have reconsidered the 5-point 
1-loop amplitude gg —>■ 'yyg. While this amplitude only had been extracted indirectly from the 5-gluon 
amplitude [289] by replacing gluons by photons until recently, we present a direct calculation [290]. 

For convenience we define all particles as incoming. 

7(Pl,^l) +7(P2, A2) +5(P3,A3,C3) +5r(p4, A4,C4) +5(P5,A5,C5) (63) 

Out-States can be obtained by crossing rules. In hadronic collisions this amplitude is relevant for the 
production of photon pairs in association with a jet and as such a contribution of the background to the 
Higgs boson search chaimel 77 -I- jet. A phenomenological analysis has already been provided 

in [265,273]. The colour structure of this amplitude is simple. It can be written as 

n2„3 

p{Aj},{Cj}r ^ 0] = yC3C4C5^AlA2A3A4A5 (-^-j 

^AiA2A3A4A5 helicity dependent linear combinations of scalar integrals and a constant term which is a 
remnant of two-point functions with coefficients of order (D — 4). Six independent helicity components 
exist: -,- — -i-i-i- - -, - -i-i-i- -. As the ampUtude is finite one expects that all 

3-point functions which carry spurious infrared poles cancel. The function basis of the problem is thus 

reduced to 2-point functions = ^^''^r(^2-2e)"^^ — ^ 4-point functions in D -|- 2 dimensions 

written as [286] 

1 g 

i^jlj2 ' ^^2^3 ' ^jijs ) ~ T -^4 {Pjl ' Pj2 ' Pjs ' Pj4 + Pj5 ) (^^) 

^j4jr, ^jlj2 ^j2j-i 

and constant terms. From unitarity one expects that the +++++,++++-,-++++ amplitudes should be 
polynomial. The other helicity ampUtudes will also contain non-polynomial functions Uke logarithms 
and dilogarithms contained in I2 and F^. 

To give an example of each case we show here the results for and A only, the 

remaining ones which have also compact representations can be found in [290]. The result is expressed 
in terms of field strength tensors T'^^ = p^e^ — p^ej which satisfy the relations 

TV(^±^±) = 2pi.efpj-ef-Sijef-ef 

Pi- ^ -Pk = {sij Pk-^f - SjkPi ■ ef)/2 (66) 
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where are the polarization vectors of the gluons and photons. 

For A'^^'^^ which is polynomial we find the following result 



_ Tr(.^+^+)Tr(^3'-^4' 



„2 2 
* 12*34 



C++++-pi.jr-.p3_(3^4) 



(67) 



The coefficient is given by 



2 2 2 

^++++_ _ ^45^13^14 _ ^13^45 _^ ^45 _ ^12 + ^45 ~ ^12^45 _|_ ^13^15 _|_ ^13 ~ ^34 
S35S15S24 S15S35 Sl5'S24 S35S15 S23S45 S23 

S34S45 , Si5 - S25 S23 + S35 S23S25 , S34 + S12 , 
1 1 (OS) 

S23S15 S45 Si3 S13S45 Si5 

We have checked that the corresponding amplitude has a S2'S> S2 Bose symmetry when the photons and 
the gluons with equal helicities are interchanged. 

For A ^'^'^ we split the result into three pieces with indices F,B,l, which belong to the part 
proportional to 6-dimensional boxes Fi, a part containing bubble graphs I2, and a constant term, re- 
spectively. 



We find 



A—+++ = ^;-+++ + ^5-+++ + A^-+++ (69) 



A — +++ 



Al 



— +++ _ 



„2 2 
''12''34 



ri — +++ 



Pi- -pS- i^-^ 4) -Pl(si3, Si4, S25) 



-(4 ^ 5) - (5 ^ 3) + (1 ^ 2) - (1 ^ 2, 4 ^ 5) - (1 ^ 2, 5 ^ 3) 
Tr(.Fr.F2")Tr(^/^4+ 



„2 „2 
* 12*34 



r< — +++ 



Pl-^5^-P3-(3^4) /f (Si5 



-(4 ^ 5) - (5 ^ 3) + (1 ^ 2) - (1 ^ 2, 4 ^ 5) - (1 ^ 2, 5 ^ 3) 

2 S34S45S35 



(70) 



The indicated permutations have to be applied to the indices of the field strength tensors, momenta and 
Mandelstam variables. The coefficients are 



— +++ 



1 s2 
i *12 



2 S35.515 
S45 



2gl3.Sl4 _ £14 _ £14 
S34 S35 



Sl5 



gl3 + ^35 _^ ^14 + S45 
Si4 + S45 Si3 + S35 



*45*13 , *14.S35 „ (£15_+S45) 



Sl3 + S35 S14S45 



'45 



Sl5 
2 

^23^^5 



S15S35 S35S15 



S15S35 



_^ Sli + g24 _^ •S12 



S15S45 035 

■Sli — -^35 ^ ■Sl4(.Si,- 



■■^15) 



S45 



^35 



S35(«13 + S35] 



2^45 + Sl5 _ „ (£15 + S45)£23 _ (2^45 + ^is) ^13(2^45 + ^is) 



S13 + S35 535(513 + 535) 



S35 



'35 



In the given expressions the 5*2 S3 symmetry under exchange of the two photons and the three gluons 
is manifest after taking into account the omitted colour factor. 

The result indicates that with our approach indeed a compact representation of comphcated loop 
amplitudes can be obtained. The application of our approach to relevant 6-point amplitudes is presently 
under study. 
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24.4 Numerical Evaluation of Multi-Leg Integrals 

As already mentioned, our final aim is a complete automatisation of one-loop calculations. The bottle- 
neck for this goal is mainly given by the calculation of the virtual ampUtudes: As the number of external 
legs increases, the growing number of invariants renders the analytical expressions more and more com- 
plicated. If in addition massive particles are involved, the complexity of the resulting expressions rapidly 
approaches a limit where the analytical evaluation of the amplitude becomes unfeasible. Therefore, a 
numerical approach seems to be more appropriate to tackle different types of one-loop ampUtudes in a 
unified and efficient way. Of course, it has to be stated what "numerical" means, as any method should fi- 
nally lead to "numbers" to be compared to data. The important question is at what stage of the calculation 
the transition from analytical to numerical evaluation should be made. 

The conventional approach to calculate NLO cross sections seeks to keep analytical expressions 
in the course of the calculation of a cross section as far as possible. Of course, there are good reasons to 
do so: If infrared (and/or ultraviolet) poles are present, they have to be isolated and canceled before any 
numerical evaluation can be attempted. Further, analytic expressions are flexible in the sense that they 
can be used in "crossed" processes with different kinematics by analytic continuation. On the other hand, 
the analytical approach may be troublesome if the calculation of differential cross sections da/dO for 
some (infrared safe) observable O and/or the implementation of experimental cuts requires modifications 
of the analytic expressions. In addition, concerning the virtual integrals, it is well known that even if a 
closed form exists, the implementation into a Monte Carlo program may lead to numerical instabiUties 
because the expressions are not appropriate for every phase space region. 

These drawbacks of the "maximally analytical" approach suggest to make the transition analytical 
numerical at an earlier stage of the calculation. A completely numerical approach has been suggested 
by D. Soper [291,292], where the sum over cuts for a given graph is performed before the numerical in- 
tegration over the loop momenta, in this way exploiting unitarity to cancel soft and coUinear divergences. 
This method is very elegant, but choosing appropriate integration contours in the multi-dimensional pa- 
rameter space is far from straightforward and therefore might be hard to automate. 

In [293], a different approach has been suggested, where the calculation is split into "virtual" 
(loop) and "real" contributions as in the conventional approach, but a subtraction formalism for the UV 
and soft/colUnear divergences of the one-loop graphs has been worked out, such that the subtracted 
integrals can be performed numerically in four dimensions. While the subtractions act on a graph by 
graph basis, the subtraction terms are added back (in analytically integrated form) after having been 
summed over the graphs, as only the summation leads to expressions which are simple enough to be 
integrated analytically. 

Another promising approach in this context is the one of [294], which tries to tackle the problem 
of infrared divergences by its very root: Starting from the observation that the clash between the long- 
range nature of the interactions in a massless gauge theory and the assumption of asymptotically free 
external states causes the appearance of IR singularities in the "conventional" amplitudes, they show that 
an appropriate redefinition of the external states, which includes the long-range interactions, leads to 
S-matrix elements which are IR finite and apply it to the case e+e~ 2 jets at NLO. 

The method we suggest here to calculate one-loop amplitudes is oriented at the aim of automati- 
sation as a main guideline. It follows to a large extent the "analytical road" in treating virtual and real 
emission corrections separately. This is feasible as the isolation of IR/UV divergences is straightforward 
if the discussed reduction formaUsm is used. UV renormaUsation is well understood, and systematic 
methods for the combination of the IR divergences from the virtual corrections with their counterparts 
from the real emission contribution also exist [269, 281, 295-299]. The main problem consists in cal- 
culating the remaining finite terms of the TV-point one-loop amplitudes, especially for AT = 5,6, and 
we will concentrate on this point in what follows. By making the transition from analytical to numerical 
evaluation for these terms at an earlier stage than in the conventional approach, complicated cancellations 
between numerous dilogarithms can be avoided. 
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We should mention that another approach to calculate loop integrals numerically is suggested by 
Passarino et al. [300, 301], based on the Bemstein-Tkachov relation [302]. 



24.41 Reduction to basic building blocks 

A one-loop A^-point amplitude involving particles with arbitrary masses (including the case m = 0) 
will be reduced to basic building blocks using the method of [288], as outlined in Section 24.2 As basic 
building blocks, we choose scalar 2-point functions I2 and 3-point functions in D dimensions and 
D + 2 dimensional box functions if^"^. The latter are infrared finite. Possible UV singularities are only 
contained in the 2-point functions and their subtraction is straightforward. The (soft and coUinear) IR 
singularities are, as a result of the reduction, only contained in 2-point functions and 3-point functions 
with one or two light-like legs. In this form, they are easy to isolate and to subtract from the amplitude. 
However, the resulting expression in general still has integrable (threshold) singularities which hinder a 
successful numerical evaluation. For example, the general 6-point integral (where all internal lines have 



different masses m-i and all external legs are off-shell, pj 



1 ... 6) depends on 21 kinematic 



invariants with one non-linear constraint among them, and its analytic form contains hundreds of diloga- 
rithms. A numerical evaluation of the latter leads to large cancellations in certain kinematic regions and 
such a representation is therefore inappropriate. 



24.42 Parameter representation of basic building blocks 

After reduction and separation of the divergent parts, we are left with finite integrals and I^^'^. To 
evaluate them we first use standard Feyimian parametrisation and then perform a sector decomposition'*^ 

1 = 6(xi > X2, . . . , Xn) + Q{X2 > Xi, . . . , Xn) H 1- Q>{XN > Xi, . . . , XjV-l) (71) 

for the integration over parameters {N = 3 for the triangle, iV = 4 for the box). Now, we carry out 
only one parameter integration. We obtain 

tD t 2 2 2\ 

h {S1,S2, S3, mi, 1712,1713) = 

^Trii^^, S3, Si, ml, ml, ml) + ^^^(sa, si, S2, m|, mj, m|) -|- S^riisi, S2, S3, ml,ml, mf) 



1 

STri^{si,S2,S3,ml,ml,ml) = J dtidt2j^^-^-^ 



+ t2) Atl + Bt2 + C - id 



with 

A 
C 
R 



, 2A 
dti—= 

Vr 



log(2^ + B- Vr) - \og{B - Vr) - log(2A + B + T)+ \og{B + T) 



t + Vr 



_ log(2yl + B + VR)- log{B + VR)- log(2A + B + T) + logjB + T) 

T-^/R 

m2 , B = [mi + m2 — S2)ti + m2 + m^ — S3 
mltl + {ml + mg - Si)ti + m| 
B^-4AC + iS, T = 2A{l + ti) - B 



(72) 



We show the expUcit expressions only for the triangle, the ones for the box are analogous and can be 
found in [303]. In the case of vanishing masses or invariants, as long as the functions remain IR finite, 

analogous expressions can be derived. 



°We define the step function to be 1 if its argument is true, and else. 
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24.43 Singularity structure 

In order to analyse the singularity structure of the integrands, we explicitly separate imaginary and real 
part. One obtains 



oD=A/ 'I 'A 'A\ 

Sth (•si,S2,S3,rni,m2,m3) 



dU 



T^-R 



\og{2A + B + T)- \og{B + T) 



+e{R < 0) 
T 



log{C) - logiA + B + C) 



= ( arctan 



-R 



+e{R > 0) 

T+y/R 
2/R 



T 



arctan 



-R 



2/R 



2A + B 
log(|2^ + 5- VR\ 



lo 



+ 7r G(S <0<2A + B)^ 
log (^\B - - iTTQ{B <Vr<2A + B 

g [\2A + B + VR\^ - log [\B + VR\) + m@{B < -VR < 2A + B))] } (73) 



Three regions which lead to an imaginary part can be distinguished: 

Region I: A + B + C >0, -2A < B <0,C > 0^ {B < ±VR <2A + B). 

Region II: A + B + C > 0,C < <^ {B < VR < 2A + B) mdnot{B < -y/R < 2A + B). 

Region III: A + B + C <Q,C > {) ^ {B < -y/R <2A + B)imdnot{B < \fR <2A + B). 

Region I is an overlap region where the imaginary part has two contributions. In regions II and III only 
one of the 0-functions in (73) contributes. All critical regions are shown in Fig. 57, which illustrates the 
analytic structure of the integrand. Line segment (a) corresponds to the integration region of a triangle 




Fig. 57: Analytic regions of the box and triangle integrands. Within regions I,II and III the integrand has an imaginary part. 
The integrable square-root and logarithmic singularities are located at the boundaries of these regions. 

function where only logarithmic singularities are present, while for line segment (b) square-root and 
logarithmic singularities are present. Note that the box function I^^^ has the same singularity structure. 
As /3^=^ and /f =6 are the basic building blocks, this analysis of the singularity structure is done once 
and for all. 
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24.44 Numerical integration 

We now discuss the numerical evaluation of the analytic expressions derived above."^^ So far, we showed 
that any finite scalar A?^-point function can be written as a linear combination of the basic building block 
'S'tvT^ of Eq. (73) and a similar 2-dimensional integral, S^^^, with coefficients that are rational functions 
of the kinematical invariants and masses. To avoid numerical instabilities, the residual dimensions of 
•^TrT^ ^b7x "^^"^ integrated over numerically. Since scalar function expressions can contain 
dozens of building blocks and an amplitude has to be evaluated many times to calculate a cross section, a 
fast method to evaluate the basic building blocks is called for. However, the naive application of standard 
numerical integration techniques is not sufficient to achieve this objective. This is due to the presence 
of integrable singularities and step discontinuities in the integrands, which the detailed analysis above 
revealed. It prevents the naive application of deterministic, integration-rule based algorithms that are 
better suited for fast evaluation than the more robust, but significantly slower Monte Carlo techniques 
commonly used to evaluate multi-particle cross sections in high energy physics. 

Several approaches can be pursued to achieve a sufficiently fast and accurate numerical integration 
of the basic building blocks. A first direction are automatic methods, i.e. methods that do not require 
knowledge of the exact location and type of the discontinuities. The key to success here are adaptive 
algorithms that iteratively divide the integration volume into non-uniform subvolumes and apply basic 
numerical integration methods to each subvolume until an optimal partition of the integration volume 
minimizes the total error. Using this approach the 1 -dimensional integral of S^^'^ can be integrated 
with negligible time requirements (fractions of a second). The 2-dimensional integral of S^^^ is much 
more challenging, but can be tackled in the same spirit by combining deterministic and Monte Carlo 
integration techniques (see Ref. [304] and references therein). We note that the time required to integrate 
all S^^^ building blocks of the scalar hexagon function using this approach depends on the kinematical 
configuration and, while sufficiently short at this stage, is no longer negfigible. If, at a subsequent stage, 
i.e. for the calculation of a certain cross section, the S^^^ building blocks had to be computed in a time 
comparable to the one for the S^^^"^ functions, a second direction could be pursued. Since the location of 
all singularities and step discontinuities is known analytically, one can identify regions with continuous 
integrand and in each region flatten the integrand either by transforming integration variables or sub- 
tracting singular approximations. The resulting bounded integrands could then be integrated numerically 
with standard deterministic methods at a speed that would facihtate milhons of amphtude evaluations in 
a reasonable amount of time. 

To demonstrate the practicaUty of our method to evaluate multi-leg integrals, we show in Fig. 58 
a scan of the 2mt = 350 GeV threshold of the 4-dimensional scalar hexagon function for a kinematical 
configuration appropriate for the Feynman diagram to the right in Fig. 56. 

24.5 Conclusion 

We have outlined new algebraic/numerical approaches for 1-loop calculations. We have shown that our 
algebraic formaUsm leads to compact representations of complicated 1-loop amplitudes. Furthermore we 
have constructed numerical methods for general hexagon kinematics. The presented methods are parts 
of a project to describe multi-particle/jet production at TeV colUders with 2 — > 4 kinematics at the 1-loop 
level. 
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25. INFRARED DIVERGENCES AT NNLO 
25.1 Introduction 

Fully differential NNLO calculations are needed to increase the accuracy of theoretical predictions and 
are relevant to high-energy collider experiments, in particular for the Tevatron and the LHC. They involve 
a variety of technical complications, such as the calculation of two-loop amplitudes, a method for the 
cancellation of infrared divergences and stable and efficient numerical methods for the final computer 
program. There has been a significant progress in the calculation of two-loop amplitudes in the past 
years [305-317]. Here we review the state of the art for the cancellation of infrared divergences at 
NNLO. 

Infrared divergences occur already at next-to-leading order. At NLO real and virtual corrections 
contribute. The virtual corrections contain the loop integrals and can have, in addition to ultraviolet 
divergences, infrared divergences. If loop amplitudes are calculated in dimensional regularisation, the 
IR divergences manifest themselves as exphcit poles in the dimensional regularisation parameter e = 
2 — D/2. These poles cancel with similar poles arising from amplitudes with additional partons but less 
internal loops, when integrated over phase space regions where two (or more) partons become "close" 
to each other. In general, the Kinoshita-Lee-Nauenberg theorem [318,319] guarantees that any infrared- 
safe observable, when summed over all states degenerate according to some resolution criteria, will be 
finite. However, the cancellation occurs only after the integration over the unresolved phase space has 
been performed and prevents thus a naive Monte Carlo approach for a fully exclusive calculation. It is 
therefore necessary to cancel first analytically all infrared divergences and to use Monte Carlo methods 
only after this step has been performed. 

At NLO, general methods to circumvent this problem are known. This is possible due to the 
universality of the singular behaviour of the amplitudes in soft and coUinear limits. Examples are the 
phase-space shcing method [281,295,296] and the subtraction method [269,298,299,320,321]. It is 
worth to examine a simple NLO example in detail to understand the basic concepts which are currently 
under discussion for an extension to NNLO. We consider the NLO corrections to 7* — > 2 jets. The real 
corrections are given by the matrix element for 7* — *• qgq and read, up to colour and coupling factors 



I.43I' = 8(1 -e) 



2 2 \Xi ^ 
+ (1 - e)— + (1 - e)— - 2e 



X1X2 Xi X2 Xi X2 



(74) 



^Contributed by: G. Heinrich, S. Weinzieri 
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where xi = S12/S123 and X2 = S23/S123. This term is integrated over the three particle phase space 



4-4+3£^-5/2+2e r / ^ \ 

#3 = rsjoo^^ / (fxS I l-yxi\ x7^x:7^x7^. (75) 

r(i-e)r(|-£) y \ h J 

Singularities occur at the boundaries of the integration region at Xi =0 and X2 = 0. Historically, 

phase space slicing [281,295,296] has been the first systematic method to treat the infrared singularities. 
Here, one splits the integration region into different parts, shown in fig. 59: A soft region, given by 
Xl < Xmin and X2 < Xmin^ two coUinear regions, corresponding to x\ < Xmin,X2 > Xmin and 
Xl > Xmin,X2 < Xmin ^nd a hard region region xi > Xmin,X2 > Xmin- The hard region is free 
of singularities and the integration can be performed numerically there. In the remaining regions the 
matrix element is approximated by the soft or collinear factorisation formulae and the integration over a 
one-parton phase space can then be performed analytically. Phase space shcing has the advantage, that 
different factorisation formulae may be used in different regions of phase space. However, there are also 
some disadvantages: The method introduces a systematic error of order 3;,„j„, it becomes rather intricate 
for colour-subleading terms and it poses a numerical problem: The hard region gives a contribution of 
the form 

a It? Xmin + & In Xmin + c. (76) 

The logarithms In^ Xmin and In Xmin cancel against the contributions from the other regions, but this 
cancelation happens only numerically. 

Within the subtraction method [269,298,299,320,321] one subtracts a suitable approximation term 
da^ from the real corrections da^. This approximation term must have the same singularity structure as 
the real corrections. If in addition the approximation term is simple enough, such that it can be integrated 
analytically over a one-parton subspace, then the result can be added back to the virtual corrections da^. 

a^^O = y" da^ + y"da^ = y (da^^-da^) + J La^ + J daA . (77) 

n+l n n+1 n \ 1 / 

Since by definition da"^ has the same singular behaviour as da^, da"^ acts as a local counter-term and 
the combination {da^ — da^) is integrable and can be evaluated numerically. Secondly, the analytic 
integration of dcr^ over the one-parton subspace will yield the explicit poles in e needed to cancel the 
corresponding poles in dcr^. The subtraction method overcomes the short-comings of the slicing method, 
but there is a price to pay: The approximation term is subtracted over the complete phase space and has 
to interpolate between different singular regions. At NLO this requires an interpolation between soft and 
colhnear regions. For the example discussed above the approximation term can be taken as a sum of two 
(dipole) subtraction terms: 



da^ = \A2{p[,p's)\'' ^ 



+ \A2{p'l,P^)\ 



Sl23 
2 1 



Sl23 



2 2 , ,X2 

Xl[Xl+X2) Xl Xl 



X2{xi+X2) X2 X2 



(78) 



The momenta p'-^, P3, p'l and p'^ are linear combinations of the original momenta pi, p2 and p^. The first 
term is an approximation for xi 0, whereas the second term is an approximation for X2 0. Note that 
the soft singularity is shared between the two dipole terms and that in general the Bom amplitudes A2 are 
evaluated with different momenta. The subtraction terms can be derived by working in the axial gauge. In 
this gauge only diagrams where the emission occurs from external lines are relevant for the subtraction 
terms. Alternatively, they can be obtained from off-shell currents. Antenna factorisation [322-325] 
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S23 S23 




Fig. 59: Partitioning of the Dalitz plot for three final-state particles: (a) phase space slicing, (b) sector decomposition. 



allows to reduce the number of subtraction terms needed and interpolates smoothly between the xi ^ 
and X2 ^ regions. We briefly sketch how to obtain the antenna function from off-shell currents. The 
amplitude for 7* —>■ qgq is constructed from the quark and antiquark currents J and J as follows: 



^3(Pl,P2,P3) = J{Pl,P2) J{P3) + J{Pl) S-y-l J{P2,P3)- 



(79) 



e-y is the polarisation vector of the photon, contracted into the Dirac matrix 7^. Singular contributions 
can only arise from the currents J{pi,P2) and J{p2,P3)- With a suitable pair of reconstruction functions 
Pi{piiP2,P3) and P2{pi,P2,P3) one achieves that the Born amplitude A2 is evaluated with the same set 
of momenta and one can approximate the amplitude ^3 by 



Ant(pi , P2 , P3 ; Pi , P2) A2 {pi , P2) 



(80) 



Note that in general the momenta pi and p2 are non-linear functions of the original momenta pi, p2 and 

P3- 

Once suitable subtraction terms are found, they have to be integrated over the unresolved phase 
space. Here, one faces integrals with overlapping divergences, as one can already see from our simple 
example: 



2 2 , T2 
— + 1-e — 

Xl{Xl+X2) Xl Xl 



(81) 



The term 1 / {xi + X2) is an overlapping singularity. Sector decomposition [326-328] is a convenient tool 
to disentangle overlapping singularities. Here one splits the integration region into two sectors xi > X2 
and Xl < X2, as shown in fig. 59. In the fist sector one rescales X2 as = X2/X1, while in the second 
sector one rescales x[ = xi/x2. Sector decomposition is discussed in detail in sect. 25.3. 



25.2 The subtraction method at NNLO 

The following terms contribute at NNLO: 



(0) 
n+2 

(1) 
n+1 

dap) 



da, 



da, 



AiO) *A0) 
•^n+2 •^n+2 



On+2-, 



n + 1 '^'n-Ul \ ^^n-Ll 



^(1) ■ 



(0) 
n+1 



m * a(2) , a(2) * m , ^(1) * ^(1) 



(82) 



where An^ denotes an amplitude with n external partons and I loops, dcpn is the phase space measure for 
n partons. Taken separately, each of these contributions is divergent. Only the sum of all contributions 
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is finite. We would like to construct a numerical program for an arbitrary infrared safe observable O. 
Infrared safety implies that whenever a. n + I parton configuration pi,...,pn+i becomes kinematically 
degenerate with a n parton configuration p[,...,p'^ we must have 

On+l{Pl,-,Pn+l) On{p'i,...,Pn)- (83) 

To render the individual contributions finite, one adds and subtracts suitable pieces [329,330]: 

+ 1 (On+i dai^l^ + On+i o dai°;V - On o da^^) 

+ j (On dai?^ + On o dQi°'2) + On o da^^ 



,1) 



(84) 



Here da^^^ is a subtraction term for single unresolved configurations of Bom amplitudes. This term 
is already known from NLO calculations. The term dan'"^^ is a subtraction term for double unresolved 
configurations. Finally, don'^^ is a subtraction term for single unresolved configurations involving one- 
loop amplitudes. 

To construct these terms the universal factorisation properties of QCD amplitudes in unresolved 
limits are essential. QCD amplitudes factorise if they are decomposed into primitive amplitudes. Prim- 
itive amplitudes are defined by a fixed cyclic ordering of the QCD partons, a definite routing of the 
external fermion lines through the diagram and the particle content circulating in the loop. One-loop 
amplitudes factorise in single unresolved limits as [324,331-336] 

4') = Sing(0'i).4^1, + Sing(i'i).4°l,. (85) 

Tree amplitudes factorise in the double unresolved limits as [323, 337-342] 

= Sing(0'2).4°2,. (86) 

To discuss the term dan''^^ let us consider as an example the Born leading-colour contributions to 
e^e~ —>■ qggq, which contribute to the NNLO corrections to e+e~ —>■ 2 jets. The subtraction term 
has to match all double and single unresolved configurations. It is convenient to construct dun''^^ as a 
sum over several pieces, 

da(°.^) = Yl ^^n''\n (87) 

topologies T 

Each piece is labelled by a splitting topology. An example is shown in fig. 60. The term Vn'"^^ (T) 
corresponding to the topology shown in fig. 60 approximates singularities in l/si2, 1/ (S12S123) and part 
of the singularities in l/sf23- taken to disentangle correctly overlapping singularities like 

l/(si2S23)- Details can be found in [329]. 
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25.21 One-loop amplitudes with one unresolved parton 

Apart from dan'"^^ also the term dan'^\ which approximates one-loop amplitudes with one unresolved 
parton, is needed at NNLO. If we recall the factorisation formula (85), this requires as a new feature 
the approximation of the one-loop singular function Sing*^^'^\ The corresponding subtraction term is 
proportional to the one-loop 1 — > 2 spUtting function V^f^'^^ a-*bc' ^ example is the leading-colour part 
for the sphtting q ^ qg [330]: 



V, 



(1,1) 

(1,0) q^qg,lc,corr 



,lip(o,i) +s-^cr 



\9l,corr{y, Z) Vf4g + h — -i>e [1 - pe{l - y){l - z)]] . 

( Sijk y ) 



(88) 



This term depends on the correlations among the remaining hard partons. If only two hard partons are 
correlated, gi is given by 



gi,intr{y,z 
1 



r(i + e)r(i-e) 



1 - z 



+ 1 - (1 - vYz' 2F1 (£, e,l+e;{l- y){l - z)) 



(89) 



Here, y = sij/sijk, z = Sik/{sik + Sjfc) and /2 = (1 - pe)/2/(l - e)/(l - 2e). The parameter 
p specifies the variant of dimensional regularisation: p = 1 in the conventional or 't Hooft-Veltman 
scheme and p = in a four-dimensional scheme. For the integration of the subtraction terms over the 
unresolved phase space all occuring integrals are reduced to standard integrals of the form 



j dy y«(l - y)i+^+'' j dz z\l - zf [1 - z{l - y)Y 2F1 (e, e; 1 + e; (1 - y)z) = (90) 

r{j + e)r{j + 8)r{j + i + c) 





r(l + a)r(l + d)T{2 + a + d + e)r(l + e) 



r(2 + a + d)T{e)r{e) 



j=0 



T{j + l)T{j + 1 + e)T{j + 3 + a + c + d + e)' 



The result is proportional to hyper-geometric functions 4F3 with unit argument and can be expanded into 
a Laurent series in s with the techniques of [343, 344]. For the example discussed above one finds after 
integration [330]: 



V, 



(1,1) 



1 



31 



(1,0) q-*qg,lc:intr ^^4 I2g3 

151 55 



+ 



+ 



145 2 15 . 
6"" 24^+72"^ +T^^ 



1 



356^ 



187 

'2A 



-L 



-t-ivr 



1 



4e3 4e2 
where L = hi{sijk/ ii'^). 



29 1 



51 1 

y-4^ 


5 2 

+ 12^ 


6 y £2 






1 1663 




12 ; 


£ 16 


1 TT^L 

72 


8 


36 


7r2\ 1 

" 3 j £ " 


139 

~ ~8~ " 


11 2 



233 107 2 5 2 
21^+16"^ +12^^ 



15 



(91) 



25.3 Isolation of infrared poles by sector decomposition 

As is well known, for ultraviolet divergences a general subtraction scheme to all orders can be established 
[326, 345, 346]. For infrared poles, i.e. soft and colUnear poles in Minkowski space, the situation is less 
settled. Of course, general cancellation theorems like the KLN theorem [318, 319] exist, but a local 
subtraction scheme acting on a graph by graph basis and being valid to all orders is not available. 
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On the other hand, the use of dimensional regularisation in combination with Feynman (or alpha 
[347]) parameters allows in principle to isolate the infrared poles as powers val/e for an arbitrary graph. 
The problem is that the corresponding parameter integrals get extremely complicated the more loops 
and scales are involved, in particular they exhibit an overlapping structure. A simple example of an 
overlapping singularity already has been given in eq. (81). 

To address this problem, an automated algorithm'*^ presented in [328] and sketched below has 
been constructed which disentangles the overlapping regions in parameter space by decomposing them 
iteratively into subsectors until the divergent contributions factorise. Arrived at the factorised form, 
subtractions can be implemented easily. The resulting parameter integrals for the pole coefficients are 
in general too complicated to be integrated analytically, but as they are finite, they can be integrated 
numerically. 



25.31 Multi-loop integrals 

The straightforward automatisation of the algorithm is one of its virtues, and this is in particular true 
when appUed to virtual loop integrals, because the latter, after Feynman parametrisation, have a "standard 
form": An L-loop graph G 'vaD dimensions with propagators is, after momentum integration, of the 
form 

°r JL Tj(z:\N-{L+l)D/2 



'=1 



The functions U and F can be straightforwardly derived from the momentum representation, or they can 
be constructed from the topology of the corresponding Feynman graph [347,348]. is a polynomial of 
degree L in the Feynman parameters, T is of degree L + 1 and also depends on the kinematic invariants 
of the diagram. The sector decomposition uses representation (92) as a starting point and proceeds as 
follows: 

O The integration domain is split into N parts, using the identity 

<<oo poo ^ 

/ d^x = y2 TT ^(^i - ^j) ' (93) 

^0 Jo 

such that G becomes a sum over N integrals Gi, where in each "primary sector" I the variable xi 
is the largest one. 
O The variables are transformed in each primary sector I as follows: 





Xltj 


3<l 


Xj = < 


XI 


3 = 1 




^ Xltj-1 


3>l 



O By construction, xi factorises from U and T. We eliminate xi in each Gi using 

dxi/xi (5(1 - x^l + 51 = 1 • 



k=l 

By applying the sector decomposition iteratively, one finally arrives at a form where all singularities are 
factorised explicitly in terms of factors of Feynman parameters like t'j^~'^- Subtractions of the form 

j\tjtf-'''J^{tj,ti^j) = J^iO,ti^j) + dtjtj'-'^' !^J^itj,ti^j)-J^iO,ti^j)^ 



'''The method of sector decomposition has been used first in [326] for overlapping UV divergences, and applied in a different 
context in [327]. 
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for each j, where lim^^^o J^itj,ti^j) is finite by construction, and subsequent expansion in e leads to a 
representation of the graph G as a Laurent series in e: 



G 



2L 

E 



Ckix,{sij,mi}) 



(94) 



The pole coefficients Ck{x, {sjj, rrij}) are sums of integrals over functions of Feynman parameters. Of 
course they also contain the kinematic invariants {sij,mi} defined by the graph. For the numerical 
integration of those functions, two cases can be distinguished: If the diagrams depend only on a single 
scale, this scale can be factored out and the pole coefficients are just numbers which can be calculated 
once and for all. For diagrams depending on several scales, Uke for example the Mandelstam variables s 
and t in the case of the massless L-loop box, the kinematic invariants have to be fixed to certain values at 
which the diagram is evaluated. These have to be in the Euclidean region in order to avoid that thresholds 
spoil the numerical integration. 

As an example for a one-scale problem, we give the result for a five-loop propagator diagram, 
shown in Fig. 61. Examples of diagrams depending on four scales are given in Fig. 62, where the two 



G[5] = i-s) 

Fig. 61: A 5-loop propagator graph 



-4-5£ 



r(4 + 5e) 40.53 



straight lines flowing through the graphs denote massive propagators which can have different masses. 
These graphs are the most complicated ones occurring in the calculation of two-loop Bhabha scattering. 

5^ rxAAAAAAAA/i 



(a) 




(c) 



Fig. 62: The two-loop four point master topologies relevant for Bhabha scattering. The wavy lines are massless and the straight 
lines are massive, with external legs on-shell. The topologies from left to right are denoted by Ga, Gb, Gc- 

Only topology Ga, for the case of one single mass scale, has been calculated analytically so far by 
Smirnov [349]. A calculation of all three graphs only exists in the massless approximation [305]. Our 

2 

numerical results at two different points for Ga,h,c{s, t, u, m?,M^) = r^(l + e) ^ -Pi/e* are given in 

i=0 

Table 6. 



25.32 Phase space integrals 

As already explained in detail above, solving the problem of isolating and subtracting the infrared poles 
occurring in NNLO phase space integrals, and of integrating over the divergent subtraction terms, is a 
major step towards a (partonic) Monte Carlo program calculating processes like for example e+e^ -^2 
or 3 jets at NNLO. The method of sector decomposition can also be very useful for this task, as it is 
a general method to isolate (overlapping) poles in parameter space. Its first application to phase space 
integrals can be found in [350], further developments of the method will be sketched in the following. 

In order to be able to use the automated sector decomposition procedure for phase space integrals, 
the phase space has to be cast into a "standard form", similar to (92) for loop integrals. For example, a 
1 — > 4 parton phase space is most conveniently written as 



/ 



dPSi 



dxj 0(a 



) - \{xiXQ,X2X5,XzX4) 0(-A) , 
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{-s,-t,-u, m'^,M'^) 


(1/5,3/10,7/2,1,1) 


(5/3,4/3,5,1,3) 




Ga 


Gb 


Gc 


Ga 


Gb 


Gc 


P2 


-1.561 


-0.5255 


-1.152 


-0.08622 


-0.03483 


-0.05832 


Pi 


-5.335 


-0.2024 


-3.690 


-0.04195 


0.07556 


0.05389 


Po 


1.421 


3.606 


1.555 


0.7323 


0.1073 


0.6847 



Table 6: Results for the double box graphs for Bhabha scattering 



where the parameters Xi are rescaled Mandelstam invariants, defined by 



and Kr = (27r)4-3-Dy(D _ i)v{D - 2)ViD - 3) 2^-2^ , V{D) = /T{D /2). 

A(xiX6, X2X^^ X3X4) is the Kallen function, \{x, y, z) = x"^ + + — 2{xy + xz + yz). 

In order to calculate the real emission part where up to two particles can be unresolved using the 
method of sector decomposition, there are several ways to go, ranging from the use of this method only 
as a check for integrals over subtraction terms calculated analytically to an almost completely numerical 
approach relying largely on this method. The former has been employed in [351], where it has been 
shown that any term appearing in a phase space integral of a 1 ^ 4 matrix element in massless QCD 
can be expressed as a linear combination of only four master integrals, where one of them is the "trivial" 
integral over the phase space alone. In the notation introduced above, the remaining three master integrals 



are 



i?6 



dPSi 



,2\-4 



= (9 



dPSA 



dPSA 



{X2 + X4^ + X6)(X3 -I- a;5 -I- Xq) 

1 



X2XSX4X5 



1 



X2X3(X2 + .X4 + .T6)(a;3 -I- X5 -I- Xe) 



(95) 
(96) 
(97) 



Their analytical evaluation could be achieved by calculating iig „ explicitly and deriving the others from 
unitarity relations involving known results for three-loop two-point functions. The numerical results 
obtained by the sector decomposition algorithm. 



-^8,0 
-^8,6 

St 



Sriq"^)'^" [0.64498 + 7.0423e + A0.507e^ + 0(e' 
5.0003 0.0013 65.832 151.53 



-2e 

SriqY'-'' 

J2\-2-2e 



Sriq' 

(47r)3= 



0.74986 0.00009 14.001 



£ 

52.911 



,2^-£ 



+ 37.552 + 0{e 
-99.031 + 0(£ 



2ii;r5 r(l -e)r(2-2e) ' 
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agree with the analytical results within a numerical precision better than 1%. 

The algorithm can also be employed to avoid complicated analytical integrations over subtraction 
terms completely. The singularities can be extracted as outlined above, and the pole coefficients can be 
calculated numerically to a high precision to check their cancellation against the double virtual respec- 
tively one-loop virtual plus single-real-emission counterparts. The remaining functions are finite, and 
the combination with an arbitrary (infrared safe) measurement function is straightforward as it does not 
hamper the numerical integration. In this way, fully differential Monte Carlo programs for 1 — >^ n parti- 
cle/jet processes can be constructed. A first step in this direction already has been undertaken in [352], 
where the contribution proportional to Nf of e'^e~ — > 2, 3 or 4 jets has been calculated using sector 
decomposition techniques. 

Certainly, the automated sector decomposition algorithm applied to phase space integrals can also 
be useful in cases where massive particles are involved. As has been already proven by explicit examples 

in the case of loop integrals [353], masses do not present a principle problem for the method, but of 
course care has to be taken that thresholds do not destroy the numerical stabiUty. 
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